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Abstract 

This  study  used  advanced  laser  Doppler  velocimetry  techniques  to  measure  the  turbulence 
intensities  and  Reynolds  shear  stresses  in  Mach  2.8  turbulent  flat  plate  and  Mach  2.9  favorable 
pressure  gradient  (FPG)  boundary  layers.  The  FPG  was  generated  using  a  convex  curved  wall  and 
had  a  strength  of  /?  =  1.3 ,  where  P  is  Clauser’s  equilibrium  parameter.  The  maximum  magnitude 
of  the  “extra”  strain  rates  normalized  by  the  main  strain  rates  was  0. 1,  which  meant  the  FPG  was 
considered  to  be  a  strong  pressure  gradient.  The  flat  plate  results  indicated  that  the  LDV 
procedures  used  in  this  experiment  prevented  angular  biasing  of  the  velocity  measurements 
reported  in  the  literature.  Analysis  of  the  LDV  system  settings  also  showed  that  this  biasing,  which 
has  been  attributed  in  the  past  to  the  angular  alignment  of  the  lasers,  may  have  actually  been 
caused,  at  least  in  part,  by  the  choice  of  record  interval  used  during  data  collection.  Measurements 
in  the  FPG  test  section  demonstrated  that  the  stabilizing  effect  of  the  FPG  reduced  the  turbulence 
intensities  below  the  location  yj5  <  0.5  but  left  them  unchanged  above  that  location.  Near  the 
wall,  the  M-turbulence  intensity  was  found  to  be  reduced  to  70%  of  the  flat  plate  value.  In  addition, 
the  FPG  reduced  the  magnitude  of  the  incompressible  Reynolds  shear  stresses  {-'pii'v')  by 
approximately  75%.  Due  to  the  reduction  in  the  turbulent  shear  stress,  the  wall  shear  stress  was 
found  to  be  reduced  by  37%.  Comparison  of  the  LDV  data  to  hot-wire  data  collected  in  the  same 
facilities  showed  that  the  assumption  of  />' «  0 ,  used  in  the  reduction  of  the  hot-wire  data,  was 
valid  in  the  flat  plate  region  but  not  the  FPG  region.  The  increase  in  the  magnitude  of  p'  was 
likely  due  to  the  streamline  curvature  associated  with  the  generation  of  the  FPG  and  the  resulting 
pressure  difference  across  the  boundary  layer. 


EXPERIMENTAL  INVESTIGATION  OF  A 


SUPERSONIC  TURBULENT  BOUNDARY  LAYER 

INCLUDING 

FAVORABLE  PRESSURE  GRADIENT  EFFECTS 


1.  Introduction 


1.1  Motivation 

In  the  fields  of  aerospace  engineering  and  fluid  mechanics,  conceptual  design  projects  are 
relying  more  and  more  on  computational  fluid  dynamics  (CFD)  simulations  to  provide  required 
data  and  test  results  as  the  cost  of  wind  tunnel  and  laboratory  testing  continues  to  increase.  In 
order  to  ensure  the  validity  of  these  designs  and  the  data  acquired  from  the  numerical  simulations, 
the  computer  codes  need  to  accurately  reflect  the  physics  of  the  problems  they  are  attempting  to 
solve. 

In  problems  where  turbulence  is  involved,  today’s  computers  do  not  have  sufficient  speed 
or  memory  to  resolve  the  length  and  time  scales  associated  with  turbulence.  For  instance,  a  Mach 
2.9  flow  with  a  freestream  Reynolds  number  of  Re^  « 1.6  x  10®  and  a  boundary  layer  thickness  of 

5  «  0.009  m  has  a  Kolmogoroff  turbulent  length  scale  of  approximately  2x10”’  m.  [6]  In  a 
control  volume  6.35  cm  x  6.35  cm  x  80  cm  (the  size  of  the  test  section  in  this  experiment)  there 
would  have  to  be  approximately  4  x  lO”  grid  points  in  order  to  accurately  calculate  a  solution  to 
the  problem.  This  would  require  approximately  4xl0’*  calculations  to  obtain  a  numerical 
solution  of  the  Navier-Stokes  equations  [5];  assuming  an  optimistic  rate  of  one  compulation  per 
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nanosecond,  the  computations  would  take  12  x  lO'^  years  to  complete!  Trying  to  compute  the 
solution  for  flow  over  an  object  of  engineering  interest,  such  as  an  airplane,  would  take  even 
longer.  This  being  the  case,  engineers  and  scientists  must  rely  on  an  approximate,  time-averaged 
form  of  the  governing  equations  in  order  to  generate  solutions  to  problems  of  interest.  Due  to  the 
non-linear  nature  of  the  Navier-Stokes  equations,  the  time-averaging  process  introduces  cross¬ 
correlation  fluctuation  terms  into  the  governing  equations  (see  Chapter  2).  This  leads  to  the 
problem  of  turbulent  closure.  In  short,  there  are  more  unknowns  than  equations,  and  the  system 
cannot  be  solved.  Turbulence  modeling  attempts  to  relate  these  extra  variables  to  other  flow 
quantities  in  order  to  close  the  system. 

1. 2  Current  Turbulence  Models 

Turbulence  model  development  has  taken  the  same  path  as  most  engineering  problems: 
solutions  for  simplified  cases  were  investigated  first  and  then,  with  time,  solutions  for  more 
complicated  problems  were  developed.  The  first  turbulence  models  were  developed  for  low-speed, 
incompressible  flow  over  a  flat  plate;  these  turbulence  models  were  found  to  work  reasonably  well 
for  mild  pressure  gradient  flows,  and  thus  their  use  for  calculating  the  effects  of  pressure  gradient 
was  generally  accepted.  The  next  step  in  the  evolution  of  turbulence  modeling  was  their  extension 
to  compressible  flows. 

Compressible  (and  especially  supersonic)  flows  have  different  physical  characteristics  than 
incompressible  flows,  but  incompressible  turbulence  models  are  used  in  compressible  calculations 
on  an  ad-hoc  basis.  The  incompressible  models  provide  fairly  accurate  results  when  computing 
flat  plate  flows,  since  most  of  the  differences  between  compressible  and  incompressible  flows  can 
be  accounted  for  by  differences  in  density  and  temperature  which  occur  across  the  boundary  layer 
[32,  33].  However,  in  supersonic  flows  pressure  gradient  effects  can  cause  large  changes  in  the 
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turbulence  quantities.  In  this  case,  the  ad-hoc  extension  of  subsonic  models  to  supersonic  flow  is 
not  valid,  and  new  turbulenee  models  need  to  be  developed. 

Turbulence  models  take  on  many  forms,  but  the  one  thing  all  turbulence  models  have  in 
common  is  that  they  presently  need  to  be  validated  against  empirical  data  gathered  from 
experimentation.  As  a  result,  it  becomes  necessary  to  have  a  broad  database  of  experimental 
results  for  generating  and  validating  turbulence  models.  A  broad  database  of  turbulence 
information  is  also  needed  in  order  to  determine  what  physical  processes  control  the  turbulence; 
hopefully,  knowledge  of  how  the  turbulence  changes  in  different  physical  environments  will  allow 
the  extension  of  computer  models  to  new  situations  without  requiring  expensive  testing  to  validate 
the  results. 

i.  3  Status  of  Current  Turbulence  Database 

Although  there  have  been  many  experiments  performed  in  the  past  to  measure  quantities  in 
compressible  flows,  the  available  turbulence  database  is  still  fairly  scarce.  Many  of  the  past 
experiments  have  focused  on  determining  flow  characteristics  over  a  flat  plate,  and  many  of  the 
experiments  attempting  to  measure  turbulence  quantities  have  been  found  to  be  flawed.  According 
to  Spina  et  al.  [33], 

The  reason  for  the  scarcity  of  measurements  and  their  generally  poor 
quality  is  simple:  the  measurement  of  turbulence  quantities  in  supersonic 
boundary  layers  is  exceedingly  difficult,  with  the  level  of  difficulty  increasing  with 
flow  complexity  and  Mach  number. 

In  addition  to  the  problems  that  exist  with  data  acquisition  in  supersonic  flows,  the  physical 
characteristics  of  supersonic  flows  provide  several  means  for  generating  a  favorable  pressure 
gradient  (FPG).  The  FPG  can  be  generated  by  imposing  an  expansion  shoek  wave  onto  a  flat  plate 
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boundary  layer  (see  Lewis  et  al.  [25]),  expanding  the  flow  suddenly  around  a  comer,  or  expanding 
the  flow  slowly  over  a  gradual  curve.  The  effect  of  the  FPG  on  the  turbulence  will  vary  depending 
on  the  method  with  which  the  FPG  is  generated.  [32]  Spina  et  al.  [33]  noted  that 

The  behavior  of  supersonic  boundary  layers  subject  to  convex  curvature 
has  not  been  studied  extensively... No  data  have  been  taken  on  smooth  convex 
comers  in  supersonic  flow  leaving  the  expansion-comer  studies. .  .as  the  only  flows 
that  provide  insight  into  the  effects  of  convex  streamline  curvature... However, 

[these]  studies  actually  provide  more  insight  into  the  effect  of  a  rapid  application 
of  bulk  dilatation  than  the  effect  of  streamline  curvature. 

In  a  review  of  the  literature,  only  one  study  was  found  which  attempted  to  analyze  the 
effects  of  a  FPG  generated  by  a  convex  curved  wall.  However,  this  study  (see  Ref  [19])  was  of 
limited  use  because  the  FPG  was  immediately  preceded  by  an  adverse  pressure  gradient  (APG) 
region.  The  APG  caused  the  turbulence  quantities  to  be  disturbed  from  their  flat  plate  values, 
limiting  the  effectiveness  of  the  data.  [32] 

1. 4  Overview  of  Current  Experiment 

The  current  experiment  was  intended  to  expand  the  available  turbulence  database  by 
providing  much-needed  data  regarding  the  effects  of  streamline  curvature  and  FPG  effects.  The 
data  were  collected  using  advanced  laser  Doppler  velocimetry  techniques  to  measure  the  flow 
veloeities  and  turbulent  quantities  over  a  flat  plate  and  a  convex  curved  surface  (see  Figure  1.1  for 
a  sketch  of  the  convex  surface).  In  a  supersonic  flow,  the  convex  curved  surface  resulted  in  the 
generation  of  a  FPG.  Comparison  of  the  flat  plate  data  and  the  FPG  data  allowed  the  effects  of  the 
FPG  to  be  examined.  In  addition,  comparison  of  the  LDV  data  with  hot-wire  data  collected  by 
Miller  [26]  in  the  same  facilities  allowed  for  examination  of  the  assumptions  used  in  the  hot-wire 
data  reduction.  This  examination  was  possible  due  to  the  advantages  offered  by  the  LDV  system. 
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Figure  1.1:  Sketch  of  FPG  Test  Section 


1.5  Advantages  of  LDV 

LDV  has  been  in  common  use  as  a  flow  measurement  device  since  the  early  1970’s.  It 
offers  the  advantage  of  a  “non-intrusive”  means  of  measuring  flow  velocities  and  turbulence 
quantities.  The  use  of  LDV  systems  to  perform  data  acquisition  is  considered  non-intrusive 
because  no  flow-obstructive  measurement  device  needs  to  be  placed  inside  the  test  section; 
however  (as  will  be  discussed  in  Chapter  3),  the  flow  needs  to  be  seeded  with  particles  which 
reflect  the  laser  light,  and  careless  seeding  can  cause  the  flow  to  be  disturbed  from  its  original 
state. 

Aside  from  providing  a  non-intrusive  method  of  measurement,  LDV  has  the  advantage  of 
measuring  the  velocities  and  velocity  fluctuations  directly.  The  hot-wire  measurements  collected 
by  Miller  and  other  researchers  actually  measure  density-weighted  velocities;  in  order  to  extract 
the  pure  velocities  fi’om  the  collected  data,  several  assumptions  must  be  made.  Comparison 
between  hot-wire  and  LDV  data  has  proven  these  assumptions  valid  for  the  flat  plate  case,  but  the 
lack  of  turbulence  data  in  FPG  regions  has  not  allowed  the  assumptions  made  in  hot-wire  data 
reduction  to  be  validated  for  this  case.  One  of  the  goals  of  the  current  experiment  was  to  compare 
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hot-wire  data  to  LDV  data  in  order  to  verify  or  reject  the  assumptions  made  in  the  case  of  a 


favorable  pressure  gradient. 


1.6  Objectives 

The  primary  objective  of  this  research  was  to  obtain  turbulence  measurements  and  gain 
insight  into  the  important  physical  mechanisms  in  a  turbulent  boundary  layer  subjected  to  a 
favorable  pressure  gradient.  In  addition,  the  LDV  data  was  to  be  compared  to  hot-wire  data  taken 
in  the  same  experimental  configuration  to  determine  if  the  assumptions  used  in  reducing  the  hot¬ 
wire  data  were  valid.  Finally,  this  study  was  intended  to  increase  the  amount  of  reliable  turbulence 
information  available.  In  order  to  make  the  presented  results  useful  to  other  researchers,  the 
criteria  established  by  Settles  and  Dodson  [31]  was  used  to  provide  a  minimum  standard  by 
which  the  results  were  presented; 

1.  Baseline  Applicability:  All  candidate  studies  for  use  must  be  experiments 
involving  turbulent  flows  in  either  supersonic  or  hypersonic  Mach  number 
range  (i.e.,  M  «  3.0  or  higher). 

2.  Simplicity:  Experimental  geometries  must  be  sufficiently  simple  that  they  may 
be  modeled  by  CFD  methods  “without  enormous  difficulty.” 

3.  Specijic  Applicability:  All  candidate  studies  passing  this  criterion  must  be 
capable  of  providing  some  useful  test  of  turbulence  modeling. 

4.  Well-defined  experimental  boundary  conditions:  All  incoming  conditions 
(especially  the  state  of  the  incoming  boundary  layer)  must  be  carefully 
documented.  For  studies  claiming  “two-dimensional”  flow,  data  indicating  the 
extent  of  the  spanwise  flow  variations  should  be  provided. 

5.  Well-defined  experimental  error  bounds:  The  experimentor  must  provide  an 
analysis  of  the  accuracy  and  repeatability  of  the  data,  or  error  bars  on  the  data 
themselves.  Further,  error  bounds  on  the  data  must  be  substantiated  in  a 
quantifiable  manner. 

6.  Adequate  documentation  of  data:  Data  must  be  documented  and  tabulated  in 
a  machine-readable  form. 
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7.  Adequate  spatial  resolution  of  data:  Experiments  must  present  data  of 
sufficiently  high  resolution,  compared  with  the  scaled  flow  in  question,  such 
that  the  key  features  of  the  flow  are  clearly  resolved. 

In  addition  to  the  criteria  established  by  Settles  and  Dodson,  an  attempt  was  made  to 

provide  adequate  data  for  converting  reported  values.  That  is,  other  researchers  may  wish  to  scale 

the  data  in  a  different  manner  than  chosen  within  this  report;  an  effort  was  made  to  provide  the 

data  required  to  convert  the  reported  values  into  any  common  non-dimensional  form. 

1. 7  Synopsis  and  Methodology  of  Current  Research 

In  order  to  determine  what  measurements  were  required  in  the  development  of  a  turbulence 
model,  it  was  first  necessary  to  examine  the  governing  equations  of  fluid  motion;  this  examination 
is  carried  out  in  Chapter  2.  The  analysis  of  the  governing  equations  showed  the  necessity  of 
measuring  the  velocity  cross-correlation  term  u'v' .  The  cross-correlation  term  could  be  measured 
directly  using  an  LDV  system;  in  order  to  gain  an  understanding  of  the  LDV  system,  the  basic 
principles  behind  LDV  were  studied  and  are  presented  briefly  in  Chapter  3.  The  actual  LDV 
system,  along  with  all  other  equipment  used  in  the  research,  is  described  in  Chapter  4.  Chapter  5 
details  how  this  equipment  was  put  to  use  in  order  to  obtain  velocity  measurements  in  the  flat  plate 
and  FPG  flowfields.  Chapter  6  describes  how  these  measurements  were  manipulated  to  calculate 
other  variables  of  interest,  while  Chapter  7  presents  and  analyzes  the  results  from  the 
measurements  and  manipulations.  Finally,  Chapter  8  summarizes  the  findings,  draws  conclusions 
from  the  analysis  of  the  data,  and  presents  recommendations  for  future  research. 
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2.  Development  of  Governing  Equations 


When  studying  turbulent  compressible  flows,  it  is  necessary  to  provide  a  detailed 
development  of  the  governing  equations  for  two  reasons.  First,  the  development  provides  insight 
into  which  terms  in  the  governing  equations  contain  the  turbulent  characteristics  of  the  flow,  as 
well  as  showing  how  turbulent  compressible  flows  differ  from  incompressible  ones.  Second  (and 
probably  most  important),  the  development  shows  how  and  why  certain  variables  are  defined 
differently  for  compressible  and  incompressible  flows.  These  definitions  are  important  for 
comparing  data  collected  by  different  researchers,  as  the  definitions  used  throughout  the  literature 
are  not  always  consistent. 


2. 1  Compressible  Navier-Stokes  Equations 

The  governing  equations  for  fluid  flows  are  the  Navier-Stokes  equations.  They  are  given 
in  compressible  form  [1]  as 
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2. 2  Time  Averaging  of  the  Governing  Equations 

When  turbulence  is  present  in  fluid  flows,  variables  in  the  governing  equations  fluctuate 
very  rapidly,  leading  to  minute  time  scales  associated  with  the  flow.  As  the  time  scale  associated 
with  a  flow  decreases,  the  amount  of  computational  work  required  to  calculate  a  solution  to  the 
governing  equations  increases.  Therefore,  instead  of  trying  to  capture  the  minute  time  scales  using 
the  unsteady  form  of  the  governing  equations,  the  equations  are  simplified  using  some  sort  of  time¬ 
averaging  (Reynolds  or  Favre  averaging)  process.  In  this  way,  only  the  time  dependence  of  the 
mean  flow  needs  to  be  modeled  with  the  time-dependent  terms  of  the  Navier-Stokes  equations,  and 
the  associated  time  scales  can  be  dramatically  increased.  This  increase  in  the  length  of  the  time 
scales  allows  the  computation  time  to  be  decreased  to  a  reasonable  level. 

In  any  of  the  various  types  of  time  averaging,  an  instantaneous  variable  g  is  assumed  to 
consist  of  a  time  averaged  part,  g,  and  a  fluctuation  from  the  mean,  g' : 

g(x,y,z,t)  =  g(x,y,z,t)  +  g’(x,y,z,t)  (2.3) 

It  should  be  noted  that,  in  the  above  equation,  g  is  a  time  average  but  is  still  allowed  to  vary  with 
time.  This  allows  the  averaging  process  to  capture  possible  time-dependence  of  the  mean  flow. 
However,  it  requires  that  the  period  over  which  quantities  are  averaged,  7] ,  is  long  compared  to 
the  turbulent  time  scale  but  short  compared  to  the  overall  time  scale,  T2 .  [37]  See  Figure  2. 1  for 
a  graphical  representation  of  this  concept. 

As  mentioned  above,  there  are  traditionally  two  types  of  time-averaging  performed  on  the 
Navier-Stokes  equations:  Reynolds  averaging  and  Favre  averaging.  The  set  of  equations  resulting 
from  Reynolds  averaging  is  referred  to  as  the  Reynolds-Averaged  Navier-Stokes  Equations 
(RANS);  similarly,  use  of  Favre'  averaging  leads  to  the  Favre'-Averaged  Navier-Stokes  Equations, 
or  FANS.  Although  a  description  of  Reynolds  averaging  would  be  sufficient  for  explaining  the 
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measurements  made  by  LDV  systems,  a  description  of  the  Favre  averaging  process  is  also 
presented  in  order  to  compare  the  LDV  data  with  hot-wire  and  CFD  results. 

2.2.1  Reynolds  Averaging.  Reynolds  averaging  consists  of  replacing  all  turbulent 
variables  in  the  Navier-Stokes  equations  with  a  time-average  plus  fluctuation  term,  and  then  taking 
the  time  average  of  the  whole  equation.  Assuming  steady  mean  flow  and  no  body  forces,  the 
RANS  equations  become  [1,6] 

dt  dXj 

dj^Uj  +p'u'^  +r,^) 

dt  dXj  dx^  dxj 

d(peo  +p'h^)  dj^pUj)  _  +  ^/r'-  -  qj  -  q] ) 

dt  dXj  dxj 

where 
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3  * 


Tjj  =  -pujUj  -  UjP'Uj  -  Ujp'u;  -  p'ulUj 


q,  =  -k 


ffT 

dXi 


qj  =  phl^ul  +  h^pX  +  p'hl,  -  p%ul 
The  quantities  and  are  the  stagnation  conditions. 


(2.5) 


2.2.2  Favre  Averaging.  Favre  averaging  is  similar  to  Reynolds  averaging,  except  that 
the  mass-dependent  terms  in  the  governing  equations  are  replaced  by  a  mass-weighted  time 
average,  (p ,  plus  a  fluctuation,  (p"  : 

(p  =  p-\-<p"  (2.6) 

where 

P 

With  this  definition,  the  FANS  equations  are  given  as  [1] 

Pt  dx. 

^  Pp 

dt  Pxj  dx^  Pxj 

^pgp)  ^  ^  P^u.f..  +  M/y..  -  q.  -  qj ) 

Pt  Pxj  dXj 

In  the  Favre'  -averaged  equations,  the  definitions  of  ?,y ,  r,y  and  are  somewhat  different: 
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■ij 
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qj  ='^';K. 


2. 3  Comments  on  the  RANS,  FANS,  and  Variable  Definitions 

The  term  r,y  as  given  in  the  RANS  is  known  as  the  Reynolds  shear  stress;  the  differences 

which  arise  in  this  term  due  to  the  effects  of  compressibility,  as  well  as  the  differences  between  the 
definitions  of  t]j  in  the  RANS  and  FANS  result  in  inconsistent  nomenclature  throughout  the 

literature.  For  instance,  the  term  -fiu'v'  is  the  only  term  appearing  in  the  Reynolds-averaged 
form  of  when  the  flow  is  incompressible.  As  a  result,  many  researchers  refer  to  this  term  alone 
as  the  Reynolds  shear  stress,  neglecting  the  other  three  terms  in  the  definition  of  Ty  .  In  order  to 

avoid  confusion  within  this  report,  the  term  -Jni'V  will  be  referred  to  as  the  incompressible 
Reynolds  shear  stress. 

Another  source  of  confusion  is  the  use  of  r,y  to  represent  the  turbulent  shear  stress  terms 

in  the  FANS;  many  researchers  consider  this  term  as  the  Reynolds  shear,  even  though  it  is  not 
equivalent  to  either  the  incompressible  Reynolds  shear  or  the  entire  Reynolds  shear  stress.  The 
differences  between  r,y  in  the  RANS  and  FANS  can  be  derived  by  examining  the  differences 

between  the  definitions  used  in  the  time-averaging  process.  In  both  cases,  the  instantaneous 
velocity  is  replaced  by  a  mean  plus  a  fluctuating  term; 

«  =  u  =  u+u"  (2.10) 

Setting  the  two  equations  equal  to  each  other  and  solving  for  u'  leads  to 

u'  =  u-u  +  u"  (2.11) 
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Now,  taking  the  time-average  of  the  Favre  expansion  (the  second  expression  in  Eq.  (2.10)), 

u=u+u"  =>  u-u=-u"  (212) 

Insertion  of  Eq.  (2.12)  into  Eq.  (2. 1 1)  results  in  the  expression 

u'  =  u"-u"  (213) 

A  similar  relation  can  be  developed  for  the  v-component  of  velocity;  multipl5dng  the  two 
expressions  together  and  taking  the  time-average  leads  to 

u'v'  =  u"v"  —  —  v"u"  +  u"v"  ~  u"v"  —  u"v"  (214) 

Note  that  the  time-average  of  a  Favre  fluctuation  term,  ,  is  not  equal  to  zero.  In  fact. 


this  time  average  is  related  to  the  Reynolds-averaged  variables  by  the  relation 


P  P 


(2.15) 


p  p  p 


P 

Multiplying  Eq.  (2.14)  by  the  mean  density,  the  incompressible  Reynolds  stress  is  related 
to  the  Favre-averaged  stress  by 

'pu'v'  =  ~pu"v"  -'pu"v"  (2- 16) 


The  term  pu"v"  is  the  term  usually  reported  as  the  Reynolds  stress  by  researchers  who  prefer  to 

use  the  Favre-averaged  equations.  However,  this  term  differs  from  the  incompressible  Reynolds 
stress  by  the  amount 

(2.17) 


P'W  P’v' 


pu"v"  -  p- 

p  p 

Since  this  term  is  a  fourth-order  fluctuation,  it  is  usually  considered  negligible  so  that 


pu'v'  =  pu"v" 


(2.18) 
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2.4  Hot-Wire  and  CFD  Results  vs.  LDV Measurements 


When  comparing  hot-wire  and  CFD  measurements  to  data  collected  by  LDV,  it  is 
necessary  to  understand  the  differences  between  the  measurements,  the  assumptions  used  to  reduce 
the  data,  and  the  possible  sources  of  inconsistency  between  the  results. 

2.4.1  Hot-Wire.  Hot-wire  techniques  measure  the  change  in  voltage  required  to  keep  a 
wire  at  constant  temperature;  this  change  in  voltage  is  sensitive  to  the  heat  transfer  from  the  wire, 
which  is,  in  turn,  sensitive  to  a  mass-weighted  velocity.  Since  the  heat  transfer  to  the  wire  is 
related  to  the  flow  Reynolds  number  ,  the  hot-wire  process  is  sensitive  to  the  mass  flux  across  the 
wire,  pW .  [8]  In  order  to  separate  the  density  from  the  velocity  components,  the  assumption  that 
p'  =  0  must  be  made.  [7]  This  assumption  is  somewhat  controversial,  and  comparison  between 
hot-wire  results  and  LDV  data  should  hopefully  prove  or  disprove  its  validity. 

2.4.2  CFD.  CFD  attempts  to  model  ;  in  order  to  reduce  the  number  of  unknowns  in 
the  model,  most  CFD  codes  use  Favre'-averaging.  Comparison  of  tJj  in  equations  (2.5)  and  (2.9) 
shows  the  use  of  Favre'-averaging  reduces  the  number  of  unknowns  in  from  four  to  one. 

However,  LDV  measures  the  product  u'v'  in  the  incompressible  Reynolds  shear,  not  the  Favre- 
averaged  shear.  As  shown  in  section  2.3,  the  difference  between  the  Reynolds-averaged  results 
obtained  by  LDV  and  the  Favre'-averaged  results  from  CFD  is  usually  considered  negligible; 
however,  in  cases  where  the  fluctuating  velocity  is  large,  discrepancies  between  the  two  results  may 
arise. 
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2.5  Comments  on  Coordinate  System  Labeling 

Throughout  this  report,  four  different  coordinate  systems  will  be  used.  In  order  to  preempt 
any  confusion,  the  four  coordinate  systems,  the  labeling  scheme  for  each,  and  the  differences 
between  them  are  presented  here. 

2.5.1  Wind  Tunnel  Coordinate  System.  The  wind  tunnel  coordinate  system,  denoted  by 

{x^t  ^  Cartesian  system  with  the  tunnel  nozzle  at  =  0 .  The  y-coordinate  is  pointed 

in  the  vertical  direction  and  the  z-coordinate  is  directed  accordingly  (see  Figure  2.2). 

2.5.2  Test  Section  Coordinate  System.  The  test  section  coordinate  system  is  labeled 
and  is  used  in  the  computation  of  the  wall  coordinates.  The  origin  of  the  test  section 

system  is  placed  at  =  0.60  m ;  the  system  orientation  is  shown  in  Figure  2.2. 


Ywt 


2.5.3  Traverse  Coordinate  System.  The  traverse  coordinate  system  dictates 

the  coordinates  used  to  control  motion  of  the  traverse.  Its  origin  is  placed  on  the  wall  of  the  test 
section  with  the  orientation  shown  in  Figure  2.3. 
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2.5.4  Body-Intrinsic  Coordinate  System.  The  final  coordinate  system  is  the  body- 
intrinsic  system,  labeled  (x,y,zj .  The  body-intrinsic  system  and  the  traverse  system  have  the  same 
origin  but  different  orientations;  its  orientation  is  also  shown  in  Figure  2.3. 


Zj  z 


Figure  2.3:  Traverse  and  Body-Intrinsic  Coordinate  Systems 

Finally,  it  should  be  noted  that  the  definitions  given  in  section  2.2  are  for  a  Cartesian 
coordinate  system.  However,  it  can  be  shown  that,  for  two-dimensional  flows  where  d  «/?(/?  is 
the  radius  of  curvature  of  the  wall),  the  definitions  given  in  section  2.2  are  still  valid  in  a  body- 
intrinsic  coordinate  system.  [4] 
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3.  Principles  of  Laser  Doppler  Velocimetry 


Laser  Doppler  velocimetry  (LDV)  provides  a  “non-intrusive”  means  of  measuring  flow 
velocities  and  turbulence  levels.  There  are  many  factors  which  influence  the  measurements  taken 
when  using  an  LDV  system;  in  order  to  understand  the  potential  bias  errors,  it  is  necessary  to 
provide  a  description  of  how  the  LDV  system  works. 


3. 1  Laser  Interference  in  the  Measurement  Control  Volume 

When  two  laser  beams  of  like  wavelength  intersect,  the  wave-like  properties  of  the  light 
interact  to  create  an  interference  pattern.  If  the  intersection  occurs  at  the  waist  of  the  two  beams, 
the  interference  pattern  is  a  set  of  fringes  which  are  nearly  parallel  [11].  The  spacing  between  the 
fringes  is  a  function  of  the  laser  beam  intersection  angle,  6,  and  the  wavelength  of  the  laser  beams, 
X,  and  is  given  as  [10,  20] 


Xjl 

^  sin(0/2) 


(3.1) 


When  a  particle  passes  through  the  area  of  intersection,  it  will  scatter  light  from  the 
fringes,  creating  an  alternating  “on-ofF’  signal  as  it  travels  through  the  light  and  dark  portions  of 
the  fringe  pattern.  The  set  of  on-off  signals  associated  with  a  single  particle  is  referred  to  as  a 
Doppler  burst.  A  typical  Doppler  Burst  is  shown  in  Figure  3.1.  As  the  particle  crosses  a  fringe, 
the  intensity  of  scattered  light  increases,  resulting  in  a  peak  in  the  Doppler  burst;  as  it  passes 
through  the  dark  space  between  the  fringes,  a  valley  is  created  in  the  Doppler  burst  pattern. 

The  intensity  distribution  across  a  laser  beam  is  not  uniform,  resulting  in  the  large  intensity 
near  the  center  of  the  Doppler  burst.  The  non-uniform  intensity  distribution  also  confines  the 
actual  measurement  volume  to  the  space  where  the  scattered  light  would  be  strong  enough  to  be 
detected;  this  results  in  the  control  volume  and  fringe  bands  being  restrieted  to  an  ellipsoid  in  the 
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Particle  Crosses 
Fringe 


center  of  the  intersection  region  and  fringe  disks  within  the  ellipsoid,  respectively  (see  Figure  3.2). 

[11] 

The  Doppler  burst  signals  are  collected  by  photomultipliers  (also  known  as  PM  tubes)  and 
processed  by  a  Burst  Spectrum  Analyzer  (BSA).  The  BSA  processes  the  received  Doppler  burst 
patterns  in  order  to  determine  the  velocities  of  the  particles  traveling  through  the  control  volume. 


Figure  3.2:  Fringe  Disks  and  Measurement  Volume 
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3.2  Velocity  Calculation 


The  time  taken  to  travel  between  two  successive  fringes  corresponds  to  the  time  between 
two  successive  peaks  in  the  Doppler  burst  and  is  labeled  (see  Figure  3.3);  the  time  has  a 
corresponding  frequency,  /q.  The  relationship  between  Tj^,  and  the  velocity  is  [10] 


U  =-^=  JnXf  =  fn - 7 - T 

-/D  /  Jd 


(3.2) 


The  ratio  [/l/2]/[sin(0/2)]  is  termed  the  calibration  factor  and  is  denoted  by  the  symbol  . 


The  calibration  factor  is  dependent  on  the  optical  setup  alone  and  is  known  before  any 
measurements  are  made;  therefore,  in  order  to  calculate  the  velocity  it  is  only  necessary  to 
determine  the  Doppler  frequency,  . 

The  Doppler  frequency  is  determined  using  N  samples  (referred  to  as  the  record  length)  to 
model  the  Doppler  burst.  The  first  sample  is  taken  when  the  intensity  of  the  Doppler  signal  is 
larger  than  a  certain  threshold  (25  mV  in  Figure  3.3).  Additional  samples  are  then  taken  every  Tg 


BSA 

Triggered 


25  mV 


I  I  I  I  I  I  I 

X 


N  Samples  taken 
timeTs  apart 


Figure  3.3:  Doppler  and  Sampling  Times 
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seconds.  If  the  Doppler  burst  signal  passes  a  second  threshold  (50  mV)  the  count  of  N  samples  is 
restarted.  That  is,  the  samples  already  collected  are  deleted  and  the  sample  counter  is  reset  to  one. 
The  sample  count  can  be  reset  at  two  more  intensity  levels,  100  and  150  mV.  This  procedure 
ensures  the  part  of  the  Doppler  burst  with  the  highest  intensity  is  used  for  the  frequency 
calculation.  [10]  If  the  particle  passes  through  each  fringe  at  a  given  time  interval,  the 
Doppler  burst  will  be  a  periodic  signal  with  period  Tj^ .  Fourier  theory  shows  that  the  frequency  of 
a  periodic  function  can  be  determined  by  sampling  the  function  at  a  fixed  rate,  Tg .  Thus,  the 
Doppler  frequency  can  be  calculated  using  the  N  samples  and  a  Fast-Fourier  Transform  (FFT). 
The  accuracy  of  the  FFT  (and  thus  the  velocity  calculation)  increases  as  N  is  increased.  [10] 

The  sampling  time,  Tg ,  is  limited  by  the  bandwidth  over  which  the  BSA  is  trying  to  detect 
frequencies.  The  sampling  frequency,  fg  =  IfTg  ,  must  be  greater  than  the  bandwidth  and  is 
usually  set  proportional  to  the  bandwidth,  fg  =  n{BW) ,  where  n  is  some  number  greater  than  one 

and  BW  is  the  bandwidth.  It  should  be  noted  that  n  is  not  an  adjustable  variable;  it  is  set 
internally  within  the  BSA.  The  limitation  on  Tg  means  that,  as  the  record  length  is  increased,  the 
total  sampling  time  (also  known  as  the  record  interval)  is  increased.  In  high-speed  flows  where  a 
particle  only  takes  a  few  microseconds  to  traverse  the  measurement  volume,  the  accuracy  of  the 
FFT  becomes  limited  by  the  record  length,  as  the  total  sampling  time  must  be  less  than  the  time 
taken  to  traverse  the  control  volume. 

3.3  Particle  Validation 

The  velocity  calculation  assumes  the  particles  are  traveling  normal  to  the  plane  of  the 
fnnges  at  a  constant  velocity;  i.e.,  the  Doppler  burst  is  required  to  be  a  periodic  function  in  order 
for  the  FFT  to  accurately  calculate  the  velocity.  Therefore,  it  is  desirable  to  neglect  particles 
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Figure  3.4:  “Good”  vs.  “Bad”  Particles 


which  do  not  enter  the  control  volume  at  an  angle  normal  to  the  fringe  disks  or  which  may  be 
accelerating  through  the  control  volume.  It  is  also  desirable  to  filter  out  particles  which  may  only 
cross  the  control  volume  near  the  edges  (see  Figure  3.4),  as  well  as  bursts  received  from  two 
particles  crossing  the  control  volume  at  the  same  time.  In  order  to  filter  out  these  effects,  different 
types  of  particle  and  signal  validation  are  employed  by  the  BSAs. 


3.3.1  Fringe  Count.  In  order  to  filter  out  the  particles  that  traverse  near  the  edges  of  the 
control  volume  and  those  which  are  traveling  almost  parallel  to  the  plane  of  the  fiinge  disks,  a 
minimum  number  of  fiinge  crossings  is  required  for  the  particle  to  be  considered  valid.  This 
restricts  valid  particles  to  those  traveling  through  the  center  portion  of  the  control  volume  and 
nearly  normal  to  the  fringe  plane. 

3.3.2  Acceleration  Detection.  The  velocity  of  each  particle  is  calculated  assuming  the 
particle  is  not  accelerating  through  the  control  volume.  This  assumption  is  justified  by  filtering  out 
accelerating  particles;  this  filtering  is  accomplished  by  comparing  the  time  taken  to  cross  a  certain 
number  of  fringes.  For  instance,  the  time  taken  to  cross  three  fringes  is  compared  to  the  time  taken 
to  cross  five  fringes;  if  they  are  the  same,  the  particle  is  non-accelerating  and  is  accepted  as  valid. 
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3.3.3  Local  Maxima  Check.  When  the  BSA  is  performing  the  time  sampling,  it  has  no 
way  of  knowing  if  the  scattered  light  is  coming  off  one  or  more  particles.  For  instance,  if  two 
particles  are  traversing  the  control  volume  nearly  simultaneously,  they  will  produce  the  burst 
pattern  shown  in  Figure  3.5.  In  order  to  minimize  the  chances  of  the  BSA  tracking  two  particles  at 
once,  the  BSA  compares  the  intensity  pattern  of  light  gathered  at  each  sample.  More  specifically, 
it  compares  the  two  largest  local  maxima  to  the  global  maximum;  the  global  maximum  must  be  at 
least  four  times  the  local  maxima  in  order  for  the  burst  to  be  accepted  as  valid  (see  Figure  3.6).  A 
burst  pattern  like  the  one  shown  in  Figure  3.5  would  produce  a  signal  violating  the  local  maximum 
criteria,  causing  the  BSA  to  reject  the  burst. 
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Figure  3.6:  Local  Maxima  Validation  (from  Ref  [10]) 
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3.3.4  Envelope,  Pedestal,  and  Quality  Factor.  The  envelope  and  pedestal  are  properties 
of  the  Doppler  burst  received  by  the  BSA.  The  pedestal  is  generated  by  running  the  burst  signal 
through  a  baseline  clamp  and  a  lowpass  filter.  The  envelope  is  generated  by  running  the  burst 
signal  through  a  bandpass  filter,  a  rectifier,  and  finally  a  lowpass  filter  (see  Figure  3.7).  The 
quality  factor  is  the  ratio  between  the  envelope  and  the  pedestal;  it  sets  a  limit  on  the  amplitude  of 
the  received  signal,  below  which  the  burst  is  rejected.  Adjustment  of  the  envelope,  pedestal,  and 
quality  factor  allows  for  filtering  of  random  noise,  such  as  laser  light  which  is  scattered  by  defects 


Figure  3.7:  Envelope  and  Pedestal  Generation  (from  Ref  [10]) 
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in  the  wind  tunnel  windows.  Proper  adjustment  assures  only  bursts  with  high  signal-to-noise  ratios 
(SNRs)  are  used  for  velocity  calculations. 

3. 4  Directional  Ambiguity  and  Frequency  Shift 

When  a  particle  traveling  through  the  control  volume  produces  a  burst  signal,  the  BSA 
cannot  tell  whether  the  particle  is  moving  “up”  or  “down”  (referenced  to  the  orientation  in  Figure 
3.2)  through  the  control  volume;  the  Doppler  frequency  depends  only  on  the  spacing  between  the 
fringes  and  not  on  the  particle’s  direction  of  motion.  This  inability  to  distinguish  the  direction  of 
motion  is  referred  to  as  directional  ambiguity. 

The  problem  of  directional  ambiguity  is  solved  by  passing  one  of  the  two  laser  beams 
through  a  Bragg  cell,  which  introduces  a  frequency  shift  of  magnitude  f  into  the  beam.  This 
frequency  shift  causes  the  fringes  to  sweep  across  the  control  volume  at  a  constant  velocity, 
Vf  =  Xff  [20].  The  motion  of  the  fringes  causes  a  frequency  shift  in  the  burst  signal  that  is 

dependent  on  the  direction  of  the  particle’s  motion,  overcoming  the  problem  of  directional 
ambiguity. 

3.5  Particle  Correlation  and  Coincidence  Window 

When  using  a  two-component  (four-laser  beam)  LDV  system  to  simultaneously  measure 
particle  velocities  in  two  dimensions,  the  velocity  in  one  dimension  is  measured  by  one  frequency  of 
light,  while  a  different  frequency  of  laser  light  is  used  for  the  other  direction.  A  particle  will  then 
produce  two  Doppler  bursts,  one  corresponding  to  each  frequency  of  light.  As  a  result,  two  BSAs 
are  needed  to  process  the  burst  signals,  one  for  each  frequency.  The  signals  at  the  two  frequencies 
are  independent  of  one  another,  so  it  becomes  necessary  to  find  a  way  of  correlating  the  two.  That 
is,  a  burst  on  one  BSA  needs  to  be  paired  with  the  burst  on  the  other  BSA  which  was  produced  by 
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the  same  particle.  This  correlation  is  accomplished  using  a  coincidence  window.  The  coincidence 
window  is  a  time  interval  set  by  the  user  of  the  LDV  system  to  determine  if  two  bursts  are 
produced  by  the  same  particle. 

As  a  particle  enters  the  control  volume,  one  of  the  two  BSAs  will  detect  the  particle  first. 
Because  different  frequencies  of  laser  light  are  used  for  each  velocity  component,  the  intensity  of 
scattered  light  from  one  of  the  beams  will  generally  be  higher  than  that  from  the  other  beam.  The 
BSA  detecting  the  higher-intensity  beam  is  used  as  a  “master”  BSA:  when  it  detects  a  particle,  it 
signals  the  other  BSA  to  begin  looking  for  the  same  particle.  When  (and  if)  the  second  BSA  finds 
the  particle,  it  records  its  own  arrival  time.  The  coincidence  window  is  the  maximum  allowable 
time  between  the  arrival  time  on  one  BSA  and  the  arrival  time  on  the  second  BSA.  The 
coincidence  window  needs  to  be  set  large  enough  that  the  two  bursts  produced  by  the  same  particle 
are  correlated,  but  yet  small  enough  that  the  burst  recorded  on  one  BSA  is  not  correlated  with  more 
than  one  burst  on  the  other  BSA.  The  size  of  the  coincidence  window  is  also  limited  by  the  time 
taken  by  a  particle  to  cross  the  control  volume:  if  the  coincidence  window  is  too  large,  there  exists 
the  possibility  of  correlating  a  particle  inside  the  control  volume  with  one  outside  the  control 
volume.  This  artificial  correlation  of  particle  velocities  can  result  in  skewed  measurements  of  the 
product  u'  v" . 

3. 6  Angular  Bias  of  Turbulence  Measurements 

Measurements  of  d  v*  (as  well  as  the  turbulence  intensities)  made  with  the  LDV  system 
can  be  skewed  by  a  phenomenon  known  as  angular  bias.  Angular  bias  is  the  result  of  operating 
the  LDV  system  with  the  fringe  plane  aligned  at  a  substantial  angle  relative  to  the  mean  flow 
(normally  45°).  [15]  Traditionally,  this  angular  orientation  was  required  to  resolve  the  problem  of 
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equipment  limitations  when  measuring  high-speed  flows.  When  the  lasers  are  aligned  in  this 
matmer,  the  velocity  measured  by  the  BSAs  is  mcos(45),  which  is  considerably  lower  than 
measuring  u  directly,  resulting  in  easier  and  more  accurate  detection  by  the  BSAs. 

The  effect  of  angular  bias  comes  into  play  in  highly  turbulent  flows.  If  the  magnitude  of 
the  v-component  of  velocity  is  the  same  order  of  magnitude  as  u  and  the  control  volume  is  oriented 
at  an  angle  to  the  mean  flow,  there  is  the  possibility  that  the  particles  may  run  parallel  to  the 
fringes  (see  Figure  3.8).  In  this  case,  the  particles  do  not  cross  enough  fringes  to  be  validated  and 
the  particle’s  signal  is  discarded.  Discarding  these  particles  results  in  erroneous  measurements  of 
the  turbulence  intensities  as  well  as  the  correlation  «'  v*  .  In  general,  angular  biasing  will  cause  a 
decrease  in  the  measured  values  of  the  turbulence  intensities  and  the  magnitude  of  «'  v*  [15]. 

Direction  of 


Figure  3.8:  Turbulence  Effect  on  Particle  Direction 

The  effect  on  the  product  u'  v"  is  most  evident  in  early  LDV  experiments  where  the  u  and 
V  velocities  were  not  measured  simultaneously  {cf.  Ref  [21]),  but  were  measured  separately  with  a 
single-component  system,  first  measuring  u  and  then  v.  The  two  results  were  correlated  in  a 
maimer  similar  to  that  which  is  used  to  determine  li  v*  from  slanted-wire  hot-wire  measurements. 

Later  LDV  measurements  partially  corrected  this  problem  using  Bragg  cells  and  frequency 
shifting  {cf.  Ref  [30]).  Instead  of  the  particle  crossing  the  fringes,  the  fringes  swept  across  the 
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path  of  the  particle,  resulting  in  enough  crossings  for  the  particle  to  be  validated.  However,  in 
high-speed  flows,  if  the  fringe  velocity  is  not  fast  enough  the  particle  may  still  cross  the  control 
volume  without  crossing  enough  fringes  to  be  validated. 

Although  the  fringe  velocity  could  be  increased  by  increasing  the  frequency  shift  supplied 
by  the  Bragg  cell,  the  speed  with  which  the  fringes  can  sweep  across  the  control  volume  is  limited 
by  the  sensitivity  of  the  receiving  optics  and  the  processing  capabilities  of  the  BSAs.  According  to 
Elena  [15],  “to  retain  a  high  precision  of  measurements  with  the  presently  available  counters,  a 
frequency  of  45  MHz  must  not  be  exceeded.” 

Finally,  the  most  recent  experiments  have  used  two-component  LDV  systems,  where  the 
fringes  from  the  second  set  of  lasers  are  perpendicular  to  the  fringes  from  the  first  set  (cf.  Ref 
[16]).  Thus,  a  particle  rutming  parallel  to  one  set  of  fringes  and  not  validated  by  that  BSA  should 
be  detected  and  validated  by  the  other  BSA.  However,  if  the  lasers  are  still  placed  at  a  45°  angle  to 
the  mean  flow,  the  light  scattered  from  the  low-intensity  laser  beam  may  not  be  powerful  enough 
for  its  BSA  to  detect,  and  the  influence  on  ti  V  would  once  again  be  lost. 

The  effect  of  angular  bias  can  be  resolved  by  aligning  the  fnnge  disks  so  the  fringes  from 
the  highest-intensity  laser  beams  run  normal  to  the  mean  flow.  Thus,  even  in  highly-turbulent 
flows,  the  fastest  particles  will  cross  the  fringes  generated  by  the  high-intensity  laser  at  an  angle  no 
less  than  45°.  Any  particle  running  nearly  parallel  to  the  high-powered  beams  should  be  moving 
slow  enough  for  the  lower-intensity  beams  to  detect. 

3. 7  Intensity  Pattern  of  Scattered  Light 

In  order  to  minimize  the  number  of  bursts  rejected  because  the  intensity  of  scattered  light  is 
too  low,  it  is  necessary  to  place  the  receiving  optics  where  the  intensity  is  the  highest.  Figure  3.9 
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dp<A,  dp~X  dp>A. 

Figure  3.9:  Intensity  of  Scattered  Light  [Dantec,  Inc.] 


shows  the  theoretical  scattering  diagram  for  the  intensity  of  light  scattered  from  a  spherical 
particle.  Note  that  the  intensity  distribution  depends  on  the  wavelength  of  light  as  well  as  the  size 
of  the  particle  in  Figure  3.9  represents  the  particle  diameter). 

It  is  clear  from  Figure  3.9  that  the  intensity  distribution  is  generally  the  highest  on  the 
“front”  side  of  the  particle.  Durst  etal.  [14]  presented  the  scattered  light  intensity  on  a  radial  log 
scale;  this  graphic  is  reproduced  in  Figure  3. 10.  Note  that  the  intensity  distribution  is  plotted  on  a 
logarithmic  scale,  and  as  such  the  intensity  on  the  front  side  of  the  particle  can  be  orders  of 
magnitude  higher  than  that  on  the  back  side.  Therefore,  in  high-speed  flows  where  the  scattered 
light  has  relatively  low  intensities  in  all  directions,  the  receiving  optics  are  placed  on  the  far  side  of 
the  measuring  volume  from  the  transmitting  optics,  an  alignment  referred  to  as  forward  scattering 
mode. 


240®  280®  S20® 


Figure  3.10:  Magnitude  of  Intensity  of  Scattered  Light  [14] 
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While  use  of  forward  scattering  mode  enables  the  collection  of  more  scattered  light,  it  also 
increases  the  level  of  difficulty  associated  with  data  collection.  The  receiving  optics  are  contained 
in  a  separate  apparatus  from  the  transmitting  optics  and  must  be  precisely  aligned  with  the  control 
volume  (when  using  backscatter  mode  with  the  current  LDV  system,  the  transmitting  and  receiving 
optics  are  contained  in  the  same  probe  and  pre-aligned  to  be  focused  on  the  same  location:  the 
process  becomes  a  matter  of  “point-and-shoot”).  When  the  transmitting  optics  are  re-positioned, 
the  receiving  optics  need  to  be  realigned  as  well. 

3.8  Seeding  the  Flow:  Intrusive  Aspects  of  LDV  and  the  Problem  of  Velocity  Slip 

To  this  point,  it  has  been  assumed  that  the  flow  contains  particles  which  will  scatter  the 
light  from  the  laser  beam;  however,  it  generally  necessary  to  introduce  these  particles  into  the 
flow,  a  process  known  as  seeding.  Even  when  there  are  naturally-occurring  particles  in  the  flow,  it 
is  better  to  filter  the  natural  particles  and  introduce  particles  of  known  size  and  material  properties 
[15].  The  necessity  of  introducing  these  particles  means  that  LDV  is  not  really  non-intrusive,  as  is 
often  suggested.  In  order  to  reduce  the  effects  of  seeding,  several  precautions  must  be  taken.  The 
most  important  aspects  are  efficient  injector  design,  proper  particle  sizing  and  particle  material 
considerations. 

3.8.1  Injector  Design.  Ideally,  an  injector  would  not  have  to  be  placed  inside  the  flow;  if 
possible,  an  opening  should  be  placed  in  the  wall  of  the  wind  tunnel,  through  which  seed  particles 
are  added  to  the  flow.  Carefixl  attention  needs  to  be  placed  on  the  sizing  of  the  hole,  as  the  hole 
diameter  can  affect  the  flow  characteristics  [15].  If  the  design  of  the  wind  tunnel  does  not  permit 
injection  from  an  opening  in  the  wall,  some  type  of  injector  needs  to  be  placed  in  the  flowfield. 
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Any  injector  which  is  placed  in  the  flowfield  will  have  a  wake  associated  with  its  presence. 
Therefore,  it  is  desirable  to  place  the  injector  far  upstream  from  the  test  section  and,  if  possible,  in 
a  region  where  the  flow  velocities  are  low  (to  minimize  the  wake  produced).  If  the  injector  were 
placed  in  a  region  of  supersonic  flow,  its  presence  would  induce  a  shock  wave.  Thus,  in  a 
supersonic  wind  tunnel,  it  is  necessary  to  place  the  injector  before  the  throat,  when  the  flow  is  still 
subsonic.  In  all  cases  the  injector  should  be  streamlined  to  reduce  the  wake  produced. 

No  matter  what  kind  of  injector  is  used,  it  is  necessary  to  inject  the  particles  at  a  speed 
equal  to  the  local  mean  flow  velocity  in  order  to  minimize  pressure-difference  effects  [15].  It  is 
also  necessary  to  use  the  lowest  seeding  density  possible  in  order  to  minimize  the  effects  of  the 
particles  on  the  flow  being  studied. 

3.8.2  Particle  Sizing.  Careful  attention  must  be  paid  to  the  size  of  the  particles  injected 
into  the  flow;  particles  must  be  small  enough  to  ensure  that  they  follow  the  smallest  motions  of  the 
fluid.  This  becomes  especially  important  in  the  study  of  high-speed  turbulence  and  flows  near 
shock  waves.  Particles  which  do  not  follow  the  flow  are  said  to  experience  velocity  slip. 

3.8.3  Particle  Material  Considerations.  Particles  should  also  be  chosen  according  to 
material  properties.  When  studying  flows  through  shock  waves  it  is  desirable  to  choose  particles 
with  a  fast  response  time;  other  times  it  may  be  advantageous  to  choose  particles  which  will  not 
cling  to  and  cloud  up  test  section  windows.  In  other  experiments,  such  as  low-density  flows,  the 
density  of  seed  particles  may  be  high  enough  to  actually  change  the  material  properties  of  the  flow. 
Therefore,  it  is  necessary  to  choose  particles  which  will  have  the  least  effect  on  the  flow  properties 
and  still  follow  the  flow  as  closely  as  possible. 
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4.  Equipment  Description 


4.1  Wind  Tunnel  Facilities 

Data  were  collected  in  the  AFIT  Mach  3.0  wind  tunnel.  The  tunnel  is  a  combination  draw- 
down/blow-down  facility;  see  Figure  4. 1  for  a  diagram  of  the  facilities.  A  Leslie  brand  dome 
regulator  was  used  to  provide  manual  control  of  the  stagnation  pressure;  the  upstream  Pitot 
pressure  was  read  using  a  Endevco  brand,  0-0.69  MPa  pressure  transducer  and  was  set  at  a 
constant  value  of  2.13 x  10^  +690 Pa.  The  total  temperature  was  monitored  with  an  Omega 
Engineering  type  K  thermocouple  placed  just  upstream  of  the  wind  tunnel  nozzle.  The  total 
temperature  was  constant  during  a  single  run  to  approximately  ±0.1  K . 

4.1.1  High-Pressure  Air  Supply.  The  high  pressure  was  supplied  by  two  Atlas  Compco 
compressors  which  could  provide  a  combined  continuous  mass  flow  rate  of  approximately  0.45 
kg/s  at  0.69  MPa.  The  high  pressure  system  supplied  air  to  other  facilities  in  AFIT  besides  the 
Mach  3.0  system,  resulting  in  slight  variations  in  the  available  pressure.  The  pressure  was 
monitored  during  each  traverse  across  the  boundary  layer  to  ensure  it  remained  constant. 

High  Pressure 


Figure  4. 1 :  Schematic  of  Wind  Tunnel  Facilities 
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The  high-pressure  air  was  dried  using  two  Pioneer  Air  Systems,  Inc.  model  R500A 
refrigeration-type  dryers.  It  was  then  passed  through  a  centrifiigal  moisture  and  particle  separator, 
and  filtered  through  a  layer  of  filament-reinforced  industrial  paper  filters.  The  filters  were 
expected  to  remove  any  naturally-occurring  particles  in  the  flow. 

4.1.2  Vacuum  System.  The  vacuum  was  created  using  three  Stokes  Penwalt  model  212- 
1 1  MicroVac  pumps.  Each  pump  was  powered  by  a  7.6  metric  hp  Reliance  motor.  The  vacuum 
tank  capacity  of  the  system  was  approximately  16  cubic  meters.  The  vacuum  pumps  were  able  to 
reduce  the  pressure  in  the  tanks  to  10  mm  of  Hg  in  approximately  six  minutes. 

4.1.3  Wind  Tunnel  Operation.  Flow  through  the  tunnel  was  initiated  by  opening  the 
Kinney  Corp.  GP-8  gate  valve  in  order  to  expose  the  tuimel  to  the  vacuum.  The  high  pressure  was 
then  released  by  opening  the  El-O-Matic  ED  10  ball  valve,  starting  the  airflow  through  the  tunnel. 
Once  airflow  was  initiated,  the  tunnel  took  approximately  eight  seconds  to  reach  steady  state, 
followed  by  a  25-second  steady-state  period  in  which  it  was  possible  to  collect  data. 

4.1.4  Test  Section.  The  test  section  had  a  cross-sectional  area  6.35  cm  x  6.35  cm.  The 
floor  and  wall  sections  were  constructed  of  1.905  cm  aluminum  alloy  and  were  modular  to  allow 
for  interchange  of  various  models.  The  test  section  was  actually  comprised  of  two  sub-sections 
(see  Figure  4. 1);  each  sub-section  was  33.02  cm  in  length.  All  seams  in  the  test  section  walls  were 
fitted  with  rubber  gaskets  which  sealed  the  wall  sections  when  the  vacuum  system  was  engaged. 
The  sections  were  also  fitted  with  adjustable  flanges  which  allowed  each  section  to  be  shifted 
relative  to  the  adjacent  section  in  order  to  minimize  shocks  and  expansions  that  might  have  formed 
off  the  seams.  The  FPG  model  was  integrated  into  the  ceiling  of  the  second  test  section  as  shown 
in  Figure  4.1;  the  specifics  of  the  design  are  discussed  below. 
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4. 1.4.1  Model  Design.  The  FPG  model  consists  of  a  curved  ceiling  for  the  test 


seetion;  the  equation  for  the  curve  is  given  in  test  section  coordinates  by  [26] 


yts  = 


0 


0<x^<4.32 


fl-o  +  +  a^xl  +  4.32  <  <  12.7 

-0.65 


(4.1) 


12.7<x^<  16.51 


where  the  units  are  in  cm  and  the  coefficients  are  given  in  Table  4.1.  It  should  be  noted  that  the 
coordinates  for  this  curve  were  also  reported  in  Ref  [27],  but  the  coordinates  in  Ref  [27]  were 
incorrect.  The  accuracy  of  the  wall  position  equation  is  ±0.0008  cm .  Figure  4.2  shows  a  scale 
drawing  of  the  test  section  ceiling. 


Table  4.1:  Test  Section  Curve  Coefficients 


«0 

«2 

cfj  X  1 000 

-0.2078 

0.0897 

-0.0095 

-0.0360 

Figure  4.2;  Scale  Schematic  of  Test  Section  Ceiling 

With  the  equation  for  the  curve  known  in  test  section  coordinates,  the  relationship  between 
the  body-intrinsic  coordinate  x  and  the  test  section  coordinate  x^^  was  given  by 

^  I - 

x  =  64.32  ±  jyjl  +  {dy^/dx,^f  dx,^  (4.2) 

4.32 

A  flat  plate  model  was  also  used  in  order  to  collect  data  in  the  same  location  as  the  FPG;  the  flat 
plate  and  favorable  pressure  gradient  data  were  collected  at  the  location  x,;^  =  1 1.5  cm  {x^,  =  71.5 
cm).  Flat  plate  data  were  also  collected  at  x^^  =  44  cm,  in  the  center  of  the  upstream  sub-section, 
in  order  to  compare  with  the  data  collected  by  Miller.  [26] 
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4. 1.4. 2  Window  Placement.  Use  of  the  LDV  system  required  the  tunnel  test 
section  walls  to  be  fitted  with  optical-grade  glass  windows.  The  windows  chosen  were  made  of 
1.27  cm  thick  optical  grade  glass.  Each  window  was  circular,  with  a  diameter  of  7.62  cm.  The 
windows  were  placed  in  the  test  section  wall  as  shown  in  Figure  4.3.  Note  that  the  test  seetion 
floor/ceiling  had  a  length  of  33.02  cm,  while  the  wall  was  only  29.21  cm  long.  The  difference  was 
due  to  the  fact  that  the  adjustable  flanges  at  the  tunnel  seams  were  designed  slightly  differently  for 
the  floor  and  ceiling  sections  than  for  the  wall  sections. 


-  29.21  - 

Figure  4.3:  Optical  Window  Placement  Schematic  (Dimensions  in  cm) 


Note  that  the  left  window  in  Figure  4.3  was  slightly  off-center  in  the  vertical  direction. 
The  displacement  was  incorporated  to  allow  measurements  over  a  larger  portion  of  the  test  section; 
more  importantly,  it  allowed  the  window  to  expose  the  section  of  the  FPG  model  where  it  curved 
back  “into”  the  test  section  wall.  Figure  4.4  shows  the  total  possible  area  covered  by  the  test 
section  windows  and  the  location  of  the  test  section. 
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Probe  Position 
(x,y,z)  =  (0AO) 


4.2  LDV System 

The  laser  Doppler  velocimetry  system  was  manufactured  by  Dantec,  Inc.  Many  parts  of 
the  system  were  interchangeable,  so  the  details  and  specifications  of  the  equipment  used  in  this 
study  are  presented  below. 

4.2.1  Laser  Optics.  The  laser  optic  system  consisted  of  the  laser  itself,  a  transmitter 
which  sent  the  laser  beam  to  the  transmitting  optics,  the  transmitting  optics,  the  receiving  optics, 
and  the  photomultiplier  (PM)  tubes.  Each  component  is  described  in  detail  below. 

The  laser  was  an  Ion  Laser  Technology,  Model  5500A-00  Argon-ion  laser  with  a 
maximum  (rated)  output  capability  of  300  mW.  Due  to  the  deterioration  of  efficiency  which 
occurs  with  age,  the  laser  was  found  to  have  a  maximum  output  power  of  approximately  275  mW. 
The  laser  was  supplied  with  power  using  an  Ion  Laser  Technology  Model  5405A  power  supply. 
The  power  supply  operated  at  210V  and  20 A,  requiring  a  220  VAC  power  source.  The  laser  was 
designed  to  operate  in  the  TEMgg  mode,  where  the  laser  beam  is  centered  around  the  optical  axis 

with  the  maximum  power  in  the  center.  The  l/e^  beam  diameter  of  the  laser  was  0.82  mm. 
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The  transmitter  was  a  Fiber  Flow  60X41  transmitter.  A  Bragg  cell  in  the  transmitter  split 
the  laser  into  three  wavelengths  of  light:  514.5  nm,  488.0  nm„  and  476.5  nm.  These  frequencies 
corresponded  to  light  beams  which  appeared  green,  blue,  and  violet,  respectively.  The  Bragg  cell 
also  split  each  wavelength  into  two  separate  beams  and  supplied  one  of  the  two  with  a  40  MHz 
frequency  shift.  As  described  in  Chapter  3,  the  frequency  shift  caused  a  motion  of  the  fringes 
across  the  control  volume.  Note  that  the  frequency  shift  supplied  by  the  Bragg  cell  is  less  than  the 
maximum  limit  (45  MHz)  given  in  Ref  [15]. 

The  transmitter  was  equipped  with  six  DANTEC  60X24  fiber  manipulators  which  were 
used  to  attach  the  fiber-optic  cables  to  the  transmitter.  The  manipulators  allowed  the  ends  of  the 
fiber  optic  cables  to  be  precisely  adjusted  in  order  to  ensure  the  laser  beam  was  centered  on  the 
cable,  providing  maximum  transmitting  efficiency.  The  laser  beams  were  then  passed  through  the 
fiber  optic  cables  to  the  transmitting  and  receiving  optics.  The  laser  and  transmitter  are  shown  in 
Figure  4.5. 

The  transmitting  optics  were  contained  in  a  DANTEC  60X61  optic  probe,  while  the 
receiving  optics  were  in  a  model  60X60  optic  probe  (shown  in  Figure  4.6).  Each  probe  contained 
an  optic  transducer,  a  distributor  unit,  and  transmitting  and  receiving  fibers.  The  transmitting 
probe  was  a  two-component  probe  which  projected  the  514.5  and  488  nm  laser  beams,  while  the 
receiving  probe  was  single-component  and  projected  the  476.5  nm  beam.  Each  probe  had  two 
transmitting  fibers  for  each  component;  therefore,  the  transmitting  optics  were  projecting  four 
laser  beams  while  the  receiving  optics  were  projecting  two.  Each  probe  was  attached  to  the 
traverse  mounting  bench  using  DANTEC  probe  supports.  The  transmitting  optics  were  placed  on 
a  60X361 1  support  which  allowed  for  angular  adjustments  in  the  Xj,  -y-j.  plane,  while  the  receiving 
optics  were  mounted  on  a  60X3631  probe  support  which  used  micrometers  to  allow  precise 
adjustment  in  all  three  dimensions. 
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On  the  front  of  the  transmitting  probe  was  a  600  mm  focal-length  lens;  the  laser  beams 
were  separated  by  38  mm  on  the  front  face  of  the  lens  (see  Figure  4.7).  The  diameter  of  the  lasers 
beams,  the  focal  length  of  the  lens  on  the  transmitting  optics  and  the  beam  separation  on  the  face  of 
the  transmitting  probe  interacted  to  produce  a  measurement  volume  9.0  mm  in  length  (along  the  z- 
axis)  and  0.276  mm  in  diameter.  The  receiving  optics  were  fitted  with  a  55X12  beam  expander  as 
well  as  the  600  mm  focal  length  lens.  The  beam  expander  was  used  to  increase  the  amount  of 
reflected  light  gathered  by  the  optics. 


Figure  4.7;  Beam  Separation  on  Face  of  Optic  Probe 


The  light  collected  by  the  receiving  optics  was  sent  through  fiber  optic  cables  to  a 
DANTEC  55X35  color  separator.  The  color  separator  split  the  received  light  into  514.5  and  488.0 
nm  wavelengths  and  transferred  each  wavelength  to  a  DANTEC  55X35  photomultiplier  (PM) 
tube.  The  PM  tube  converted  the  received  light  into  an  electrical  signal  associated  with  the 
frequency  of  the  light.  This  electrical  signal  was  then  sent  to  the  BSA  for  processing. 


4.2.2  Data  Acquisition  Equipment.  The  Doppler  bursts  detected  by  the  receiving  optics 
were  processed  by  one  Dantec  57N20  BSA  Enhanced  (w-component  of  velocity)  and  one  Dantec 
57N35  BSA  Enhanced  Model  S  (v-component  of  velocity).  The  57N20  BSA  was  used  as  the 
master  BSA  for  coincidence  filtering  calculations;  the  interior  clock  had  a  frequency  of  1  MHz. 

The  software  used  in  the  data  collection  was  the  Dantec  Burstware®  2.0  package;  the 
software  was  run  on  a  Gateway  2000  486DX/33  personal  computer. 
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4.2.3  Traverse.  A  Dantec  3D  Traverse  System  was  used  to  control  the  position  of  the 
measurement  volume.  The  traverse  was  powered  by  three  Cleveland  Machine  Controls  stepper 
motors  which  were  regulated  by  Dynapar  encoders,  model  M20100003331.  The  encoders 
converted  1000  pulses  into  one  revolution  of  the  load  screw,  resulting  in  a  2  mm  displacement  of 
the  traverse.  The  entire  system  was  rated  at  an  accuracy  of  ±80  pm  over  a  600  mm  range  (0.13 
pm/mm)  and  could  support  a  maximum  load  of  445  N.  The  DANTEC  traverse  mounting  bench 
was  replaced  by  a  1.83  m  long  DANTEC  laser  mounting  bench  for  reasons  which  will  be  described 
in  Chapter  5. 

4.2.4  Seeding  Apparatus .  The  seeding  apparatus  consisted  of  a  TSI,  Inc.  six-jet  atomizer 
attached  to  an  injector  which  was  placed  inside  the  settling  chamber  of  the  wind  tunnel.  The 
connection  from  the  atomizer  to  the  injector  was  accomplished  using  Tygon-brand  plastic  tubing 
(type  R-3603,  1.27  cm  inner  diameter). 

4.2.4. 1  Atomizer.  The  atomizer  was  manufactured  by  TSI,  Inc.  and  is  shown  in 
Figure  4.8.  It  was  a  variable-flow  atomizer,  allowing  for  particle-density  control  through  two 
options.  The  first  was  to  use  a  pressure  control  valve  to  regulate  the  pressure  difference  through 
the  atomizer.  The  second  method  for  particle  density  control  was  manual  control  over  the  number 
of  jets  in  use. 

4. 2. 4. 2  Seeding  Material.  Due  to  cost  and  availability  concerns,  the  seeding 
material  chosen  for  use  in  this  experiment  was  Bertoli  brand  extra  light  olive  oil.  The  atomizer 
manual  [34]  gave  the  mean  particle  size  for  olive  oil  particles  as  d^  «  0.6  pm. 
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Figure  4.8:  Six-Jet  Atomizer 
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4.2.4. 3  Injector  Design.  The  injector  was  designed  from  copper  tubing  with  an 
inner  diameter  of  1.27  cm.  In  order  to  minimize  the  effect  of  the  injector  on  the  flow,  the  copper 
tubing  was  streamlined  by  flattening  it  in  a  vise.  Openings  were  cut  into  the  injector  such  that 
particles  would  only  be  injected  into  the  flow  over  the  upper  half  of  the  wind  tunnel.  The  injector 
was  placed  in  the  settling  chamber,  about  3  cm  upstream  of  the  flow  straightener.  A  schematic  of 
the  injector  and  its  attachment  to  the  atomizer  is  shown  in  Figure  4.9. 


Side  View  Top  View 


Seed  From  Atomizer 

Figure  4.9:  Injector  Design  and  Location  in  Wind  Tunnel 

4. 3  Flow  Visualization  Equipment 

The  shadowgraphs  and  schlierens  were  taken  using  a  Xenon  Corp.  43 7B  Nanopulser 
which  fed  power  to  a  De  Leone  Corp.  N-787B  broad-spectrum  spark  generator.  The  light  from  the 
spark  was  focused  using  one  or  two  14.4  cm  diameter  focal  mirrors,  each  with  a  focal  distance  of 
152.4  cm.  In  the  case  of  the  schlieren  photographs  (see  Chapter  5  for  details  of  the  arrangement), 
the  light  beam  was  also  passed  through  a  1 1.2  cm  diameter  focusing  lens.  The  photographs  were 
taken  with  a  Polaroid  545/  camera  using  Polaroid  type  57  high-speed,  ISO  3000/36°  film. 
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5.  Experimental  Procedures 


5. 1  Wind  Tunnel  Operation 

Airflow  was  initiated  in  the  wind  tunnel  as  described  in  Section  4.3.1.  During  the  eight- 
second  span  when  the  wind  tunnel  was  approaching  steady-state  operation,  two  jets  on  the  atomizer 
were  opened  to  begin  injection  of  seed  into  the  flow.  Once  the  steady-state  point  was  reached,  data 
were  collected  for  12  seconds,  collecting  as  many  bursts  as  possible.  Operation  of  the  wind  tuimel 
was  stopped  by  closing  the  high-pressure  supply  and  then  shutting  the  gate  valve  to  seal  off  the 
vacuum  tanks. 

5. 2  LDV  System  Preparation 

Due  to  the  relatively  low  intensity  of  the  laser  beams,  the  high  flow  velocities  and  the 
associated  low  SNR,  the  LDV  system  was  operated  in  forward-scatter  mode.  The  two-component 
probe  was  used  as  the  transmitter,  white  the  single-component  probe  functioned  as  the  receiver.  In 
order  to  simplify  data  collection,  it  was  desired  to  have  both  the  transmitting  and  receiving  optics 
attached  to  the  traverse  system;  otherwise,  whenever  the  transmitting  optics  were  moved  to  a  new 
location,  the  receiving  optics  would  need  to  be  realigned.  In  order  to  overcome  this  problem,  both 
the  transmitting  and  receiving  optics  needed  to  be  placed  on  the  traverse  mounting  bench,  but  the 
focal  distance  of  the  optics  was  too  long  to  allow  this  to  happen.  The  problem  was  solved  by 
interchanging  the  mounting  bench  on  the  Dantec  traverse  system  with  the  laser  mounting  bench; 
the  laser  mounting  bench  was  long  enough  to  allow  both  the  transmitting  and  receiving  probes  to  be 
placed  on  the  traverse  and  moved  simultaneously.  Interchanging  the  two  benches  forced 
realignment  of  the  laser  transmitter  with  the  laser,  ensuring  the  laser  system  was  operating  at  peak 
efficiency. 
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After  the  optic  probes  were  attached  to  the  traverse  system,  the  traverse  needed  to  be 
aligned  with  the  wind  tunnel  and  the  optics  with  the  test  section  wall;  in  addition,  the  receiving 
optics  had  to  be  aligned  with  the  control  volume  and  the  BSA  settings  had  to  be  optimized. 

5. 2. 1  Traverse  Alignment,  Laser  Orientation  and  Origin  Placement 
In  order  to  prepare  for  data  collection,  the  traverse  had  to  be  aligned  with  the  test  section. 
Once  the  traverse  was  aligned,  the  optics  needed  to  be  aligned  with  the  test  section  wall  in  order  to 
collect  data  in  a  body-intrinsic  coordinate  system. 

5.2.1. 1  Traverse  Alignment.  To  simplify  conversion  between  the  traverse 
coordinate  system  and  the  test  section  coordinates,  the  traverse  bench  was  aligned  perpendicular  to 
the  wind  tunnel  (see  Figure  5.2).  The  alignment  was  accomplished  by  pointing  the  lasers  parallel 
to  the  traverse  bench  (the  DANTEC  probe  supports  were  marked  at  1°  intervals,  allowing  the 
lasers  to  be  accurately  aligned  with  the  traverse  bench)  and  rotating  the  entire  traverse  assembly 
until  the  lasers  were  reflected  off  the  windows  back  onto  themselves.  Once  the  traverse  was 
aligned  normal  to  the  tuimel,  the  lasers  were  rotated  3.5°  off  the  axis  of  the  traverse  bench  (i.e., 
3.5°  off  the  z-axis,  see  9  in  Figure  5.1).  This  was  necessary  to  ensure  that  the  extremely  sensitive 
PM  tubes  in  the  receiving  optics  were  not  exposed  to  direct  laser  light  from  the  transmitting  optics. 
While  this  alignment  caused  a  slight  bias  to  the  collected  data,  it  was  small  enough  to  be  neglected 
(see  Appendix  B). 
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Figure  5.1:  Traverse  and  Laser  Alignment  Relative  to  the  Wind  Tunnel 
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An  attempt  was  made  to  align  the  traverse  in  a  way  which  would  keep  the  transmitted 
lasers  normal  to  the  turmel  wall,  eliminating  the  bias  described  above.  Due  to  the  fact  that  the 
transmitted  lasers  could  not  be  pointed  directly  at  the  receiving  optics,  this  would  force  the  traverse 
bench  to  be  placed  so  that  it  was  no  longer  perpendicular  to  the  test  section.  Orienting  the  traverse 
in  this  maimer  would  complicate  calculation  of  traverse  coordinates  in  the  space,  but 

the  measurement  biasing  would  be  eliminated.  However,  use  of  this  orientation  was  not  possible, 
because  it  would  have  forced  the  vertical  support  of  the  traverse  to  be  placed  directly  below  the 
wind  tunnel  (which  was  not  physically  possible).  The  traverse  support  could  have  been  placed  to 
the  side  of  the  wind  tunnel  by  adjusting  the  traverse  mounting  bench  off-center  from  the  vertical 
support,  but  this  option  was  also  not  physically  feasible.  Recall  that  the  mounting  bench  on  the 
traverse  was  originally  the  laser  mounting  bench  and  was  larger  (and  heavier)  than  the  normal 
traverse  mounting  bench;  placing  the  larger  bench  off-center  would  cause  a  torque  on  the  traverse 
system  which  would  result  in  failure.  The  alignment  which  was  chosen  for  use  in  data  collection 
(where  the  traverse  bench  was  normal  to  the  tunnel)  also  forced  the  bench  to  be  placed  slightly  off- 
axis,  but  not  far  enough  to  cause  a  large  torque  on  the  system. 

Finally,  note  that  the  angle  q)  in  Figure  5. 1  is  given  by  (p  =  29  =  T .  For  particles  with  a 
size  on  the  order  of  the  wavelength  of  laser  light  (600  ran  in  this  study),  the  intensity  pattern  of 
scattered  light  is  given  by  the  center  pattern  in  Figure  3.9.  The  receiving  optics  were  7°  off-axis; 
thus,  the  relative  intensity  seen  by  the  receiving  optics  was  almost  maximum,  as  shown  in  Figure 
5.3. 


CP=^ 


Figure  5.3;  Scattered  Light  Intensity  Received  by  Optics 
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5.2. 1.2  Origin  Placement  and  Streamwise  Angular  Alignment.  Once  the 
traverse  was  aligned  normal  to  the  wind  tunnel,  the  origin  of  the  data-collection  coordinate  system 
needed  to  be  set.  The  -location  was  placed  in  the  correct  location  (all  tests  were  conducted  at 
=71.5  cm,  except  for  one  flat  plate  test  at  =  44  cm)  downstream  of  the  nozzle  by  aligning 
the  lasers  with  a  mark  on  the  outside  of  the  test  section  ceiling.  When  the  lasers  were  aligned  in 
this  direetion,  the  traverse  was  moved  in  the  -Zj  direction  (recall  that  Zj.  is  positive  into  the  wall 
of  the  tunnel)  until  the  lasers  were  near  the  wall  of  the  test  section. 

The  streamwise  laser-angle  alignment  and  placement  of  the  origin  in  the  Zj  direction  was 
a  simultaneous  process.  Using  a  pair  of  argon-ion  laser  protective  goggles,  the  reflection  of  the 
laser  beams  on  the  tuimel  wall  appeared  as  a  set  of  small  dots  which  were  located  at  the  center  of 
the  laser  beam,  where  the  intensity  was  maximum.  The  optics  were  rotated  and  the  traverse  was 
moved  in  the  Zj.  direction  until  the  two  green  laser  beams  (the  514.5  nm,  highest-intensity  beams) 

were  bisected  by  the  tunnel  wall;  this  location  was  set  as  {x,y,z)  =  (0,0,0) .  Note  that,  with  the 

lasers  aligned  in  this  manner,  the  w-component  of  velocity  ran  parallel  to  the  wall  while  the  v- 
component  was  normal  to  the  wall,  with  positive  v  pointing  away  from  the  wall.  Aligning  the 
lasers  in  this  orientation  was  expected  to  reduce,  if  not  eliminate,  the  effects  of  angular  bias 
reported  by  other  researchers.  [13] 

5.2.2  Alignment  of  Receiving  Optics.  After  the  transmitting  optics  were  aligned  with  the 
test  section  wall,  the  receiving  optics  needed  to  be  aligned  with  the  control  volume.  Using  the 
lasers  projected  from  the  receiving  optics,  the  control  volumes  of  the  transmitting  and  receiving 
optics  could  be  quickly  aligned  with  a  fair  degree  of  accuracy.  The  alignment  was  then  fine-tuned 
by  lowering  the  lasers  into  the  freestream,  injecting  a  small  amount  of  seed  into  the  wind  tunnel 
(while  the  wind  tunnel  was  not  in  operation)  and  using  the  “on-line  display”  feature  of  the  software. 
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The  on-line  display  allowed  for  constant  update  of  the  validity  and  data  acquisition  rates;  as  the 
receiving  optics  were  adjusted,  the  validity  and  data  rates  would  change.  The  micrometer  adjusters 
on  the  receiving  probe  allowed  the  receiving  optics  to  be  adjusted  very  accurately  in  order  to 
maximize  the  data  acquisition  and  validity  rates.  Using  this  technique,  the  maximum  data 
acquisition  rate  seen  during  the  alignment  process  was  on  the  order  of  4  kHz;  the  maximum 
validity  rate  during  the  alignment  process  was  found  to  be  approximately  92-94%  simultaneously 
on  each  BSA  (the  BSA  attached  to  the  high-intensity  lasers  was  able  to  achieve  a  higher  validity 
rate  than  this,  but  only  at  the  expense  of  the  second  component). 

5.2.3  Software  Settings.  The  final  step  before  actual  data  collection  was  to  adjust  the 
settings  on  the  BSAs  to  maximize  the  number  of  bursts  collected  while  ensuring  only  the  valid 
bursts  were  processed.  The  proper  setting  for  most  of  the  variables  depended  on  the  characteristics 
of  the  flow  field:  high  vs.  low  speed,  the  amount  of  turbulence  in  the  flow,  if  there  were  shocks 
present,  etc.  In  order  to  determine  the  effects  each  setting  had  on  the  collected  data  and  to  decide 
on  proper  settings  for  the  wind-tunnel  experiment,  preliminary  tests  were  conducted  using  a  high¬ 
speed  nozzle  attached  to  the  atomizer;  the  results  from  these  tests  helped  establish  the  settings  used 
in  the  wind-tunnel  tests. 

During  data  collection  in  the  wind  tunnel,  most  of  the  variables  were  kept  constant 
throughout  the  testing,  but  a  few  were  altered  from  one  data  point  to  the  next.  The  fixed  and 
adjustable  settings  are  presented  and  discussed  below.  It  should  be  noted  that  BSA  1  was 
collecting  the  data  from  the  5 14.5  nm  (green)  laser  beams,  while  BSA  2  collected  the  data  from  the 
488  nm  (blue)  laser  beams. 

5.2. 3.1  Fixed  Settings.  The  fixed  settings  were  quantities  such  as  the  envelope, 
pedestal,  and  quality  factor  which  were  kept  constant  throughout  all  tests.  The  constant  settings 
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used  in  the  data  collection  are  presented  in  Tables  5.1  to  5.3.  It  should  also  be  noted  that  the  data 
collection  mode  under  the  “Valid”  menu  was  set  to  “Valid  Data.” 


Table  5.1:  “Quick”  Menu  Constant  Software  Settings 


Category 

Bandwidth  (m/s) 

325.12 

246.7 

Shifter  Mode 

Norm 

Norm 

Signal  Gain  (dB) 

42 

42 

High  Voltage  (V) 

1304 

1504 

Max.  Anode  Current  (mA) 

1.6 

1.6 

Pedestal  Attenuation  (dB) 

6 

6 

Duty  Cycle  (%) 

100 

100 

Dead  Time 

0 

0 

Table  5.2:  “Soft”  Menu  Constant  Software  Settings 


Category 

Timer  Clock 

Master 

Slave 

Coincidence  Mode 

Master 

Slave 

Arrival  Time  Base 

Internal 

Internal 

Burst  Detection 

Both 

Both 

Oversize  Rejection 

1 

1 

Buffer  Mode 

FIFO 

FIFO 

Max.  Anode  Current  (mA) 

1.6 

1.6 

Quality  Factor  (%) 

50 

50 

Collection  Mode 

Burst 

Burst 

Table  5.3:  “BSA  Program”  Menu  Constant  Software  Settings 


Category 

Timeout 

12 

12 

Number  of  Bursts 

150000 

150000 

Velocity 

Yes 

Yes 

Transit  Time 

Yes 

Yes 

Arrival  Time 

Yes 

Yes 
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5. 2. 3. 2  Adjustable  Settings.  Settings  which  were  adjusted  during  data  collection 
were  the  center  frequency  on  BSA  1,  the  record  length,  and  the  size  of  the  FIFO  buffer.  The  center 
frequency  on  BSA  1  needed  to  be  adjusted  as  the  control  volume  was  moved  across  the  boundary 
layer  and  the  measured  velocities  increased.  The  center  frequency  was  set  as  close  as  possible  to 
the  average  velocity  at  the  given  point. 

As  mentioned  in  Chapter  4,  the  control  volume  was  0.276  mm  long  in  the  streamwise 
direction;  a  particle  traveling  normal  to  the  control  volume  (as  was  required  by  the  validation 
processes  described  in  Chapter  3)  at  650  m/s  (the  highest  velocity  measured  in  the  current 
experiment  when  considering  fluctuations  from  the  mean)  would  take  0.42  ps  to  cross  from  one 
side  to  the  other.  Thus,  the  record  interval  had  to  be  less  than  0.42  ps.  Due  to  the  finite  sampling 
rates  used  by  the  BSAs,  the  limit  on  the  record  interval  placed  a  limit  on  the  record  length, 
restricting  it  to  8  or  16  samples.  Only  a  few  comparison  test  runs  were  taken  with  the  record 
length  set  at  8  samples;  unless  otherwise  noted,  it  should  be  assumed  that  the  record  length  was  set 
to  16  samples. 

The  size  of  the  FIFO  buffer  determined  how  many  bursts  were  to  be  collected  before  the 
burst  data  were  stored.  For  instance,  if  the  FIFO  buffer  was  set  at  4000,  the  BSAs  would  collect 
4000  bursts  and  then  send  them  to  memory;  4000  more  bursts  would  be  collected  and  then  sent  to 
memory,  etc.  A  difficulty  with  the  FIFO  buffer  size  was  that  the  last  group  of  bursts  collected  did 
not  get  sent  to  the  memory.  For  instance,  if  the  BSA  recorded  5234  bursts  and  the  FIFO  buffer 
was  set  to  4000,  the  last  1234  bursts  would  be  lost.  This  became  problematic  in  the  FPG  testing, 
when  the  streamlines  diverged,  resulting  in  a  lower  seeding  density  and  fewer  recorded  bursts  near 
the  wall.  In  order  to  send  as  many  bursts  as  possible  to  the  memory,  the  FIFO  buffer  needed  to  be 
reduced,  resulting  in  loss  of  a  smaller  number  of  burst  records.  However,  if  the  FIFO  buffer  was 
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set  too  small,  the  BSA  spent  too  much  time  transferring  data  to  memory,  resulting  in  fewer  total 
bursts  collected. 

J.3  Data  Collection  Techniques 

After  the  optics  were  aligned  and  the  software  was  prepared,  the  LDV  system  was  ready  to 
begin  data  acquisition.  The  data  were  collected  along  a  traverse  path  normal  to  the  wall;  the 
closest  distance  to  the  wall  at  which  data  collection  was  possible  was  limited  by  the  optics.  Aside 
from  taking  data  at  the  same  locations  as  the  data  taken  by  Miller  [26],  data  were  also  collected  at 
several  locations  in  the  flow  in  order  to  measure  the  flow  strain  rates. 

5. 3. 1  Traverse  Path  Across  the  Boundary  Layer.  The  traverse  path  ran  normal  to  the 
wall  at  the  location  =  0.715  m .  At  this  location,  the  equation  for  the  wall  position  provides  an 
equation  for  the  normal  line,  given  (in  traverse  coordinates)  as 

Zj.  =7.0169yj.  (5.1) 

The  error  associated  with  the  traverse  position  was  found  to  increase  linearly  with  each 
data  point  taken  across  the  boundary  layer  (see  Appendix  B).  In  addition,  the  total  temperature 
rose  slightly  as  the  compressors  heated  from  constant  use.  Therefore,  in  order  to  minimize  the 
error  and  maintain  conditions  as  constant  as  possible,  each  traverse  across  the  boundary  layer  was 
limited  to  collection  of  data  at  approximately  15  points.  Successive  traverses  across  the  boundary 
layer  collected  data  at  different  values  oi  yjd  in  order  to  create  a  more  robust  velocity  profile. 

5. 3. 2  Wall  Proximity  Limitations.  The  closest  position  to  the  wall  at  which  data  could  be 
collected  was  limited  by  the  configuration  of  the  optics,  as  shown  in  Figure  5.4.  Note  that  Figure 
5.4  gives  the  theoretical  minimum  distance  from  the  wall  as  0.71  mm,  based  on  an  infinitely  thin 
laser  beam  and  ideal  conditions.  In  reality,  the  minimum  distance  was  found  to  be  approximately 


5-9 


Figure  5.4;  Wall-Proximity  Limitation  for  Two-Component  Measurements 


1.13  mm  due  to  laser  beams  of  finite  thickness  and  reflections  off  the  wall.  Because  the  5 14.5  nm 
beams  were  oriented  parallel  to  the  wall,  it  was  possible  to  obtain  measurements  of  the  w-velocity 
component  down  to  a  distance  of  y  «  0.07  mm  from  the  wall. 

5.3.3  Finite-Difference  Grid.  Strain  rates  across  the  boundary  layer  were  computed  by 
collecting  data  at  several  locations  in  order  to  build  a  finite-difference  grid.  The  finite-difference 
approximations  used  were  central  differences  in  the  lateral  direction  and  backward  differences  in 
the  streamwise  direction.  A  backward  difference  was  chosen  because  of  the  termination  of  the  wall 
curvature  and  the  resulting  compression  shock  slightly  upstream  of  the  measurement  position.  The 
finite  difference  grid  used  is  shown  in  Figure  5.5. 


Figure  5.5:  Finite-Difference  Grid 
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The  data  were  collected  at  the  points  =71.5  cm,  =  IQ. 96  cm  and  x.^^  =  69.43  cm. 
Converting  these  into  the  body-intrinsic  coordinate  system  led  to  x  =  71 .52  cm,  x  =  70.98  cm,  and  x 
=  70.44  cm,  respectively. 

5.3.4  Atomizer  Pressure.  A  pressure-control  valve  on  the  atomizer  and  one  on  the  air 
supply  to  the  atomizer  allowed  for  fine-tuning  of  the  atomizer  pressure.  In  order  to  reduce  the 
seeding  effects,  the  pressure  was  kept  as  low  as  possible  to  reduce  the  amount  of  seed  in  the 
boundary  layer  but  still  allow  sufficient  data  collection.  The  lowest  acceptable  pressure  was  found 
to  be  about  3.035  x  10^  N/m^  .  The  effects  of  seeding  density  on  the  flow  were  investigated 
experimentally  by  collecting  data  with  the  pressure  set  at  different  levels;  the  results  from  these 
tests  are  presented  in  Appendix  A. 

5.4  Coincidence  Filtering 

Coincidence  filtering  of  the  data  was  performed  by  the  BSAs,  based  on  a  value  for  the 
coincidence  window  which  was  set  by  the  user.  According  to  Ref  [11],  best  results  would  be 
obtained  when  the  coincidence  window  was  set  at  least  twice  as  large  as  the  record  interval. 
Recalling  fi-om  Section  5. 2. 3.2  that  the  fastest  particles  crossed  the  control  volume  in 
approximately  0.4  ps,  compliance  with  this  advice  restricted  the  record  length  to  8  samples. 
However,  when  the  record  length  was  set  to  16  samples,  the  record  interval  was  0.333  ps,  slightly 
lower  than  the  minimum  traverse  time.  Following  the  recommendations  of  Ref  [11]  would  have 
presented  possible  biasing  problems,  as  the  coincidence  window  would  have  been  set  to  0.666  ps. 
This  setting  would  provide  a  high  probability  that  a  particle  inside  the  control  volume  would  be 
correlated  with  one  already  outside  the  control  volume.  To  eliminate  this  possibility,  when  the 
record  length  was  set  to  16  samples  the  coincidence  window  was  set  equal  to  the  record  interval. 
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0.333  (AS.  Thus,  the  coincidence  window  was  always  set  to  0.333  |as,  independent  of  the  choice  for 
record  length.  The  effects  of  coincidence  filtering  on  the  collected  data  were  examined  in  detail; 
the  results  of  this  investigation  are  presented  in  Appendix  A. 

5.5  Flow  Visualization 

Flow  visualization  was  accomplished  using  ten-nanosecond  shadowgraph  and  schlieren 
photography.  The  schlieren  was  sensitive  to  the  gradients  in  the  density,  while  the  shadowgraph 
was  sensitive  to  the  second-derivative  of  the  density.  The  Kolmogoroff  turbulent  length  scale  for 

this  flow  was  given  in  Chapter  1  as  2  x  10~’  m.  In  the  highly-turbulent  regions  of  the  boundary 
layer  the  flow  was  moving  at  approximately  500  m/s;  this  gave  the  turbulent  time  scale  as 
approximately  0.4  ns.  Therefore,  the  ten-nanosecond  spark  generator  was  not  expected  to  capture 
the  smallest  eddies  in  the  flow.  However,  it  had  the  shortest  spark  duration  of  any  available  source 
and  proved  adequate  for  capturing  density  gradients  across  the  boundary  layer. 

5.5.1  Shadowgraph  Alignment.  The  shadowgraph  optic  alignment  is  shown  in  Figure 
5.6.  The  light  source  was  placed  slightly  off  the  z-axis  and  at  the  focal  point  of  the  parabolic 
mirror.  With  proper  alignment,  this  produced  a  14.4  cm  collimated  light  beam  that  traversed  the 
test  section  along  the  z-axis.  Proper  alignment  was  achieved  using  a  sighting  laser:  the  laser  was 
placed  at  the  same  height  as  the  test  section,  and  the  laser  beam  was  reflected  off  the  focal  mirror 
and  projected  back  upon  itself  The  mirror  was  then  turned  horizontally  until  the  light  projected 
from  the  source  was  focused  on  the  tunnel.  The  camera  was  placed  on  the  opposite  side  of  the  test 
section  from  the  light  source,  as  close  to  the  test  section  as  possible. 

5.5.2  Schlieren  Alignment.  The  schlieren  optical  setup  was  similar  to  the  shadowgraph 
setup  as  far  as  alignment  of  the  light  source  and  first  focal  mirror;  however,  as  shown  in  Figure 
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Figure  5 . 6 :  Shadowgraph  Alignment 


5.7,  the  schlieren  setup  was  more  complex.  Once  passing  through  the  test  section,  the  light  was 
reflected  off  a  second  focal  mirror  in  the  direction  of  a  knife  edge  which  was  placed  at  the  focal 
point  of  the  second  mirror.  The  light  was  then  reflected  off  a  flat  mirror,  through  a  focusing  lens, 
and  to  the  camera. 

In  order  to  produce  quality  schlieren  photographs,  there  were  three  variables  that  needed  to 
be  fine-tuned:  the  knife  edge  needed  to  be  placed  at  the  focal  point  of  the  mirror,  the  amount  of 
light  blocked  by  the  knife  edge  needed  to  be  precisely  adjusted,  and  the  camera  had  to  be  placed  at 


Focal  Mirror 


Figure  5.7:  Schlieren  Alignment 
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the  correct  location  to  ensure  the  schlieren  was  in  focus. 


Placing  the  knife  edge  at  the  focal  point  of  the  mirror  was  relatively  simple;  the  spark 
generator  was  placed  on  “continuous”  mode  and  the  position  of  the  knife  edge  was  adjusted  until 
the  image  of  the  spark  on  the  knife  edge  came  into  focus.  The  camera  was  focused  using  a 
transparency  with  lettering  on  it  which  was  placed  on  either  side  of  the  test  section.  The  position  of 
the  camera  was  adjusted  until  the  lettering  on  the  transparency  became  clear,  and  its  position  was 
marked.  The  camera  was  focused  once  with  the  transparency  in  front  of  the  test  section  and  once 
behind;  the  position  of  the  camera  was  marked  at  each  location,  and  the  final  position  was  set  half 
way  between  the  two  in  order  to  place  the  focal  point  in  the  center  of  the  test  section.  Slight 
adjustments  were  made  during  the  course  of  taking  pictures  in  order  to  fine-tune  the  focus. 

The  amount  of  light  blocked  by  the  knife  edge  was  the  hardest  variable  to  adjust.  Blocking 
too  much  light  would  (obviously)  cause  the  schlieren  to  be  too  dark;  conversely,  if  not  enough  light 
was  blocked,  the  effect  of  the  knife  edge  would  be  lost  and  the  resulting  photograph  would  look  like 
a  shadowgraph.  The  best  schlieren  photographs  seemed  to  be  the  ones  where  1/3  to  1/2  of  the  light 
was  blocked.  The  process  of  taking  schlierens  involved  trial-and  error,  attempting  to  align  the 
lighting  and  the  focus  at  the  same  time. 

Schlieren  photographs  were  taken  with  the  knife  edge  in  two  separate  orientations.  One 
orientation  had  the  knife  edge  horizontal  in  order  to  capture  the  dpjdy  density  gradients.  The 
second  orientation  had  the  knife  edge  vertical  to  capture  the  streamwise  density  gradients. 
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6.  Data  Reduction  Procedures 


6. 1  Velocity  Calculation 

The  velocity  was  calculated  using  a  FFT  to  compute  the  Doppler  frequency  as  given  in 
Chapter  3;  the  Doppler  frequency  was  then  converted  into  a  velocity  using  the  frequency-velocity 
calibration  factor,  Cj^ .  For  the  optical  setup  in  the  current  experiment,  the  calibration  factors  for 

the  u  and  v  components  of  velocity  were 


C.„=  8.128-^,  =  7.709-^ 

MHz'  MHz 


(6.1) 


Once  the  instantaneous  velocities  had  been  measured,  the  mean  and  fluctuating  velocities 
(m  ,  V  ,  u' ,  v' )  were  computed,  allowing  calculation  of  the  total  velocity,  Q.  However,  since  the 
flow  was  turbulent,  the  mean  total  velocity  had  to  be  separated  from  the  fluctuating  total  velocity. 
This  was  accomplished  by  taking  the  Reynolds  average  of  the  total  velocity: 


=  u^  +v^ 


{Q  +  Q'}  ={u+u’f  +  (v  +  v')^ 


(6.2) 


1  +  ^  =«^fl  +  ^l  +v^ +2vv'  +  (v') 

\  0)  \  uJ 


\  Q) 

Assuming  that  u'Ju  was  small,  the  binomial  expansion  could  be  applied  to  the  first  term;  taking 
the  Reynolds  average  of  the  equation  then  led  to 

Subtracting  this  result  from  Eq.  (6.2),  the  fluctuating  velocity  was  given  by 

uu'  vv'  (6.4) 

0'  =  -^  +  -=^ 

0  e 

The  data  output  by  the  BSAs  did  not  include  calculations  of  the  instantaneous  value  of 
Q' ;  the  only  information  output  were  the  root-mean  square  (RMS)  values  of  the  velocity 
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components.  Therefore,  the  fluctuating  total  velocity  had  to  be  presented  as  an  RMS  value  also. 
Taking  the  square  of  Eq.  (6.4)  and  noting  that  ujQ  s  1  and  v/g  ^  1, 


Therefore,  in  order  to  calculate  the  fluctuating  Mach  number,  it  became  necessary  to  calculate 
a' j  a  .  Manipulation  of  the  speed  of  sound  relation  led  to 
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(6.11) 


a+a'  =  ^rK{T  +  T')=^ 


1/2 

=  a 

1+^^' 

L  T  \ 

_  2T 

where  the  last  simplification  was  made  possible  by  assuming  T’jT  was  small  and  using  the 
binomial  expansion.  Subtracting  a  from  both  sides. 


a'  1  r 


+  0 


(6.12) 


a  2T 

Where  the  symbol  O  represents  “terms  of  the  order  of  magnitude  [of  the  bracketed  terms].” 

Using  the  assumption  of  small  fluctuating  Mach  number  and  the  Strong  Reynolds  Analogy 
(SRA),  Morkovin  [28]  derived  an  expression  for  the  value  of  T'jT  : 


^=:-(r-i)M^ 


Q 


(6.13) 


Q) 


The  SRA  assumes  that  the  ratio  of  the  fluctuating  total  temperature  to  the  total 
temperature,  Tq/TJ,  ,  is  zero.  This  assumption  was  verified  by  Miller,  who  found  this  ratio  to  be 
less  than  0.02  for  the  flat  plate  case.  [26] 

For  Mach  2.8  flow  over  a  flat  plate,  values  of  Q'  jQ  ^  0.1  could  be  expected.  [16,  22,  26] 
In  the  FPG  case,  the  stabilizing  effect  of  the  FPG  was  expected  to  further  reduce  the  fluctuations. 
Therefore,  Eq.  (6. 13)  could  be  used  to  provide  an  estimate  for  the  largest  values  of  T'jT  : 


r 


<  (1.4 -1)(2.8)^  (0.1)  =  0.31 


(6.14) 


Inserting  Eq.  (6.12)  into  Eq.  (6.10)  and  retaining  only  first-order  terms,  the  expression  for 
the  fluctuating  Mach  number  became 


M'  =  M 


£ 

Q 


l  +  |(y-l)M^ 


(6.15) 


Thus,  the  RMS  value  of  the  fluctuating  Mach  number  was 
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(6.16) 


{M')  =M 


Q 


6. 3  Mach  Number  Calculation 

The  local  Mach  number  was  calculated  using  the  measured  velocity  and  total  temperature 
in  conjunction  with  the  assumptions  of  adiabatic  flow  and  a  thermally  and  calorically  perfect  gas. 
The  adiabatic  energy  equation  was  derived  for  turbulent  flow  by  expanding  the  adiabatic 
relationship  between  temperature,  stagnation  temperature  and  Mach  number: 


T  2 

_  T  +T' 

T  +  T'  = _ ^ ^ — 


T 


1  +  ^~(m  +  M')^ 
To 

1  +  + 


(6.17) 


Where  Tq  was  neglected  based  on  the  SRA.  Again  using  the  binomial  expansion  and  retaining 
only  the  first-order  terms,  Eq.  (6.17)  could  be  written  as 


=  1  +  ^^  +2MM'+{M'f 


(6.18) 


Taking  the  time  average  of  both  sides,  the  turbulent  adiabatic  energy  equation  was  found 


to  be 


4^  =  1  +  ^^M^  «1  +  — 

T  2  2  ^  2 


(6.19) 


The  above  equation  was  inserted  into  Eq.  (6.9);  the  resulting  equation  was  manipulated  to  yield 

iV2  (6.20) 


M  = 


2Q- 


2yR,T,-Q\Y-l), 
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It  should  be  noted  that  the  Mach  number  was  calculated  assuming  constant  total 


temperature  (adiabatic  flow)  throughout  the  boundary  layer.  This  assumption  was  verified  for  the 
flat  plate  case  by  measurements  made  by  Kistler.  [22]  He  reported  values  of 

T  ~T 

■^0,00  ^aw 

where  is  the  adiabatic  wall  temperature.  The  adiabatic  wall  temperature  is  calculated  using 
the  adiabatic  recovery  factor,  r  »  ,  and  the  relationship 

(6.22) 

2o,«>  ~  2’g 

where  is  the  temperature  at  the  edge  of  the  boundary  layer.  Manipulating  Eq.  (6.21), 

while  Eq.  (6.22)  becomes 

T,.  1  +  #-1)/2K 

For  Pr  =  0.71 ,  /•  =  0.892 ;  with  =  2.9,  =  0.932  and 


—  >  0.96  «  const 


(6.25) 


^  0,oo 


and  the  Mach  number  profile  can  be  calculated.  The  result  given  in  Eq.  (6.25)  also  agrees  with 
results  presented  by  van  Driest.  [35]  In  addition.  Miller  et  al.  [27]  found  that  7J,/7J)co  was 

essentially  constant  across  the  flat  plate  boundary  layer  in  the  facilities  used  for  the  present 
experiment. 
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6. 4  Density  and  Temperature  Calculation 

Calculation  of  the  density  and  temperature  profiles  were  complicated  by  the  fact  that  the 
only  measurements  taken  in  the  test  section  were  the  wall  temperature  and  pressure.  The  density 
and  temperature  in  the  boundary  layer  were  estimated  based  on  certain  assumptions;  if  the  validity 
of  these  assumptions  broke  down,  the  reported  values  of  the  temperature  and  density  (and  thus  the 
incompressible  Reynolds  shear)  would  be  inaccurate.  Thus,  in  order  to  minimize  the  error 
associated  with  the  calculations,  different  assumptions  and  methods  of  calculation  were  used  for 
the  flat  plate  and  FPG  cases. 


6.4.1  Flat  Plate.  In  the  flat  plate  region,  the  density  was  calculated  assuming  the 
pressure  normal  to  the  surface  was  constant;  with  this  assumption,  the  modified  Crocco-Busemann 
approximation  yields  [35,  6] 


■^  =  — =  l  +  fi'(— ^ 


P 


-{Ai 


(6.26) 


(6.27) 


where  and  are  the  density  and  temperature  at  the  wall  and  A'  and  B'  are  given  by 

V  ;  j,  jj 

J,,  l  +  [fr-l)/2K  1 

T  IT 

^ e 


The  density  at  the  wall,  p^,  was  calculated  using  the  measured  wall  temperature  and 
pressure  and  the  perfect  gas  relation: 


KT. 


(6.28) 


6.4.2  FPG  Region.  The  assumption  of  constant  pressure  across  the  boundary  layer  was 
not  necessarily  valid  in  the  FPG  test  section  and  needed  to  be  evaluated.  In  the  presence  of  a 
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curved  body,  the  streamlines  will  experience  curvature,  causing  a  pressure  gradient  normal  to  the 
streamlines.  The  magnitude  of  this  gradient  was  given  by  Spina  et  al  as  [33] 


^  -piP'  _  (6-29) 

R  R 


and  the  pressure  difference  between  the  wall  and  the  outer  edge  of  the  boundary  layer  could  be 
estimated  as 


(6.30) 


If  the  ratio  SjR  was  made  small  enough  (usually  by  a  large  radius  of  curvature  and 
therefore  a  small  streamwise  pressure  gradient,  approaching  the  flat  plate  case)  the  pressure 
difference  could  be  made  negligible.  For  the  present  test  configuration,  the  radius  of  curvature  at 
the  test  location  was  48.08  cm  and  S  « 1.28  cm;  with  =  2.9 , 


Pe 


=  -0.31 


(6.31) 


Equation  (6.31)  shows  that  the  assumption  of  constant  pressure  across  the  boundary  layer 
was  not  valid  in  the  FPG  test  section.  In  order  to  avoid  large  errors  in  density  calculations  caused 
by  making  this  assumption,  an  alternate  method  for  calculating  the  temperature  and  density  was 
used. 


The  alternate  method  assumed  that  the  pressure  distribution  between  the  wall  and  the 
freestream  was  linear.  A  numerical  investigation  of  the  flowfield  in  the  current  experimental 
configuration  {cf.  Ref  [17])  found  that  the  pressure  variation  across  the  boundary  layer  was  nearly 
linear,  justifying  this  assumption.  In  order  to  perform  the  linear  interpolation,  the  wall  and  edge 
pressures  were  required.  The  wall  pressure  had  been  measured;  the  edge  pressure  was  calculated 
using  the  assumption  that  there  were  no  shocks  present  in  the  freestream,  leading  to  a  constant  total 
pressure  along  the  edge  of  the  boundary  layer  and 
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(6.32) 


Pe  =Po 


~r 

r-\ 


The  pressure  at  any  point  in  the  boundary  layer  could  then  be  estimated  using  linear 
interpolation.  The  temperature  across  the  boundary  layer  was  found  assuming  constant  total 
temperature  and  using  Eq.  (6.19);  with  the  temperature  and  pressure  known,  the  density  was 
calculated  using  the  perfect  gas  relation. 

It  is  important  to  mention  here  that  this  analysis  was  not  expected  to  provide  an  exact 
measure  of  the  density,  temperature  and  pressure  across  the  boundary  layer.  The  sole  purpose  was 
to  obtain  an  estimation  of  the  density  across  the  boundary  layer  in  order  to  present  the  direct 
measurements  of  u’V  in  a  manner  which  allowed  comparison  to  data  in  the  literature;  the  form 
most  commonly  used  was  -^pu'v'jr^  ,  necessitating  the  estimation  of  the  density. 


6.5  Boundary  Layer,  Displacement  and  Momentum  Thickness  Calculations 

Two  boundary  layer  thicknesses  were  calculated;  one  thickness  based  on  the  velocity  ( ) 

and  one  on  the  Mach  number  It  was  necessary  to  calculate  both  thicknesses  for  two 

reasons.  First,  in  order  to  accurately  compare  the  present  results  with  the  Mach  number 
measurements  of  Miller  et  al.  [27]  which  were  based  on  conventional  Pitot  and  cone-static 
pressure  profiles,  needed  to  be  calculated.  The  second  reason  for  calculating  both  thicknesses 
was  that  the  Mach-number  boundary  layer  thickness  gave  an  estimation  of  the  temperature 
boundary  layer  thickness,  the  thickness  which  was  measured  from  the  flow  visualization  results. 

The  boundary  layer  thicknesses  were  calculated  by  choosing  the  largest  velocity  or  Mach 
number  from  a  given  traverse  as  the  edge  value,  non-dimensionalizing  the  profile  by  the  edge  value, 
and  searching  outward  from  the  wall  until  a  value  of  u/u^  or  M/M^  >  0.995  was  attained.  The 
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value  of  5  was  then  calculated  by  linear  interpolation  between  the  data  points  just  above  and 


below  the  boundary  layer  edge. 

Integral  analysis  of  the  boundary  layer  equations  leads  to  the  definition  of  the  displacement 
and  momentum  thickness,  [4]  The  displacement  and  momentum  thicknesses  are  given  by 

(6.33) 


•’A  pjiJ 


© 


=  f  ^ 

*  pji„ 


Pe^eJ 
1-- 


dy 


respectively.  The  integration  was  carried  out  numerically  using  the  measured  data  and  a 
trapezoidal  integration  scheme.  In  order  to  filter  out  the  uncertainty  associated  with  the  density 
calculations,  “incompressible”  (or  kinematic)  displacement  and  momentum  thicknesses  were  also 
calculated: 

(6.34) 


A  ' 


dy 


0,..J 


1  — 

UJ 


dy 


Note  that  the  kinematic  thicknesses  were  based  solely  on  measured  data,  reducing  the  possibility  of 
error  incurred  from  inaccurate  assumptions.  Also  note  the  asymptotic  nature  of  the  integrals 
reduced  the  error  associated  with  random  fluctuations  which  may  have  affected  the  choice  of 


and  Sf^ ;  as  a  result,  these  parameters  were  less  subject  to  random  errors. 


6. 6  Turbulence  Statistics 

Calculation  of  statistical  variables  associated  with  the  turbulent  flow  was  important  for 
two  reasons.  The  first  was  that  some  of  the  statistical  variables  appeared  directly  in  the  time- 
averaged  forms  of  the  Navier-Stokes  equations.  The  second  was  that  measurement  of  the  turbulent 
statistics  helped  to  describe  various  features  of  the  turbulent  flowfield. 
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As  expected,  the  main  component  around  which  all  the  turbulent  statistics  were  developed 
is  the  fluctuation  from  the  mean,  x'  =  x-x .  The  variable  x  was  used  here  to  denote  the  fact  that 
the  statistical  calculations  were  performed  for  both  the  u  and  v  components  of  velocity.  Since  x' 
was  centered  about  the  mean,  its  average  value  was  zero  by  definition.  Thus,  in  order  to  get  an 
estimate  of  the  average  size  of  the  fluctuation  from  the  mean,  the  RMS  value  of  x'  was  employed. 
The  RMS  value  (also  known  as  the  standard  deviation)  of  x'  is  defined  as 


Xrms  =  o-x  = 

where  the  summation  is  carried  out  over  the  «  samples  collected.  It  should  be  noted  that,  to  be 
statistically  correct,  the  denominator  of  the  summation  should  be  «  -  1  as  shown  above.  However, 
the  Burstware®  software  performed  these  calculations  based  on  a  denominator  of  n .  Because  of 
the  large  number  of  samples  collected  («>2000  for  all  data  points),  the  difference  between  the 
two  was  negligible  (err <0.05%).  The  skewness  (Sk,^)  and  flatness  (Fl,^ )  were  also  calculated 
by  the  Burstware®  as 


Sk, 


aln 


Fix- 


4 

(7M 


(6.36) 


Once  again,  these  were  not  the  statistically  unbiased  definitions,  but  the  large  number  of  samples 
collected  negated  the  error  caused  by  the  above  formulations. 

The  skewness  function  is  important  to  LDV  measurements  in  helping  to  determine  whether 
or  not  the  results  were  biased  by  seeding  problems.  The  flatness  can  be  used  to  calculate  the 
intermittency  function,  .  The  intermittency  function  is  the  ratio  between  the  flatness  of  the 
Gaussian  curve  (3.0)  and  the  flatness  of  the  measurements  (i.e.,  =  3.0/Fl„ )  and  provides  a 

measure  of  the  fraction  of  time  that  the  flow  is  turbulent  at  a  given  location. 
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The  distribution  of  the  intermittency  function  was  compared  with  the  Gaussian  integral 
curve  given  by 

y^=05i\-Qrix)  (6-37) 

where 


-1 

V 

J 

(6.38) 


and 


2  ^  2 

edX  =  ^\e'*^dt 


(6.39) 


The  constants  ^  and  (p  were  adjusted  until  the  Gaussian  curve  matched  the  measured  distribution; 
“such  a  curve  shows  that  the  instantaneous  position  of  the  edge  of  the  boundary  layer  has  a  random 
character  with  a  mean  position  at  y/S^  =[<'>]  •”  [Klebanoff]  The  constant  ^  represented  the 
standard  deviation  of  the  boundary  layer  edge  from  . 

The  next  statistical  variable  to  be  calculated  was  arguably  the  most  important.  The 
velocity  correlation  u'v'  appeared  directly  in  the  RANS  and,  along  with  the  density,  comprised  the 
incompressible  Reynolds  shear  stress.  The  correlation  was  calculated  using  the  formula 

u'v'  =  (6.40) 

Note  that,  while  the  average  values  of  u'  and  v'  were  zero,  the  average  of  the  product  u'v'  may 
not  be  zero. 

The  final  statistical  value  calculated  was  the  velocity  correlation  coefficient,  ,  given  by 


u'v' 


O^u^v 

was  used  to  provide  insight  into  the  nature  of  the  turbulence, 
turbulence  was  isotropic  (no  preferred  direction  of  turbulent  structures),  R^ 


(6.41) 

For  instance,  if  the 
should  have  had  a 
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magnitude  of  1 .0.  [6]  The  Cauchy-Shwartz  inequality  indicated  that  <  1.0 .  [3]  Finally,  low- 
speed  flat  plate  data  indicated  that  =  0.4  -  0.5 .  [6] 


6. 7  Wall  Shear  Stress 

The  wall  shear  stress  was  calculated  using  the  bisection  method  to  iteratively  solve  the 
skin-friction  relationship  developed  by  van  Driest  [35]; 


0.242 

2{A')  -B' 

+  sin  ' 

(  > 

B' 

A'^Cf{TjT^) 

[^4(Af +{B’y  \ 

=  AT  +  log(Re^  C^)  -  <ylog(7;/7;) 


where  /r  =  0.41  and  (o  =  0.68  .  It  should  be  noted  that  this  relation  is  for  turbulent  flow  over  a  flat 
plate;  its  usefulness  in  computing  the  wall  shear  in  the  FPG  test  section  will  be  determined  later. 
It  is  also  important  to  point  out  that  Re^  represents  the  Reynolds  number  based  on  the  distance 
from  the  leading  edge  of  a  flat  plate  immersed  in  a  uniform  flowfield.  In  order  to  estimate  a  value 
for  Re^ ,  the  wind  tunnel  was  modeled  as  a  flat  plate  in  uniform  flow  where  the  throat  of  the  wind 
tunnel  was  taken  as  the  leading  edge  of  the  flat  plate.  The  inaccuracy  of  this  assumption  and  its 
effect  on  the  reported  data  will  also  be  discussed  later. 


6. 8  Van  Driest  Velocity  Profile 

In  order  to  compare  the  velocity  profiles  in  compressible  flows,  van  Driest  [35]  proposed 
transformation  of  the  velocity  profile  into  an  effective  velocity  given  by 


^eff  - 


(6.43) 
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The  values  of  pjp^,  were  calculated  using  the  Crocco-Busemann  relation  given  in  Eq.  (6.26); 
using  this  relation  in  Eq.  (6.43)  allows  the  effective  velocity  to  he  calculated  analytically. 
Performing  the  integration  leads  to 


sin  ^ 

2{A')\ulu,)-B' 

+  sin  * 

f  \ 

B' 

[  ^4(Af  i-{Bf  j 

+(»')’ J. 

(6.44) 


Van  Driest  also  suggested  non-dimensionalizing  the  effective  velocity  by  the  wall-friction 
velocity,  u  ,  where 

(6.45) 


The  non-dimensional  velocity  ju*  was  plotted  versus  a  scaled y-coordinate. 


(6.46) 


where  is  the  fluid  viscosity  measured  at  wall  temperature.  This  analysis  allowed  the  empirical 
relationship  between  and  u/ii*  which  has  been  developed  for  low-speed,  incompressible  flat 
plate  flow,  to  be  extended  to  high-speed,  compressible  flat  plate  flow.  [36]  Thus, 


^eff  1 ,  +  2n  .  2r^y 

=  — Iny  +C  + - sin^'  — — 

UK  K  \2  6 


(6.47) 


where  /c  =  0.41  and  C  =  5.0 .  n  is  the  Coles  wake  parameter,  given  by  White  [36]  as 

n«0.8[y9  +  0.5f''  (6-48) 

The  variable  P  in  Eq.  (6.48)  is  Clauser’s  equilibrium  parameter,  defined  as 

S*  dp^ 


(6.49) 


dx 


It  should  he  noted  that  Eq.  (6.47)  is  valid  in  the  “logarithmic”  region  (80<y^  <500)  as  well  as 
the  outer  portion  of  the  boundary  layer.  However,  for  flow  very  near  the  wall  (y^  <  10 )  the  flow 
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is  moving  slow  enough  that  it  is  laminar.  This  region  is  referred  to  as  the  laminar  sublayer,  in  the 
laminar  sublayer,  the  velocity  correlation  is  given  zs  uju*  -y^  .  [6]  For  the  current  experiment, 
the  laminar  sublayer  was  only  (approximately)  5  pm  thick;  due  to  the  wall-proximity  limitations 
on  the  LDV  system,  data  collection  was  not  possible  in  the  laminar  sublayer. 


6. 9  Turbulent  Kinetic  Energy 

The  turbulent  kinetic  energy  (TKE)  provides  a  measure  of  the  energy  contained  in  the 
turbulent  motion  of  the  flow.  The  TKE  is  defined  as 


2 

For  a  two-dimensional  flow,  even  though  w  =  0,  the  fluctuation  in  the  lateral  direction,  w' ,  will 
not  necessarily  be  zero.  The  current  experimental  setup  did  not  permit  measurements  of  w' ;  in 
order  to  calculate  the  TKE,  an  assumption  was  made  to  estimate  its  size.  It  was  assumed  that  the 
fluctuation  w'  is  of  the  same  order  as  v' ,  and  the  estimated  TKE  was  calculated  as 


2  2 


+ 


(v')^ 


(6.51) 


6.10  Strain  Rate  Calculation 

As  mentioned  in  Chapter  5,  data  were  collected  at  the  points  in  a  finite-difference  grid. 
Using  the  notation  shown  in  Figure  6.1,  the  central-derivatives  in  the  y-  and  z-directions  at  the 
point  {i,j)  were  estimated  as  [1] 


du 

dy 

ij,k  yij^\,k  yij-\,k 

du 

dz 

i,J,k  “  ^i,j,k-\ 

(6.52) 
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y 


Figure  6.1:  Finite-Difference  Grid  Notation 


The  backward-difference  approximation  was  complicated  by  the  fact  that  the  data  were  not 
taken  at  evenly-spaced  locations.  Using  a  Taylor’s  series,  the  velocity  at  locations  ,  Xf^_2 
was  given  as 


_ /  \  _ /  \  du  /  \  d  u 

«(x_,)  =  «(xo)  +  —  (x_,-Xo)  +  ^ 


du\ 


d^u 


U{X_2)  =  «(Xo)  +  —  {X_2  -  Xo)  +  ^^2 

Xq 


{X-2-Xof 


(6.53) 


+ 


+ 


'*0 


where  the  subscripts  were  understood  to  represent  the  value  of  i.  Multiplying  the  top  equation  by 
-(x_2  -  XqY  ,  the  lower  equation  by  (x__,  -  x^)^ ,  adding,  and  solving  for  the  derivative  at  Xg , 

(6.54) 


du 

dx^^u(x_2)- 

<fX2^?/(x_j)  + 

dxY 

-dx^ 

M(Xg) 

dx 

Xo 

dx2dx^  —dx-^dx-Y 

where  dx^  =  x.,  -  Xg  and  dx2  =  x_2  -  Xg . 
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6.11  Streamwise  Pressure  Gradient 

The  streamwise  pressure  gradient,  dp^Jdx,  was  needed  in  the  calculation  of  Clauser’s 
equilibrium  parameter  given  in  Eq.  (6,49).  The  pressure  gradient  was  calculated  assuming  the 
streamwise  rate  of  change  of  pressure  was  constant  across  the  boundary  layer,  i.e., 

^^dp^  (6.55) 

dx  dx 

With  this  assumption,  the  measured  values  of  wall  pressure  were  inserted  into  Eq.  (6.54) 
in  place  of  u  ,  resulting  in  a  value  for  dp^jdx . 
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7.  Results  and  Discussion 


Recall  that  the  underlying  goal  of  the  present  study  was  to  use  advanced  laser  Doppler 
velocimeter  techniques  to  investigate  the  turbulent  structures  in  a  supersonic  boundary  layer  and 
the  effect  of  a  pressure  gradient  on  those  structures.  In  order  to  present  and  analyze  the  measured 
LDV  data  in  a  concise  and  methodical  manner,  it  is  first  necessary  to  discuss  some  of  the 
supporting  data,  such  as  measurements  of  the  wind  tunnel  operating  conditions  and  the  flow 
visualization  results.  Sections  7. 1  and  7.2  discuss  these  supporting  measurements;  the  remaining 
sections  of  this  chapter  present  and  analyze  the  flat  plate  and  FPG  region  results.  Since  the  present 
LDV  system  had  never  been  operated  in  the  facilities  used  for  this  experiment,  several  tests  were 
performed  to  determine  the  susceptibility  of  the  measurements  to  such  variables  as  the  record 
length,  seeding  density,  and  coincidence  window.  The  results  from  these  tests  are  presented  in 
Appendix  A.  A  detailed  error  analysis  for  all  the  measured  and  computed  variables  is  presented  in 
Appendix  B.  Finally,  to  be  useful  as  a  basis  for  the  development  of  turbulence  models,  the  criteria 
of  Settles  and  Dodson  [31]  (see  Chapter  1)  needed  to  be  met.  Any  elements  of  the  criteria  which 
had  not  been  previously  addressed  (data  repeatability,  flow  two-dimensionality,  and  tabulated  data) 
are  presented  in  Appendix  C. 

7. 1  Pressure  and  Temperature  Measurements 

Measurements  of  the  stagnation  pressure  showed  that  the  wind  tunnel  was  operated  at  a 

stagnation  pressure  of  2.13  x  10^  ±  1.38  x  10^  Pa.  The  average  stagnation  temperature  was 
approximately  298  K  with  a  +0.1  K  variation  during  each  run.  The  mean  value  was  found  to 
increase  slightly  as  testing  progressed  due  to  constant  operation  of  the  compressors.  The  mean 
temperature  never  rose  above  303  K.  Because  of  this  variation  in  temperature,  the  stagnation 
temperature  and  pressure  were  recorded  for  each  data  point  and  the  data  was  reduced  accordingly. 
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The  results  from  the  wall  pressure  measurements  in  the  FPG  region  are  summarized  in 
Table  7.1.  Note  from  Eq.  (4.2)  that  the  location  x  =  63.22  cm  was  still  in  the  flat  plate  region, 
upstream  of  the  wall  curvature.  The  pressure  change  from  this  station  to  the  location  x  =  71 .52  cm 
was  negligible  in  the  flat  plate  case,  allowing  the  pressure  measurement  at  x  =  63.22  cm  to  be  used 
for  the  flat  plate  wall  pressure  as  well.  Using  the  data  from  Table  7.1  and  Eq.  (6.54),  the 
streamwise  derivative  of  the  wall  pressure,  dp^/dx,  was  estimated.  Assuming  the  rate  of  change 
of  the  edge  pressure  was  equal  to  the  rate  of  change  of  wall  pressure  (this  assumption  was 
qualitatively  verified  by  comparison  to  CFD  results  presented  in  Ref  [17]  and  shown  in  Figure 
7. 1),  the  pressure  gradient  in  the  FPG  region  was 

(7.1) 

dx  dx  m 
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Table  7.1:  Wall  Pressure  Measurements  In  the  FPG  Test  Seetion 


X  (cm) 

63.22 

68.92 

71.52 

P.  (Pa) 

8690 

6000 

5310 

Wall  temperature  measurements  were  restricted  to  the  flat  plate  and  produced  a  result  of 
=  295  K .  For  Mach  2.8  flow  over  a  flat  plate,  the  ratio  of  adiabatic  wall  temperature  to  total 
temperature  was  given  in  Chapter  6  as  =  0.93 .  With  7J, „  « 298  278  K  (a  6% 

difference  from  the  measured  value  of  ).  Therefore,  if  necessary,  the  assumption  of  an  adiabatic 
wall  may  be  justified. 

7.2  Flow  Visualization 

The  FPG  shadowgraph  and  schlieren  results  are  presented  in  Figures  7.2  to  7.4.  Figure 
7.2  shows  the  shadowgraph.  Figure  7.3  displays  the  schlieren  which  emphasizes  the  djdy  gradient 
of  density,  and  Figure  7.4  presents  the  schlieren  which  captures  the  djdx  gradient  of  density.  In 
all  cases,  the  airflow  was  moving  from  left  to  right.  Examination  of  Figures  7.2  to  7.4  showed  that 
the  seam  just  upstream  of  the  test  section  generated  two  weak  shocks  which  propagated 
downstream.  Since  these  waves  were  very  weak,  their  effects  on  the  flow  were  deemed  negligible. 
The  strong  shock  near  the  wall  on  the  right  side  of  the  photographs  was  due  to  the  fact  that  the  wall 
curvature  terminated  at  x^^  =  12.7  cm. 

Although  the  smallest  eddies  appeared  smeared,  as  expected  (see  Section  5.5),  Figure  7.3 
was  still  able  to  capture  the  ffpidy  gradients  in  the  flowfield,  clearly  showing  the  edge  of  the 
boundary  layer.  Since  the  schlieren  provided  density  gradient  information,  the  boundary  layer 
evident  in  the  schlieren  (and  shadowgraph)  photographs  was  5j,  the  thermal  boundary  layer 
thickness.  The  boundary  layer  thickness  estimated  from  the  photographs  was  5^.  « 12.0  mm. 
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Figure  7.4:  Schlieren  (Vertical  Knife  Edge)  in  FPG  Region 
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Due  to  the  slight  unsteadiness  in  the  10  ns  light  source,  the  background  shading  varied 
from  photograph  to  photograph.  Therefore,  the  schlierens  in  Figure  7.4  were  presented  separately 
in  order  to  refrain  from  inadvertently  introducing  artificial  gradients  into  a  composite  photograph. 
The  changes  in  the  shading  in  Figure  7.4a  showed  that,  as  the  flow  sped  up  through  the  expansion 
area,  the  density  changed.  Examination  of  the  adiabatic  relationships  given  in  Eqs.  (6.19)  and 
(6.32)  showed  that,  as  the  Mach  number  increased,  the  density  decreased.  Therefore,  the  darkest 
sections  of  the  photograph  corresponded  to  the  highest  density  regions  in  that  photograph.  With 
this  knowledge.  Figure  7.4b  showed  that  the  density  began  to  increase  once  again  after  the  flow 
reached  the  end  of  the  curved  wall.  Detailed  inspection  of  Figure  7.4a  also  revealed  the  presence  of 
organized  structures  within  the  boundary  layer  region.  The  angle  these  structures  made  with  the 
test  section  wall  appeared  to  be  approximately  40-60°. 

7.3  Flat  Plate  Results 

Investigation  of  the  flat  plate  (zero-pressure  gradient,  or  ZPG)  results  was  important  for 
two  reasons:  first,  the  flat  plate  data  were  used  as  a  basis  for  comparison  against  which  the  effects 
of  the  FPG  were  assessed.  Second,  the  flat  plate  data  allowed  the  present  measurement  techniques 
to  be  evaluated  against  accepted  experimental  data  and  correlations  from  the  literature. 

7.3.1  Mean  Flow  Characteristics.  Mean  flow  characteristics  which  were  compared  to 
data  in  the  literature  included  the  velocity,  Mach  number,  and  density  profiles.  In  addition,  the 
boundary  layer  thicknesses  (S,S* ,S* ,6,6^)  were  compared  to  accepted  empirical  correlations  in 
order  to  assess  the  aceuracy  of  the  measurements. 

7. 3. 1.1  Velocity  Profiles  and  Boundary  Layer  Thickness.  In  order  to  determine 
the  boundary  layer  thickness,  the  mean  velocity  was  scaled  by  the  edge  velocity  (m^  «  599  m/s  at  x 
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=  71.5  cm)  and  plotted  versus  distance  from  the  wall  (see  Figure  7.5).  Scaling  the  velocity  by  the 
edge  value  was  necessary  because  the  wind  tunnel  conditions  changed  slightly  between  traverses 
across  the  boundary  layer,  leading  to  a  slight  fluctuation  in  the  freestream  velocity.  The  value  of 
delta  was  taken  as  the  distance  from  the  wall  at  which  uju^  =  0.995 .  From  Figure  7.5,  the 
boundary  layer  thickness  over  the  flat  plate  was  found  to  be  approximately  9.0  mm.  The  mean 
velocity  profile  over  the  flat  plate  was  plotted  versus  yjS  in  Figure  7.6  and  compared  to  Pitot  and 
cone-static  pressure  probe  data  presented  in  Ref.  [26].  It  should  be  noted  that,  in  Figure  7.6  and 
all  subsequent  presentations  of  hot-wire  and  conventional  probe  (CP)  results,  the  coordinate 

associated  with  hot-wire  and  conventional  probe  measurements  has  been  multiplied  by  the  ratio 
^mI^u  calculated  from  the  present  LDV  data  in  order  to  present  an  accurate  comparison.  For  the 
flat  plate  case  - 1.11  while  in  the  FPG  region  -  1.09  .  As  can  be  seen  in  Figure 

7.6,  the  agreement  of  the  present  results  to  those  in  Ref  [26]  was  excellent. 


Figure  7.5:  Flat  Plate  Velocity  Profile  versus  Distance  From  the  Wall 
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Figure  7.6:  Flat  Plate  Mean  Velocity  Profile 

7. 3. 1.2  Mach  Number.  The  flat  plate  Mach  number  profile  was  calculated  as 
discussed  in  Section  6.3  and  is  presented  in  Figure  7.7.  The  data  presented  in  Ref.  [26]  was  taken 
at  the  location  x  =  44  cm,  where  »  2.9,  as  opposed  to  «  2.8  measured  downstream  at  x  = 
71.5  cm  (the  reduced  fireestream  Mach  number  at  x  =  71.5  cm  was  due  to  the  continued  boundary- 
layer  growth  within  the  wind  tunnel).  Since  the  Mach  number  was  dependent  on  the  location  in  the 
wind  tunnel,  LDV  data  needed  to  be  collected  at  x  =  44  cm  in  order  to  compare  the  present  LDV 
data  to  the  conventional  data.  The  x  =  44  cm  location  was  in  the  center  of  the  first  wind  tunnel  test 
section,  hence  in  a  position  not  covered  by  the  optical  windows  (see  Figure  4.4).  In  order  to  collect 
data  at  this  location,  the  tunnel  walls  were  replaced  with  optical-grade  Plexiglas.  The  Plexiglas 
allowed  the  lasers  to  make  measurements  of  the  mean  velocity,  but  the  quality  of  the  Plexiglas  was 
not  good  enough  to  provide  for  accurate  measurements  of  the  higher-order  velocity  correlations. 
The  results  from  the  Mach  number  calculations  at  x  =  44  cm  showed  excellent  agreement  with  the 
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Figure  7.7;  Flat  Plate  Mach  Number  Profile 

conventional  data,  proving  the  assumptions  used  to  extract  the  Mach  number  from  the  LDV  data 
were  valid. 
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7.3. 1.3  Density.  For  flow  over  a  flat  plate,  the  density  profile  was  calculated 
using  the  modified  Crocco-Busemann  relation;  results  from  these  calculations  are  presented  in 
Figure  7.8.  The  assumptions  used  to  extract  the  density  from  the  LDV  data  were  the  most 
questionable  and  affected  calculation  of  other  variables  such  as  the  wall  shear  stress.  Therefore,  in 
order  to  validate  the  density  calculations,  the  density  profile  was  compared  to  results  from  the 
literature.  Specifically,  the  computed  density  was  compared  to  conventional  probe  measurements 
taken  by  Dotter  [13]  in  the  facilities  used  for  the  current  experiment  and  CFD  calculations  given 
in  Ref.  [17].  The  agreement  was  once  again  excellent;  the  discrepancy  from  the  CFD  results  near 
the  freestream  is  probably  due  to  the  fi’eestream  turbulence  in  the  wind  tunnel. 
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Figure  7.8:  Flat  Plate  Non-Dimensional  Density  Profile 


7. 3. 1.4  Integral  Momentum  Thicknesses.  The  momentum  integral  thicknesses 
were  calculated  at  each  traverse  position  across  the  boundary  layer.  The  results  from  these 
calculations  are  presented  in  Table  7.2. 


Table  7.2:  Flat  Plate  Integral  Thicknesses  Calculated  From  LDV  Measurements 


X  (cm) 

S*  (mm) 

0  (mm) 

S*  (mm) 

0j  (mm) 

71.5 

1.31 

1.27 

0.50 

0.61 

67 

1.50 

1.09 

0.46 

0.59 

To  determine  the  accuracy  of  the  measured  integral  thicknesses,  the  results  from  Table  7.2 
were  compared  with  data  from  the  literature.  The  value  of  6  was  validated  by  comparison  with 
hot-wire  results  from  Ref  [22]  and  LDV  results  from  Ref  [16];  the  hot-wire  results  presented  in 
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Ref.  [22]  indicated  that,  for  Mach  2.8  flat  plate  flow,  1.5  mm,  while  the  LDV  results  from  Ref 
[16]  measured  ^=0.84  mm.  Obviously,  the  exact  value  of  6  depends  on  the  experimental 
facility;  the  value  measured  in  the  current  experiment  falls  between  the  values  presented  in  the 
literature  and  was  accepted  as  accurate. 

7.S.2  Turbulence  Intensities.  Turbulence  intensities  were  calculated  for  both  the  u  and  v 
components  of  velocity  and  are  presented  in  Figures  7.9  and  7. 10,  respectively.  In  Figure  7.9,  note 
the  close  agreement  to  the  LDV  data  of  Ref.  [16];  the  LDV  data  of  Ref.  [30]  did  not  compare 
quite  as  well,  but  it  was  reasonable.  In  contrast,  the  hot-wire  data  of  Refs.  [26]  and  [22]  followed 
each  other  fairly  well,  but  showed  a  significant  reduction  in  the  value  of  the  turbulence  intensity 
near  the  wall  when  compared  to  the  LDV  data.  In  order  to  extract  the  turbulence  intensities  from 
the  hot-wire  data,  the  assumption  of  negligible  static  pressure  fluctuations  (i.e.,  p’  =0)  was  made. 
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Figure  7.9:  Flat  Plate  w-Turbulence  Intensity  Profile 
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The  differences  between  the  hot-wire  and  LDV  data  may  signal  possible  errors  in  this  assumption. 
Finally,  the  LDV  data  appeared  to  follow  the  hot-wire  data  of  Ref  [23]  closely.  The  data  from 
Ref  [23]  was  hot-wire  data  taken  in  low-speed,  incompressible  flow,  where  the  p'  =  Q  assumption 
was  not  necessary.  These  results  tended  to  reinforce  Morkovin’s  hypothesis  that  the  turbulent 
structures  (at  least  to  first  order)  are  relatively  unaffected  by  compressibility.  [28] 

Figure  7.10  shows  a  plot  of  the  v-turbulence  intensity  profile.  Note  the  excellent 
agreement  in  all  the  data  except  for  the  data  from  Ref.  [26].  In  fact,  the  large  perturbation  of  the 
Ref  [26]  data  from  the  other  flat  plate  data  prompted  examination  of  the  results  from  Ref  [26].  If 
both  the  turbulence  intensities  were  plotted  on  the  same  axes,  the  v-component  data  from  Ref.  [26] 
fell  nearly  on  top  of  the  LDV  ^-component  data  and  vice-versa  (see  Figure  7. 1 1),  prompting  one  to 
believe  that  the  author  of  Ref  [26]  may  have  inadvertently  switched  the  data  from  the  hot-wires. 


Figure  7.10:  Flat  Plate  v-Turbulence  Intensity  Profile 
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In  any  event,  with  the  exception  of  the  data  presented  in  Ref.  [26],  the  present  results  agreed  well 
with  accepted  data  from  the  literature  and  were  taken  as  valid. 

7. 3. 3  Fluctuating  Mach  Number.  The  fluctuating  Mach  number  profile  is  given  in  Figure 
7.12.  Note  that,  in  all  regions  of  the  boundary  layer,  the  fluctuating  Mach  number  was  much  less 
than  1.0.  This  is  important  because  if  the  fluctuating  Mach  number  approaches  unity,  local 
“shocklets”  can  form  in  the  flow.  [33]  In  addition,  Morkovin’s  hypothesis  can  also  be  stated  as 
“the  dynamics  of  a  compressible  boundary  layer  will  follow  the  incompressible  pattern  closely,  as 
long  as  the  fluctuating  Mach  number,  M' ,  remains  small.”  [33]  Since  M' «  1  for  the  current 
flowfield,  Morkovin’s  hypothesis  should  hold  for  the  flat  plate  case. 
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Figure  7. 12;  Fluctuating  Mach  Number  over  Flat  Plate 

7.3.4  Wall  Shear  Stress.  Before  continuing  examination  of  the  velocity  measurements,  it 
was  necessary  to  calculate  the  wall  shear  stress  in  order  to  determine  the  scaling  parameters  used  in 
the  presentation  of  the  incompressible  Reynolds  stress  and  in  van  Driest’s  velocity  profiles.  The 
wall  shear  stress  was  calculated  using  the  relationship  developed  by  van  Driest  given  in  Eq.  (6.42) 
and  was  found  to  have  an  average  value  of  67.0  Pa.  Recall  that  van  Driest’s  correlation  was 
developed  for  uniform  flow  over  a  flat  plate  and  that  its  use  required  the  tunnel  to  be  modeled  as  a 
flat  plate  in  Mach  2.8  flow  whose  leading  edge  was  at  the  same  location  as  the  throat  of  the  wind 
tuimel.  In  reality,  the  wind  tunnel  had  a  finite  boundary  layer  thickness  at  the  throat,  and  as  the 
flow  passed  through  the  nozzle  it  was  exposed  to  a  substantial  favorable  pressure  gradient.  In 
addition,  from  x  =  27  cm  to  the  test  location  (x  =  71.52  cm)  the  tunnel  walls  were  flat,  but  due  to 
boundary  layer  growth  the  fi'eestream  was  under  the  influence  of  a  slight  adverse  pressure  gradient. 
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Therefore,  the  assumptions  used  in  the  derivation  of  van  Driest’ s  relation  were  not  completely  valid 
for  the  current  experimental  configuration. 

7.3.5  "Incompressible”  Reynolds  Stresses.  Recall  that  the  incompressible  Reynolds 
stress  was  defined  as  -'pu'v' .  In  order  to  present  this  term  in  a  non-dimensional  fashion,  two 
methods  of  scaling  were  employed.  The  first  method  was  the  most  common,  non-dimensionalizing 
by  the  wall  shear  stress.  Figure  7.13  presents  the  LDV  data  in  this  manner  compared  to  the  data 
from  the  literature.  Note  that  the  incompressible  Reynolds  shear  stresses  measured  in  this  study 
did  not  exhibit  a  tendency  to  roll  off  near  the  wall.  The  roll-off  experienced  by  other  researchers 
has  been  attributed,  at  least  in  part,  to  angular  biasing.  [15]  Aligning  the  lasers  so  the  fringe 
pattern  ran  perpendicular  to  the  mean  direction  of  flow  seemed  to  have  eliminated  this  bias. 


-p  u'v'  /Tw 

Figure  7.13:  Incompressible  Reynolds  Shear  Stress  Scaled  by  Wall  Shear  Stress 


7-14 


Note  that  in  Figure  7.13  the  incompressible  Reynolds  shear  stress  measured  in  the  current 
experiment  tended  to  overshoot  the  wall  shear  stress.  The  overshoot  could  have  been  the  result  of 
one  of  two  errors.  The  most  likely  error  was  that  the  wall  shear  stress  was  underestimated.  A 
larger  value  of  the  wall  shear  stress  would  shift  the  LDV  data  to  the  left,  causing  the  overshoot  to 
be  decreased  or  eliminated.  The  second  possible  explanation  for  the  overshoot  was  that  the  flat 
plate  data  was  actually  taken  in  a  slightly  adverse  pressure  gradient  (APG)  as  mentioned  above. 
Studies  in  APG  flows  (such  as  that  performed  in  Ref.  [18])  have  shown  that  an  APG  causes  the 
turbulent  shear  stress  profile  in  a  plot  of  -pw'v'/r^  to  overshoot  the  wall  shear  stress  value 
before  curling  back  towards  unity. 

It  can  be  shown  that  the  boundary  layer  shear  stress  should  tend  toward  the  value  of  the 
wall  shear  stress  by  integrating  the  boundary  layer  momentum  equation  from  the  wall  to  a  point  in 
the  boundary  layer  a  distance  y  from  the  wall  (where  y  is  still  within  the  boundary  layer).  The 
integration  leads  to  the  relation  [6] 


du  T 
M—  +  r 
dy 


+~y 

ax 


(7.2) 


For  a  flat  plate,  dp^ldx^Q,  so 


du  Y 

dy 


(7.3) 


In  the  logarithmic  region  p{duldy)  «  tly,  so  Tly=T„.  Recall  from  Chapter  2  that 


=  -pu'v'  -  u p'v’  -vp’u’  -  p'u'v' .  Neglecting  the  last  three  terms,  =  -pu'v' ,  and  the 
profile  in  Figure  7.13  was  expected  to  tend  towards  unity.  Note  that  if  the  last  three  terms  in 


were  not  negligible,  the  wall  shear  stress  would  be  underestimated. 

The  fact  that  the  profile  in  Figure  7.13  was  expected  to  tend  towards  unity  was  used  to 
estimate  a  value  for  the  “actual”  wall  shear  stress  by  extrapolating  the  curve.  The  value  of  wall 
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shear  stress  found  using  this  process  was  70  Pa,  a  4.3%  difference  from  the  value  calculated  using 
van  Driest’s  relationship.  The  measured  values  for  the  incompressible  Reynolds  shear  stress  are 
presented  in  Figure  7.14  scaled  by  this  new  value  of  the  wall  shear  stress.  Note  that,  when  scaled 
in  this  manner,  the  measured  data  fell  almost  exactly  on  top  of  KlebanofTs  curve  for 
incompressible  flow,  lending  support  to  Morkovin’s  hypothesis. 

The  second  method  of  non-dimensionalization  was  to  scale  the  correlation  -u'V  by  the 
square  of  the  local  velocity,  .  This  scaling  was  used  for  three  reasons.  First,  this  scaling 
eliminated  any  linear  biasing  which  may  have  influenced  the  velocity  measurements  (see  Appendix 
B).  Second,  this  scaling  enabled  easier  comparison  with  CFD  results.  Note  that  -u'v'ju^  is 
equivalent  to  the  incompressible  Reynolds  shear  stress  divided  by  ;  CFD  turbulence  models 


-pu'v'  /Tw 

Figure  7. 14:  Incompressible  Reynolds  Shear  Scaled  by  Extrapolated  Wall  Shear  Stress 
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cannot  resolve  the  term  u'v' ,  only  «  -pu'v' .  Therefore,  use  of  this  scaling  allowed  for  simple 

comparison  to  CFD  results.  Finally,  and  most  importantly,  scaling  by  allowed  the  presentation 
of  a  meaningful  turbulent  shear  stress  profile  that  only  contained  information  from  direct 
measurements.  Thus,  scaling  the  incompressible  Reynolds  shear  stress  by  the  local  velocity 
eliminated  the  possibility  of  incurring  errors  associated  with  inaccurate  calculation  of  other  scaling 
parameters.  LDV  measurements  of  the  incompressible  Reynolds  shear  stress  scaled  using  local 
quantities  are  presented  in  Figure  7.15.  Note  the  excellent  agreement  with  the  hot-wire  data  from 
Ref  [26].  This  agreement  showed  that  the  assumption  of  /?'  =  0  used  to  reduce  the  hot-wire  data 
was  valid  for  flat  plate  flow. 


Figure  7.15:  Incompressible  Re5^olds  Shear  Stress  Scaled  by  Local  Quantities 
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7.3.6  Van  Driest  Velocity  Profile.  The  mean  velocity  was  plotted  using  the  scaling  laws 
given  by  van  Driest  [35]  as  presented  in  Section  6.8.  This  velocity  profile  is  plotted  in  Figure 
7.16,  along  with  the  flat  plate  correlation  given  in  Eq.  (6.47).  For  a  flat  plate,  Clauser’s 
equilibrium  parameter  was  equal  to  0,  leading  to  a  value  of  IT  =  0.476  for  the  Coles  wake 
parameter. 

The  solid  line  in  Figure  7.16  represents  a  semi-empirical  curve  which  has  been  proven 
throughout  the  literature  to  be  accurate,  so  the  slight  disagreement  of  the  present  data  to  the  curve 
seemed  a  bit  disheartening.  However,  it  should  be  pointed  out  that  the  data  points  in  Figure  7.16 
were  generated  using  the  estimated  wall  shear  stress  of  70  Pa;  as  discussed  in  Section  7.3.5,  this 
value  was  probably  underestimated.  Increasing  the  wall  shear  stress  tended  to  shifl:  the  points 
closer  to  the  semi-empirical  curve  (a  value  of  =15  Pa  collapsed  the  data  the  best),  pointing  out 


Figure  7.16:  Flat  Plate  Van  Driest  Velocity  Profile 
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once  again  that  the  wall  shear  stress  was  probably  underestimated.  It  should  be  mentioned  here 
that  a  value  of  =  75  Pa  is  just  inside  the  documented  accuracy  of  van  Driest’s  skin  friction 
relation  (±10% ,  cf.  Ref  [9])  when  compared  to  the  estimated  value  of  =  70  Pa.  However,  to 
accurately  preent  the  measured  incompressible  Reynolds  shear  stresses,  it  is  highly  recommended 
that  direct  measurements  of  the  wall  shear  stress  be  obtained  experimentally. 

7. 3. 7  Turbulence  Statistics 

7.3. 7.1  Skewness,  Flatness,  and  Intermittency.  The  values  of  skewness,  flatness, 
and  intermittency  were  important  when  comparing  the  LDV  measurements  to  data  in  the  literature 
since  they  helped  determine  if  the  present  measurements  were  being  biased  by  the  seeding  process. 
For  example,  non-uniform  seeding  in  the  current  experiment  would  cause  the  measured  velocity 
profiles  to  be  skewed  differently  than  they  were  in  other  experiments. 

The  probability  of  biasing  is  also  increased  in  compressible  flows  due  to  temperature  and 
density  gradients  in  the  boundary  layer.  A  warmer,  less-dense  “clump”  of  fluid  may  be  able  to 
hold  more  seed  particles  than  a  clump  of  cold,  dense  fluid  (it  would  have  more  “free  space” 
between  the  fluid  molecules  for  the  seed  particles  to  occupy).  The  LDV  system  would  then  register 
more  bursts  for  the  warmer  clump  of  fluid  than  the  cold  clump.  This  would  skew  the 
measurements  towards  the  velocity  of  the  warmer,  less-dense  clump  of  fluid  when,  in  reality,  they 
should  be  skewed  towards  the  velocity  of  the  cold,  dense  clump  of  fluid. 

In  order  to  determine  if  the  measurements  were  being  biased,  plots  of  the  skewness  and 
flatness  were  compared  to  data  in  the  literature.  Unfortunately,  most  of  the  data  available  for  the 
skewness  and  flatness  in  compressible  flows  was  data  collected  in  other  LDV  experiments.  Thus, 
if  the  skewness  and  flatness  data  in  the  literature  were  subject  to  biasing,  comparing  the  present 
data  to  that  which  had  the  same  possibility  of  being  biased  would  be  problematic.  Obviously,  there 
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exists  a  need  for  skewness  and  flatness  data  collected  by  means  which  are  not  subject  to  the 
possible  density  biasing  discussed  here.  However,  for  lack  of  other  data  with  which  to  compare, 
the  LDV  skewness  and  flatness  data  in  the  literature  has  been  taken  as  valid.  The  skewness  data 
collected  in  the  present  study  is  shown  in  Figure  7.17  compared  to  data  from  Ref  [16].  The  data 
from  the  present  study  agrees  well  with  that  from  Ref  [16],  indicating  that  the  methods  used  in  this 
study  were  as  valid  as  those  used  by  past  researchers. 

The  flatness  data  is  presented  in  Figure  7.18  and  compared  to  data  fi'om  Ref  [16].  The 
correlation  was  once  again  very  good,  especially  considering  the  scatter  in  the  data.  The  flatness 
did  not  provide  much  insight  into  the  physics  of  the  flowfield,  but  examination  of  the  intermittency 
fiinction  (which  is  calculated  using  the  flatness)  did  provide  useful  information. 


Figure  7. 17:  Skewness  of  Velocity  Components  in  Flat  Plate  Flow 
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Figure  7.18:  Flatness  of  Velocity  Components  in  Flat  Plate  Flow 

The  profile  of  the  intermittency  function  is  shown  in  Figure  7.19.  Even  though  the  flatness 
curve  in  Figure  7. 18  agreed  well  with  the  data  from  Ref  [16],  there  was  a  noticeable  difference  in 
the  intermittency  profile.  The  curve  fit  to  the  data  from  Ref  [16]  had  ^  =  0.1 14S  and  ^  =  0.64 , 
whereas  in  the  present  experiment  they  were  ^  =  022S  and  (p  =  0.68 .  This  showed  that  the 
boundary-layer  edge  was  essentially  at  the  same  location  as  in  Ref  [16],  but  the  standard  deviation 
from  the  mean  position  was  almost  twice  as  large.  The  increase  in  the  standard  deviation  was 
probably  due  to  the  1%  freestream  turbulence  level  measured  in  the  current  experiment,  compared 
with  a  freestream  turbulence  level  of  about  0.5%  measured  in  Ref  [16]. 

7.i.  7.2  Correlation  Coefficient.  The  profile  of  the  correlation  coefficient  is 
presented  in  Figure  7.20.  Near  the  wall  -R^  «  0.42  ;  recall  from  Chapter  6  that  experimental 
data  indicated  -R^  =  0.4  -  0.5  for  low-speed,  incompressible  flow.  [6]  Thus,  the  compressibility 
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Figure  7.19:  Flat  Plate  Intermittency  Function  Profile 
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Figure  7.20:  Flat  Plate  Correlation  Coefficient 


7-22 


effects  associated  with  high-speed  flow  did  not  appear  to  have  a  large  influence  on  the  correlation 
coefficient.  In  addition,  the  present  results  showed  a  reduction  in  the  correlation  coefficient 
throughout  the  boundary  layer  when  compared  to  the  data  of  Ref.  [16],  possibly  due  to  the  higher 
freestream  turbulence  level  in  the  present  facilities. 

7. 3. 7. 3  Turbulent  Kinetic  Energy.  Results  from  the  TKE  calculations  are 
presented  in  Figure  7.21  along  with  CFD  results  from  Ref  [17].  The  CFD  calculations  showed  a 
nearly-linear  profile  through  the  boundary  layer,  while  the  measurements  tended  to  decrease  near 
the  wall.  This  roll-off  was  due  to  the  estimation  of  vr' «  v'  (see  Section  6.9):  in  the  literature,  w' 
has  been  shown  to  be  about  1.5  v'  near  the  wall,  while  it  is  approximately  equal  to  v'  in  the 
freestream.  [6]  This  under-estimation  of  w'  near  the  wall  led  to  the  inaccuracy  in  the  computed 
TKE  in  this  region.  Also  note  that  the  CFD  results  consistently  underestimated  the  measured  TKE. 
Once  again,  this  was  probably  due  to  the  non-zero  freestream  turbulence  level  in  the  wind  tunnel. 


Figure  7.21:  Flat  Plate  TKE  Profile 
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7.4  FPG  Results 


With  the  testing  procedures  validated  against  existing  data  and  producing  viable  results, 
the  data  from  the  FPG  testing  were  compared  to  the  flat  plate  data.  The  FPG  results  and  their 
implications  are  discussed  below. 

7. 4. 1  Mean  Flow  Characteristics 

7. 4. 1.1  Velocity.  In  the  FPG  test  section,  was  found  to  be  approximately  607 
m/s.  Scaling  the  measured  velocities  by  this  value,  the  non-dimensional  velocity  profile  was  used 
to  calculate  a  value  for  Following  the  same  procedure  as  that  used  for  the  flat  plate, 
=  1 1.9  mm;  the  resulting  non-dimensional  profile  is  shown  in  Figure  7.22.  The  FPG  did  not 


Figure  7.22:  FPG  Velocity  Profile  in  Non-Dimensional  Coordinates 
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have  a  large  effect  on  the  non-dimensional  velocity  profile;  in  the  outer  region  of  the  boundary 
layer  the  velocity  defect  was  increased  while  in  the  near-wall  region  it  decreased,  but  in  either  case 
the  change  was  very  slight. 

7. 4. 1.2  Mach  Number  and  Density.  The  Mach  number  profile  is  presented  in 
Figure  7.23,  white  the  density  profile  is  presented  in  Figure  7.24.  The  Mach  number  profile  shows 
good  agreement  with  the  measurements  presented  in  Ref  [26],  especially  when  considering  the 
large  errors  associated  with  the  value  of  5  (see  Appendix  B).  Recall  that  the  density  profile  in  the 
FPG  region  was  calculated  in  a  different  manner  than  for  the  flat  plate  (see  Section  6.4.2);  in  order 
to  determine  the  validity  of  the  calculations,  the  density  profile  was  compared  to  CFD  results  from 
Ref  [17].  Once  again,  the  profile  matched  extremely  well  and  the  assumptions  used  to  calculate 
the  density  (the  assumptions  of  constant  total  temperature  and  a  linear  pressure  distribution 
through  the  boundary  layer)  were  justified. 


Figure  7.23:  FPG  Mach  Number  Profile 
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Figure  7.24:  FPG  Non-Dimensional  Density  Profile 

7. 4. 1.3  Integral  Momentum  Thicknesses.  The  integral  momentum  thickness 
results  for  the  FPG  region  are  presented  in  Table  7.3.  Note  that,  when  comparing  these  results  to 
the  flat  plate  values  given  in  Table  7.2,  the  thicknesses  tended  to  increase  dramatically.  This  at 
first  seemed  counter-intuitive,  as  a  stabilizing  FPG  should  be  expected  to  decrease  the  thicknesses. 
It  appeared  that  this  may  have  been  the  case  between  x  =  70.44  and  x  =  70.98  cm,  where  the 
thicknesses  decreased  slightly.  In  this  region,  the  wall  was  only  slightly  curved,  and  the  stabilizing 
effect  of  the  FPG  was  strong  enough  to  overcome  the  effect  of  the  wall  curvature.  However,  at  the 
location  x  =  70.5  cm  the  wall  curvature  outweighed  the  effects  of  the  FPG,  leading  to  an  increase 
in  the  integral  momentum  thicknesses. 
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Table  7.3:  FPG  Integral  Thicknesses  Calculated  From  LDV  Measurements 


X  (cm) 

5*  (mm) 

e  (mm) 

S*  (mm) 

6',  (mm) 

70.44 

2.80 

1.11 

0.41 

0.67 

70.98 

2.76 

1.09 

0.41 

0.66 

71.52 

3.95 

1.43 

0.56 

0.89 

7.4.2  Turbulence  Intensities.  The  u-  and  v-component  turbulence  intensities  are 
presented  in  Figure  7.25  and  compared  to  the  flat  plate  values.  Note  that,  in  both  cases,  the  FPG 
drastically  reduced  the  turbulence  intensities  in  the  lower  half  of  the  boundary  layer  but  had 
relatively  no  effect  above 0.5. 

7.4.3  Fluctuating  Mach  Number.  The  fluctuating  Mach  number  results  for  the  FPG  test 
section  are  presented  in  Figure  7.26.  Near  the  wall,  the  FPG  tended  to  decrease  the  value  of  the 
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fluctuating  Mach  number  when  compared  to  the  flat  plate  results.  This  not  only  provided  evidence 
of  the  stabilizing  effect  of  the  FPG,  but  it  also  indicated  that  there  were  no  local  shocklets  present 
in  the  FPG  boundary  layer. 

7.4.4  Wall  Shear  Stress.  The  wall  shear  stress  in  the  FPG  region  was  calculated  using 
the  relationship  given  by  Eq.  (7.2).  From  Eq.  (7.1),  dp^Jdx  =  -2  x  10'*  Pa/m,  so  that 

+  2  X  y  (7.4) 

Using  the  value  of  -pu'v'  calculated  at  the  iimermost  point,  the  wall  shear  stress  in  the  FPG 
region  was  found  to  be  r„,  =  44  Pa  (C^  =  1.15  x  10"^).  When  compared  to  the  flat  plate  value  of 
approximately  70  Pa,  the  FPG  was  seen  to  reduce  the  wall  shear  by  approximately  40%. 
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7.4.5  “Incompressible”  Reynolds  Stresses.  The  incompressible  Reynolds  stresses 
measured  in  the  FPG  region  are  compared  to  the  flat  plate  results  in  Figure  7.27.  Note  that  the 
FPG  caused  a  significant  reduction  in  the  ratio  of  the  turbulent  shear  stress  to  the  wall  shear  stress 
when  compared  to  the  flat  plate  curve.  Noting  that  the  magnitude  of  the  wall  shear  stress  was 
reduced  by  the  FPG,  the  reduction  in  the  absolute  magnitude  of  -~pu'v'  was  higher  than  that 
implied  by  the  curves  in  Figure  7.27. 

A  better  feel  for  the  effect  of  the  FPG  on  the  absolute  magnitude  of  -Jm'v’  was  achieved 
by  scaling  the  incompressible  Reynolds  stress  by  the  local  density  and  velocity;  this  scaling  was 
used  in  the  data  presented  in  Figure  7.28.  When  scaled  by  local  parameters,  the  FPG  was  seen  to 
reduce  -'pii'v'  by  approximately  75%  near  the  wall  (as  opposed  to  the  60%  reduction  seen  when 


-p  uV  /-Cw 

Figure  7.27:  FPG  Incompressible  Reynolds  Stresses  Scaled  by  Wall  Shear  Stress 
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-uV  /  u^ 

Figure  7.28:  FPG  vs.  ZPG  Incompressible  Reynolds  Shear  Stress  Scaled  by  Local  Quantities 


scaled  by  ).  This  was  the  expected  result  based  on  the  reduced  wall  shear  stress  and  heat  flux 
values  reported  by  Spina  et  al.  for  FPG  flows.  [33]  Note  that,  above  y/S  »  0.5 ,  the  velocity 

fluctuations  were  found  to  be  practically  uncorrelated  (i.e.,  u'v'  w  0).  The  stabilizing  effect  of  the 
FPG  was  more  than  likely  responsible  for  this  effect. 

The  hot-wire  data  from  Ref  [26]  was  also  presented  in  Figure  7.28  for  comparison.  The 
hot-wire  data  followed  the  same  trends  as  the  LDV  data,  but  the  magnitudes  of  the  hot-wire  data 
were  substantially  larger.  Recall  that  the  agreement  between  the  flat  plate  hot-wire  and  LDV  data 
(see  Figure  7.15)  was  excellent.  The  large  difference  in  the  FPG  region  was  more  than  likely  the 
result  of  inaccurate  assumptions  in  the  reduction  of  the  hot-wire  data.  That  is,  the  hot  wire  data 
reduction  assumed  /?'  =  0 ;  the  results  in  Figure  7.28  indicated  that  this  assumption  was  probably 
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invalid  in  the  FPG  region.  In  order  for  the  assumption  of  p'  =  0  to  become  invalid  in  the  FPG 
region,  the  FPG  would  have  to  cause  the  pressure  fluctuations  to  increase.  In  contrast,  to  this 
point  the  FPG  has  been  seen  to  stabilize  and  reduce  all  other  fluctuating  quantities.  If  the  size  of 
the  pressure  fluctuations  was  magnified  by  the  FPG,  the  increase  was  probably  due  to  the 
streamline  curvature  and  the  associated  pressure  change  across  the  boundary  layer  rather  than  the 
streamwise  pressure  gradient.  That  is,  the  non-zero  value  of  dpjdy  across  the  FPG  boundary 
layer  may  have  provided  a  mechanism  for  generating  p' .  Using  a  Prandtl  mixing-length  model, 
p' « lp{dp/dy) ,  where  is  a  mixing  length  associated  with  the  pressure  fluctuations.  Thus,  over 

the  flat  plate  dpjdy  »  0 ,  leading  to  /?' «  0 ;  in  the  FPG  region  dpjdy  0 ,  so  p'  ^0. 

Finally,  Figures  (7.27)  and  (7.28)  showed  that  the  profile  of  the  FPG  incompressible 
Reynolds  shear  stress  actually  became  negative  around  yjS»  0.65 .  This  negative  value  could 
prove  problematic  when  trying  to  develop  a  turbulence  model  for  use  in  CFD  codes.  Due  to  the 
fact  that  CFD  turbulence  modeling  has  focused  on  flat  plate  and  adverse  pressure  gradient  flows, 
current  CFD  models  (such  as  the  model  used  in  Ref.  [17])  do  not  admit  the  possibility  of  <  0 . 

7.4.6  Van  Driest  Velocity  Profile.  The  van  Driest  velocity  profile  for  the  FPG  region  is 
presented  in  Figure  7.29.  The  FPG  was  found  to  slightly  increase  the  effective  velocity  when 
compared  to  the  calculated  flat  plate  value;  however,  recall  that  the  flat  plate  value  was  fairly 
inaccurate  due  to  the  error  in  the  computed  wall  shear  stress.  When  comparing  the  FPG  data  to  the 
flat  plate  empirical  curve  fit,  the  FPG  was  seen  to  have  a  substantial  effect  on  the  effective  velocity 
profile.  It  should  be  noted  that  in  the  FPG  region,  p--\3,  leading  to  an  undefined  value  for  the 
Coles  wake  parameter,  11 .  Therefore,  it  was  not  possible  to  calculate  a  correlation  to  use  for 
comparison  with  the  FPG  profile. 
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Figure  7.29:  FPG  van  Driest  Velocity  Profile 

7. 4. 7  Turbulence  Statistics 

7  4. 7 1  Skewness,  Flatness,  and  Intermittency.  Comparison  of  the  flat  plate  and 
FPG  skevmess  (see  Figure  7.30)  showed  that  the  FPG  reduced  the  magnitude  of  both  the  u  and  v 
components  of  skewness.  In  the  flat  plate  case,  the  high  negative  value  of  and  high  positive 
value  of  Sk^  showed  that  more  particles  were  crossing  through  the  control  volume  from  the  side 
nearest  the  wall  (“below”  the  control  volume).  The  FPG  reduced  the  skewness,  indicating  a  more 
even  balance  in  the  number  of  particles  entering  the  control  volume  from  above  and  below. 

As  in  the  flat  plate  case,  examination  of  the  flatness  in  the  FPG  region  did  not  provide  a 
great  deal  of  insight  into  the  effects  of  the  FPG  on  the  flow,  but  the  measurements  are  compared  to 
the  flat  plate  values  in  Figures  7.31  and  7.32  for  completeness.  In  contrast  to  the  flatness,  the 
intermittency  fimction  provided  some  insight  into  ±e  stabilizing  effects  of  the  FPG.  The 
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intermittency  profile  for  the  FPG  region  is  compared  to  the  flat  plate  intermittency  in  Figure  7.33. 
The  error  function  curve  fit  showed  that,  for  the  FPG  region,  (p  =  0.8  and  ^  =  0.17^ .  Recall  that 
the  flat  plate  curve  fit  resulted  in  ^  =  0.68  and  ^  =  0.22^ .  Since  (p  represented  the  random 
location  of  the  boundary  layer  edge,  the  FPG  was  seen  to  push  this  location  outward  towards  d . 
The  FPG  also  slightly  reduced  the  standard  deviation  around  (p ,  but  not  by  a  significant  amount. 
Recall  from  the  turbulence  intensity  results  in  Figure  7.25  that  the  turbulence  was  not  affected 
much  above  yfd  »  0.5 ;  since  the  random  position  of  the  boundary  layer  edge  was  outside  this 
value,  the  standard  deviation  represented  by  should  not  have  changed  significantly.  Comparison 
of  the  magnitudes  of  the  flat  plate  and  FPG  intermittency  profiles  also  showed  that,  although  the 
FPG  decreased  the  magnitudes  of  the  velocity  fluctuations,  it  caused  the  flow  to  be  turbulent  for  a 
higher  percentage  of  time:  i.e.,  the  FPG  was  distributing  the  energy  from  a  few  particles  with  large 
velocity  fluctuations  over  many  particles,  each  with  a  smaller  fluctuating  velocity. 
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Figure  7.30:  FPG  vs.  ZPG  Skewness  Comparison 
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u  -  Flatness 

Figure  7.31:  «-Coniponent  Flatness  Comparison 
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Figure  7.32:  v-Component  Flatness  Comparison 
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7.4. 7.2  Correlation  Coefficient.  The  calculated  correlation  coefficient  profile  for 
the  FPG  region  is  compared  to  the  flat  plate  profile  in  Figure  7.34.  Recall  that  above  yjd  «  0.5 

the  FPG  caused  a  reduction  in  the  velocity  correlation  u'v'  while  leaving  the  turbulence  intensities 
relatively  unaffected.  These  effects  were  reflected  in  the  fact  that  the  correlation  coefficient  was 
approximately  zero  above  y/(5«0.5.  Below  the  flow  velocity  decreased  due  to  the 

presence  of  the  wall,  causing  the  correlation  coefficient  to  tend  toward  the  low-speed  value 
associated  with  flat  plate  flow  ( -  0.4  -  0.5  ). 

7. 4.7. 3  Turbulent  Kinetic  Energy.  The  FPG  TKE  profile  is  shown  in  Figure 
7.35;  however,  due  to  the  estimation  of  w' =  v' ,  the  accuracy  of  this  profile  was  very 
questionable,  especially  near  the  wall.  Note  the  excellent  comparison  to  the  CFD  results;  the 
discrepancy  near  the  wall  was  probably  due  to  an  underestimation  of  w'  in  this  region. 
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Figure  7.34:  FPG  and  ZPG  Correlation  Coefficient  Comparison 
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Figure  7.35:  FPG  TKE  Profile 
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7.5  Strain  Rate  Profiles 

Strain  rate  profiles  were  calculated  across  the  boundary  layer  for  both  the  flat  plate  and  the 
FPG;  the  flat  plate  strain  rates  calculated  were  dii jdx ,  dujdy ,  dvjdx  and  dvjdy .  The  strain 
rates  calculated  in  the  FPG  region  also  included  dujdz  and  dvjdz .  The  results  from  these 
calculations  are  presented  in  Figures  7.36  to  7.41.  It  should  be  noted  that  the  strain  rates  were 
calculated  using  finite-difference  methods,  and  as  a  result  the  error  in  the  calculations  was 
relatively  large  (see  Appendix  B). 

Upon  examination  of  Figures  7.36  to  7.39,  it  was  apparent  that  the  FPG  did  not  have  a 
large  effect  on  the  strain  rates  in  the  direction  normal  to  the  wall.  The  comparison  of  the  flat  plate 
and  FPG  values  of  Sii  jdx  showed  that  the  FPG  did  not  affect  the  strains  except  in  the  very  near¬ 
wall  region,  where  they  went  negative.  The  FPG  was  seen  to  shift  dvjdx  from  the  flat-plate 
value  of  zero  to  a  finite  negative  value;  this  value  was  essentially  constant  across  the  boundary 
layer.  The  values  of  du  jdz  and  dv  jdz  were  also  measured  for  the  FPG  test  section;  the  values 
of  these  strain  rates  were  found  to  be  at  least  an  order  of  magnitude  below  any  other  measured 
strain  rate  and  were  considered  negligible. 

Finally,  Smith  and  Smits  [32]  pointed  out  that  Clauser’s  equilibrium  parameter  is  not 
always  sufficient  for  characterizing  the  strength  of  a  pressure  gradient  in  supersonic  flows.  For 
instance,  a  pressure  gradient  generated  by  a  curved  wall  will  affect  the  strain  rates  differently  than 
the  same  pressure  gradient  generated  over  a  flat  plate.  As  a  result,  they  suggested  classifying  a 
distortion  as  weak  if  the  ratio  of  each  extra  strain  rate,  e,  to  the  main  strain  rate,  {dujdy) ,  was 

approximately  0.01.  A  distortion  where  0.1  was  considered  strong.  The  maximum 

value  of  el{dujdy)  present  in  this  study  was  for  the  dvjdy  strain  rate;  its  magnitude  was 
approximately  0.1  and  the  distortion  created  by  the  FPG  was  classified  as  strong. 
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Figure  7.36:  du  jdx  Strain  Rates  Over  Flat  Plate  and  in  FPG  Test  Section 
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Figure  7.37:  dujdy  Strain  Rates  Over  Flat  Plate  and  in  FPG  Test  Section 
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Figure  7.39;  ffv  jdy  Strain  Rates  Over  Flat  Plate  and  in  FPG  Test  Section 


7. 6  Progression  of  Turbulent  Quantities  in  FPG  Test  Section 

The  data  collected  for  use  in  generating  the  strain  rates  allowed  a  comparison  of  the 
turbulent  quantities  as  they  developed  through  the  FPG  test  section.  The  progression  of  the  u 
turbulence  intensity  and  the  incompressible  Reynolds  shear  stresses  are  shown  in  Figures  7.42  and 
7.43,  respectively.  In  both  cases  the  intermediate  data  fell  in  between  the  flat  plate  and  FPG  test 
section  data,  as  was  expected. 


Figure  7.42:  Progression  of  u  Turbulence  Intensity  Through  Boundary  Layer 
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Figure  7.43:  Progression  of  Incompressible  Reynolds  Shear  Stress  Through  FPG  Test  Section 
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8.  Summary,  Conclusions,  and  Recommendations 


8. 1  Summary  and  Conclusions 

The  main  goal  of  the  present  study  was  to  increase  the  database  of  available  turbulence 
information  for  compressible  flows.  Specifically,  the  effects  of  a  favorable  pressure  gradient 
generated  by  a  convex  curved  wall  were  studied  in  detail.  The  data  in  the  FPG  region  showed  that 
the  FPG  had  a  stabilizing  effect  on  the  turbulent  quantities.  The  magnitudes  of  the  turbulence 
intensities  were  reduced  (when  compared  to  the  corresponding  flat  plate  data)  below  the  value 
yjS « 0.5 ;  above  this  value,  the  turbulence  intensities  were  unaffected  by  the  FPG.  The 
magnitude  of  the  incompressible  Reynolds  shear  stress  was  also  reduced  by  the  FPG;  in  the  near¬ 
wall  region  the  magnitude  was  approximately  25-30%  of  the  flat  plate  value. 

The  FPG  was  also  seen  to  affect  the  accuracy  with  which  the  hot-wire  data  could  be 
reduced.  That  is,  discrepancies  between  the  hot-wire  data  and  the  LDV  data  in  the  FPG  region 
indicated  that  the  assumption  of  p’  =  Q  used  in  the  reduction  of  the  hot-wire  measurements  was  not 
valid  in  the  FPG  test  section. 

In  addition  to  studying  the  effects  of  the  FPG  on  the  turbulence,  the  current  experiment 
used  methods  which  took  steps  to  improve  on  the  LDV  measurements  of  flat  plate  data  presented  in 
the  literature.  Specifically,  the  fringe  disks  in  the  laser  Doppler  measurement  volume  were  aligned 
perpendicular  to  the  mean  flow  direction  in  order  to  prevent  the  problems  associated  with  angular 
bias.  Comparison  and  analysis  of  raw  and  coincident  data  (see  Appendix  A)  showed  that  this 
“angular”  bias  may  also  have  been  due  to  the  choice  of  record  interval.  In  any  event,  the 
procedures  used  in  the  current  experiment  eliminated  the  near-wall  roll-off  of  the  turbulent  shear 
stresses  observed  by  other  researchers. 
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Finally,  in  order  for  this  study  to  be  useful  as  a  source  of  turbulence  information,  there 
were  a  number  of  variables  which  needed  to  be  reported  for  classifying  the  flow  and  allowing 
simple  manipulation  of  the  results.  The  variables  which  were  determined  as  being  necessary  to 
accomplish  this  purpose  are  presented  in  Tables  8.1  and  8.2. 


Table  8.1:  Flow  Field  Parameters 


Variable 

Flat  Plate 

FPG 

(kg/m^) 

0.102 

0.0626 

dpy,ldx  (Pa/m) 

0 

-2  X  10^ 

C/ 

1.63  X  10'^ 

1.15  X  10“^ 

(Pa) 

70.0 

44.0 

Re^ 

8880 

9160 

1.75  X  10^ 

1.97  X  10^ 

Re* 

1.26  X  10^ 

1.17x10^ 

Variable 

Flat  Plate 

FPG 

To,.  (K) 

298 

298 

Po.^  (Pa) 

2.13  X  10^ 

2.13x10^ 

7;(K) 

295 

295 

u,  (m/s) 

602 

607 

2.79 

2.91 

7;(K) 

116 

no 

Pe  (Pa) 

7870 

6530 

Pe  (kg/m^) 

0.0237 

0.207 

Table  8.2:  Boundary-Layer  Structure  Parameters 


Variable 

Flat  Plate 

FPG 

5*  (mm) 

1.31 

3.95 

6  (mm) 

1.27 

1.43 

(mm) 

0.615 

0.887 

el{duldy)\^^ 

0 

0.1 

Variable 

Flat  Plate 

FPG 

Su  (mm) 

9.85 

11.9 

Sm  (mm) 

10.8 

13.1 

5*  (mm) 

0.499 

0.556 

P 

0 

-1.34 

0.17 

0.22 

8.2  Recommendations 

The  data  gathered  in  the  current  experiment  has  been  proven  useful  as  a  basis  against 
which  CFD  turbulence-models  can  be  evaluated  {cf.  Ref  [17]);  however,  there  are  several  aspects 
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of  the  experiment  which  could  be  improved  in  order  to  make  the  data  even  more  valuable  in  the 
future. 

The  first  recommendation  is  to  obtain  accurate  measurements  of  the  wall  shear  stress  over 
the  flat  plate  and  in  the  FPG  region.  The  wall  shear  stress  is  the  determining  factor  in  the  heat 
transfer  to  the  surface  and  the  magnitude  of  the  friction  drag.  The  current  LDV  setup  was  limited 
to  measurements  a  finite  difference  from  the  wall,  making  it  hard,  if  not  impossible,  to  determine  a 
value  for  the  wall  shear  stress.  For  instance,  the  most  accurate  estimation  of  seemed  to  come 

from  the  relationship  given  in  Eq.  (7.2);  however,  if  the  terms  -up'v' ,  -vp'u' ,  or  -p'u'v'  are 
not  negligible,  even  this  estimation  is  not  very  accurate.  For  instance,  recall  from  Section  7.3.6 
that  Eq.  (7.2)  gave  the  flat  plate  shear  stress  as  approximately  70  Pa,  while  a  value  of  80  Pa  was 
required  to  collapse  the  data  onto  the  empirical  van  Driest  curve.  Obtaining  accurate 
measurements  of  the  wall  shear  stress  would  not  only  aid  in  the  analysis  of  the  data  in  this  report, 
but  it  would  be  useful  for  future  research  conducted  in  the  present  facilities  as  well. 

The  second  recommendation  is  to  collect  8-  and  16-sample  data  with  a  more  powerful  laser 
(one  capable  of  producing  an  adequate  SNR  in  the  freestream)  in  order  to  determine  if  the 
coincidence-filtering  bias  was  really  due  to  the  angular  alignment  of  the  lasers  or  if  it  could  be 
attributed  to  the  size  of  the  record  interval  (see  Appendix  A).  In  order  to  accomplish  this 
comparison,  the  study  should  include  (but  would  not  be  limited  to)  measurements  of  turbulence 
quantities  using  the  settings  shown  in  Table  8.3. 

The  third  recommendation  is  to  collect  skewness  measurements  via  some  means  which  are 
not  subject  to  the  possible  density  biasing  that  may  occur  in  compressible  flow  LDV 
measurements.  Comparison  of  the  skewness  profiles  would  help  confirm  or  deny  the  validity  of  the 
LDV  results  presented  in  this  report  and  throughout  the  literature. 
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Table  8.3:  Proposed  Settings  for  Angular  Bias/Record  Length  Study 


Record  Length 
(BSA  1) 

Record  Length 
(BSA  2) 

Coincidence  Window 
Size 

0° 

8  Samples 

8  Samples 

Twice  Record  Interval 

0° 

16  Samples 

16  Samples 

Equal  to  Record  Interval 

0° 

8  Samples 

16  Samples 

Twice  Record  Interval 

0° 

16  Samples 

32  Samples 

Equal  to  Record  Interval 

45° 

8  Samples 

8  Samples 

Twice  Record  Interval 

45° 

16  Samples 

8  Samples 

Equal  to  Record  Interval 

Fourth,  it  is  recommended  that  hot-wire  and  LDV  measurements  be  taken  over  a  flat  plate 
with  an  imposed  FPG  pressure  gradient  with  roughly  the  same  strength  as  the  one  generated  in  this 
study.  Comparison  of  the  hot-wire  and  LDV  data  would  show  whether  or  not  the  assumption 
/?'  =  0  is  valid  for  the  flat  plate,  imposed  pressure  gradient  case.  If  the  assumption  were  found  to 
hold,  it  would  be  obvious  that  the  discrepancies  between  the  hot-wire  and  LDV  data  collected  in 
the  current  experiment  were  the  result  of  the  streamline  curvature  and  not  the  FPG.  That  is,  it 
would  show  that  the  pressure  difference  across  the  boundary  layer  in  the  FPG  test  section  caused  in 
increase  in  the  magnitude  of  p' ,  even  though  the  FPG  tended  to  stabilize  all  other  fluctuating 
quantities. 

Finally,  although  there  is  no  known  method  for  doing  so  at  this  time,  measurements  of  the 
fluctuating  pressure  p’  should  be  achieved  in  order  to  determine  the  limits  on  the  validity  of  the 
hot-wire  data  reduction. 
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A.  Factors  Affecting  LDV Results 


While  collecting  LDV  data,  several  parameters  (record  length,  seeding  density,  coincidence 
window,  etc.)  were  found  to  affect  the  flow  and  were  varied  in  order  to  determine  their  effect  on  the 
measured  values  of  the  velocity.  In  addition,  past  researchers  have  noted  a  discrepancy  between 
the  coincidence-filtered  data  and  the  raw  data;  if  this  discrepancy  was  present  in  the  current 
experiment,  the  differences  would  have  to  be  noted  in  the  reported  results.  That  is,  the  presented 
data  in  Chapter  7  would  require  a  label  to  denote  whether  it  was  raw  or  coincidence-filtered. 
Therefore,  the  factors  which  would  affect  the  LDV  measurements  (and  thus  the  presented  results) 
are  examined  below. 

A.  1  Variation  of  Testing  Parameters 

As  mentioned  above,  several  testing  parameters  were  varied  initially  to  determine  their 
effects  on  the  measured  velocities.  The  three  most  important  parameters  were  the  record  length, 
seeding  density,  and  the  size  of  the  coincidence  window. 

A.  1.1  Record  Length.  Recall  that  the  record  length  was  confined  to  a  setting  of  either  8 
or  16  samples.  When  the  record  length  was  set  to  8  samples,  the  LDV  system  was  unable  to 
resolve  the  relatively  low  fireestream  turbulence  (see  Figure  A.l).  In  addition,  when  the  record 
length  was  set  to  8  samples  the  second-order  correlation  u'v'  showed  enormous  scatter  near  the 
wall  (Figure  A.2).  In  contrast,  the  16-sample  results  were  well-behaved  in  both  cases.  To  avoid 
the  inaccuracies  present  in  the  8-sample  measurements,  the  record  length  was  set  to  16  samples  for 
the  remainder  of  the  tests.  With  the  record  length  set  at  16  samples,  the  coincidence  window 
setting  was  restricted  to  0.333  ps,  the  length  of  the  record  interval. 
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Figure  A.  1 ;  8-SampIe  vs.  16-Sample  Turbulence  Intensity  Data  Over  a  Flat  Plate 
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Figure  A.2;  8-Sample  vs.  16-Sample  Second-Order  Correlation  Over  a  Flat  Plate 


A.  1.2  Seeding  Density.  In  accordance  with  the  adviee  of  Elena  [15],  the  seeding  density 
was  set  as  low  as  possible  to  minimize  the  effeets  of  the  seed  on  the  flow.  The  seeding  was 
controlled  via  a  pressure  regulator  on  the  atomizer;  the  optimum  (i.e.,  normal)  setting  was  found  to 
be  approximately  3  x  10^  Pa  with  two  atomizer  jets  opened.  Lower  seeding  pressures  ( ~  2.6  x  10^ 
Pa  with  two  jets  opened)  were  found  to  be  aeeeptable  in  the  flat  plate  case,  but  in  the  FPG  section 
the  seeding  density  decreased  as  the  streamlines  expanded.  This  deerease  in  the  seed  density 
(particularly  near  the  wall)  complicated  the  measurements  by  reducing  the  number  of  bursts 
collected  to  an  unaeeeptable  level.  In  order  to  eolleet  a  sufficient  number  of  bursts,  several  tests 
were  run  at  each  data  point  near  the  wall  and  the  turbulence  statistics  were  ealculated  by 
combining  the  burst  data  from  the  individual  tests.  To  avoid  these  complications,  a  seeding  density 
suitable  for  use  in  both  the  flat  plate  and  FPG  test  cases  was  chosen.  Figure  A.3  shows  the  effects 
of  the  normal  and  low  seeding  on  the  measurements  in  the  FPG  test  section.  The  change  in  the 
correlation  u'v'  between  the  two  eases  was  negligible.  Since  u'v'  was  a  second-order  term,  it 

should  have  been  the  most  sensitive  to  any  changes  in  the  flow  strueture;  u'v’  was  relatively 
unaffected  by  the  change  in  the  seeding  density,  indicating  that  use  of  the  higher  setting  did  not  bias 
the  results. 

A.  1.3  Coincidence  Window.  In  order  to  determine  the  effeet  the  size  of  the  coincidence 
window  had  on  the  coincidence-filtered  data,  several  of  the  data  points  in  the  interior  of  the 
boundary  layer  were  reproeessed  with  the  coineidence  window  set  to  0.033  ps  and  3.33  ps, 
providing  an  order-of-magnitude  difference  on  each  side  of  the  chosen  value  for  the  coincidenee 
window.  With  the  record  length  set  to  16,  there  was  absolutely  no  effeet  on  the  number  of 
coincident  partieles  when  the  window  was  varied  over  this  large  span.  Therefore,  the  chosen 
eoineidenee  window  was  proven  to  be  acceptable. 
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Figure  A.3:  Effects  of  Seeding  Density  on  Second-Order  Correlation  in  FPG  Test  Section 

A.2  Raw  vs.  Coincident  Data 

Throughout  the  literature,  researchers  have  reported  decreases  in  the  turbulence  intensities 
and  the  second-order  correlation  when  comparing  coincidence-filtered  data  to  the  raw 
measurements.  This  result  has  been  attributed,  at  least  in  part,  to  angular  biasing.  [15]  A 
decrease  in  the  turbulent  quantities  would  show  that  the  coincident  data  were  not  a  uniform  sample 
of  the  overall  data,  but  were  biased  toward  the  less-turbulent  particles.  Figures  A.4  and  A.5  show 
the  comparisons  between  the  raw  and  coincident  data.  The  agreement  was  excellent,  showing  that 
the  coincident  data  were  a  representative  sample  of  the  overall  data  collected.  The  biasing  which 
had  been  noted  by  previous  researchers  [16,19,30]  was  reduced  or  eliminated,  most  likely  by 
keeping  the  fringe  pattern  perpendicular  to  the  mean  flow,  thereby  eliminating  the  problem  of 
angular  bias. 
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Figure  A.4;  Raw  vs.  Coincident  Mean  Velocity  Profiles 
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Figure  A.5;  Raw  vs.  Coincident  Turbulence  Intensity  Profiles 


The  biasing  towards  slower  particles  which  has  been  observed  by  other  researchers  could 
also  have  been  attributed  to  the  choice  of  record  interval.  In  the  past,  the  record  interval  has  been 
set  to  approximately  half  the  traverse  time  of  the  fastest  particles  (Figure  A.6a).  In  the  present 
experiment,  with  the  record  length  set  to  16  samples,  the  record  interval  was  approximately  equal 
to  the  traverse  time  of  the  fastest  particles  (Figure  A.6b).  The  biasing  caused  by  coincidence 
filtering  has  been  observed  to  be  highest  near  the  wall,  where  the  particle  velocities  are  slower. 
Figures  A.6c  and  A.6d  present  the  burst  pattern  produced  by  particles  traveling  at  half  the  speed  of 
the  fastest  particles.  Figure  A.6c  shows  the  case  where  the  record  interval  was  set  to  half  the 
traverse  time  of  the  fastest  particles.  For  the  slow  particle,  the  record  interval  was  only  1/4  of  the 
traverse  time;  there  was  a  significant  amount  of  time  after  the  record  interval  where  the  burst 
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Figure  A.6;  Possible  Biasing  Due  to  Length  of  Record  Interval 
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pattern  was  still  increasing  and  the  BSA  could  start  tracking  the  same  particle  over  again.  The 
slower  the  particle,  the  higher  the  probability  that  the  same  particle  would  get  counted  twice.  In 
contrast.  Figure  A.6d  shows  the  case  where  the  slow  particle  was  tracked  with  the  record  interval 
set  equal  to  the  traverse  time  of  the  fastest  particles  (the  procedure  used  in  the  current  experiment). 
The  record  interval  ended  at  approximately  the  same  point  as  the  peak  of  the  burst  pattern;  since 
the  burst  intensity  was  decreasing  after  this  point,  the  BSA  would  not  begin  to  track  the  same 
particle  again. 

The  time  available  for  the  BSA  to  start  tracking  the  same  particle  over  again  was  not  only 
dependent  on  the  length  of  the  record  interval,  but  also  on  the  time  between  when  the  particle 
entered  the  control  volume  and  the  start  of  the  record  interval  (see  Figure  A.6c).  If  this  time  was 
large,  even  the  shortest  record  interval  would  extend  past  the  peak  of  the  Doppler  burst,  preventing 
the  BSA  from  tracking  the  same  particle  more  than  once.  The  time  spacing  between  when  the 
particle  entered  the  control  volume  and  the  start  of  the  record  interval  was  dependent  on  the  BSA 
burst  processing  algorithm  and  the  power  of  the  laser  beam. 

Recall  from  Chapter  3  that  the  BSAs  used  in  the  current  experiment  would  restart  the 
record  length  count  as  the  Doppler  burst  intensity  increased  past  certain  threshold  values,  ensuring 
the  peak  of  the  Doppler  burst  was  used  to  calculate  the  Doppler  fi-equency.  If  this  method  was  not 
incorporated  into  the  data  acquisition  process  (as  may  have  been  the  case  in  early  LDV 
equipment),  there  would  have  been  a  higher  probability  that  the  record  interval  did  not  extend  past 
the  peak  of  the  Doppler  burst  and  that  a  single  particle  could  have  been  tracked  more  than  once. 
Even  if  the  count-restarting  method  was  used,  its  advantages  could  have  been  negated  by  using  a 
laser  which  was  too  powerful:  as  the  laser  beam  power  increases,  the  intensity  of  the  scattered 
light  increases.  A  higher  intensity  of  scattered  light  would  enable  easier  detection  of  the  particles 
traveling  through  the  control  volume;  however,  if  the  laser  beam  was  too  powerful,  the  scattered 
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light  intensity  may  have  exceeded  the  final  BSA  threshold  at  the  first  or  second  sampling  time. 
This  would  have  led  to  a  short  time  between  when  the  particle  entered  the  control  volume  and  the 
start  of  the  record  interval,  increasing  the  possibility  of  tracking  the  same  particle  more  than  once. 
Advances  in  LDV  technology  allowed  the  current  measurements  to  be  taken  with  a  300  mW  laser, 
as  opposed  to  the  more  powerful  lasers  used  in  the  past  (the  data  from  Ref  [16]  was  collected 
using  a  9  W  laser),  helping  to  reduce  the  possibility  of  tracking  the  same  particle  more  than  once. 

Note  that,  if  the  record  interval  was  set  such  that  there  was  a  possibility  of  counting  the 
slower  particles  more  than  once,  the  raw  and  coincident  data  near  the  wall  would  both  be  biased 
towards  lower  velocities  than  were  actually  present  in  the  flow.  As  a  result,  one  might  believe  that 
the  coincident  data  would  still  follow  the  raw  data  closely,  but  this  would  not  be  the  case.  If  the 
same  particle  was  counted  twice  on  each  BSA,  the  coincidence  filtering  would  produce  not  two,  but 
four  coincident  bursts  (see  Figure  A.7),  further  biasing  the  coincident  data  toward  the  particles 
with  lower  velocities.  Note  that,  if  the  biasing  were  due  to  the  angular  alignment  alone  (as  is 
implied  in  the  literature),  both  the  raw  and  coincident  data  would  be  biased,  but  on  a  one-to  one 
ratio,  not  a  four-to-one  ratio.  That  is,  the  biasing  would  occur  because  the  BSA  fails  to  track  a 
particle  altogether,  not  because  the  effect  of  the  particle  is  altered  during  the  coincidence-filtering 
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Figure  A.7:  Number  of  Coincident  Bursts  Produced  by  a  Single  Particle 
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process.  If  there  were  a  problem  with  angular  biasing,  the  turbulent  quantities  would  exhibit  a  roll¬ 
off  near  the  wall,  but  the  roll-off  would  be  present  in  both  the  raw  and  coincident  data,  resulting  in 
a  negligible  difference  between  the  two  sets  of  data. 

In  order  to  determine  if  the  biasing  experienced  by  other  researchers  was  a  result  of  the 
record  interval  or  the  fringe  alignment,  the  8-sample  data  collected  in  the  present  experiment  was 
analyzed.  The  coincident  and  raw  turbulence  intensity  profiles  were  compared;  if  the  biasing  was 
present  in  the  coincident  data,  it  could  be  attributed  to  the  record  interval  setting  and  not  the  fringe 
aligmnent,  and  vice-versa.  Figure  A.  8  presents  the  comparison  of  8-sample  raw  and  coincident 
turbulence  intensities.  Near  the  wall,  there  was  a  slight  decrease  in  the  turbulence  intensity  due  to 
the  coincidence  filtering;  however,  contrary  to  what  was  expected,  the  decrease  became  larger  in 
the  freestream.  The  decrease  in  the  coincidence-filtered  data  near  the  outer  edge  of  the  boundary 
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Figure  A.8:  8-Sample  Raw  and  Coincident  Turbulence  Intensities 
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layer  was  probably  due  to  the  angular  alignment  of  the  optics  and  the  choice  of  record  interval  used 
for  the  v-component  measurements. 

Recall  from  Chapter  5  that  the  record  length  on  BSA  2  was  set  the  same  as  the  record 
length  on  BSA  1.  In  the  past,  when  the  lasers  were  aligned  ±45°  to  the  flow,  the  two  velocity 
components  measured  were  of  the  same  order  of  magnitude,  so  setting  the  two  record  lengths  equal 
to  one  another  was  required.  However,  with  the  lasers  oriented  so  the  fringes  were  perpendicular 
to  the  two  main  velocity  components,  the  lasers  were  measuring  velocities  of  different  magnitudes, 
so  the  record  length  should  have  been  set  accordingly  for  each  BSA.  Use  of  the  record  interval 
required  by  the  w-component  of  velocity  on  both  BSAs  meant  the  v-component  lasers  (which  were 
measuring  a  much  lower  velocity)  had  a  high  possibility  of  counting  the  same  particle  multiple 
times.  In  addition,  near  the  wall  the  turbulence  caused  fluctuating  velocities  with  magnitudes  up  to 
120  m/s.  In  contrast,  the  freestream  v-component  was  approximately  30-50  m/s  with  low 
fluctuations.  Thus,  on  the  average,  the  magnitude  of  the  measured  v-component  of  velocity  was 
lower  in  the  freestream.  Therefore,  the  record  interval  on  BSA  2  was  set  such  that  there  was  a 
high  possibility  of  tracking  the  same  particle  more  than  once,  and  the  probability  increased  towards 
the  freestream.  If  there  was  a  high  correlation  between  particles  with  lower  v  velocities  and  lower 
magnitudes  of  u'  (as  seems  reasonable),  then  the  biasing  present  in  Figure  A.8  may  have  been  due 
to  the  choice  of  record  interval  on  BSA  2. 

In  any  event,  the  biasing  near  the  wall  observed  in  Figure  A.8  was  larger  than  that 
observed  in  Figure  A.5,  indicating  that  the  majority  of  the  biasing  observed  in  this  and  past 
experiments  was  probably  caused  by  the  choice  of  record  interval  and  not  by  the  angular  alignment 
of  the  fringe  disks  in  the  LDV  measurement  volume.  In  order  to  determine  absolutely  if  this  was 
the  case,  it  is  recommended  that  a  detailed  study  into  the  possibility  of  record  length  biasing  be 
performed. 
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B.  Error  Analysis 


In  any  experiment,  there  is  a  certain  amount  of  error  associated  with  the  measurements.  It 
is  necessary  to  determine  the  magnitude  of  these  errors  in  order  to  determine  the  validity  of  the 
experiment.  If  the  actual  measured  quantities  are  used  to  compute  other  values,  the  error  of  the 
computed  values  needs  to  be  determined  as  well. 

In  order  to  determine  the  error  of  a  computed  quantity,  its  dependence  on  the  measured 
variables  needs  to  be  determined.  For  instance,  if  a  fimction  /  depends  on  x,  and  Xj ,  Xj  may 
have  a  greater  effect  on  /  than  X2 ;  as  a  result,  the  error  in  x,  would  affect  /  more  than  the 
error  in  Xj .  This  effect  is  quantified  by  differentiating  /  with  respect  to  each  of  the  dependent 
variables: 


A/  = 


<^x, 


Ax,  + 


l/o 


Axo 


(B.l) 


fo 


where  /q  is  the  mean  value  of  /  around  which  the  function  is  being  linearized.  It  is  often  more 
convenient  to  express  the  error  in  /  as  a  percent  of  the  measured  value  ;  in  this  case,  the 
percent  error  to  /  is  given  as  Sj-  ^  Af  ff^  x  100% : 


^/  = 


<?x, 


fo 


+ 

fo  ‘  0^-. 


(B.2) 


l/o 


fo 


'*2 


In  order  to  simplify  notation,  note  that  the  symbol  Af  will  be  used  to  denote  absolute  error  in  the 
quantity  / ,  while  will  be  used  to  denote  the  percent  error  in  / . 

Equations  (B.l)  and  (B.2)  give  the  maximum  possible  error  to  / ;  however,  if  the  errors 
in  X,  and  X2  are  independent,  it  may  be  reasonable  to  expect  that  they  will  not  both  be  at  their 
peak  values  at  the  same  time.  In  this  case,  the  error  to  /  may  be  given  by  a  Euclidean  norm: 
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df 


dx^ 
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f  ' 
Jo 


+ 


df 


r  X-i 

%  J 


(B.3) 


Unless  stated  otherwise,  error  calculations  performed  here  will  use  the  formulation  given  in  Eq. 
(B.3). 


In  the  following  analysis,  the  first  four  sections  examine  the  error  associated  with  the 
actual  measurements  (velocity,  stagnation  pressure,  and  stagnation  temperature);  once  these 
values  have  been  determined,  the  error  associated  with  calculated  quantities  is  evaluated. 


B.  1  Traverse  Location  Error 

Even  though  the  traverse  system  was  rated  as  being  accurate  to  ±0.13  pm/mm  (see 
Chapter  4),  the  location  error  was  greatly  increased  due  to  an  apparent  glitch  in  the  LDV  software. 
After  the  lasers  had  been  aligned  at  the  origin  location,  the  traverse  was  moved  to  the  first  data 
acquisition  location.  The  error  in  this  location  was  «0.15  pm.  However,  after  the  data  were 
collected  and  processed,  when  the  software  returned  to  the  traverse-control  window  the  location  of 
the  traverse  was  not  always  read  correctly.  For  example,  if  data  were  collected  at  the  location  y  = 
0.63  mm,  upon  returning  to  the  traverse-control  window  the  software  might  have  read  the  location 
as  y  =  0.62  mm.  If  the  error  occurred  (it  did  not  occur  every  time),  the  expected  location  was 
decreased  by  0.01  mm.  Therefore,  assuming  the  traverse  across  the  boundary  layer  started  at  the 
wall,  the  location  error  increased  as  the  distance  from  the  wall  increased.  The  data  were  collected 
along  a  line  normal  to  the  wall,  so  there  were  two  non-zero  traverse  coordinates  which  had  the 
possibility  of  being  affected  by  this  error;  however,  the  truncation  error  was  not  observed  to  affect 
more  than  one  of  the  two  coordinates  at  a  time.  Thus,  the  location  error  associated  with  the  data 
collection  point  away  from  the  wall,  Ayj ,  was 

=Ajo±  0.01(7-1)  (B.4) 
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where  Ajo  was  the  error  associated  with  placement  of  the  origin  and  units  were  in  mm.  The  error 


associated  with  the  origin  placement  and  the  error  in  the  probe  angle  alignment  are  discussed  in  the 
next  section. 


B.  2  Origin  Placement  and  Probe  Angle  Alignment  Error 

The  procedure  used  for  aligning  the  lasers  with  the  test  section  wall  was  described  in 
Chapter  5.  This  procedure  at  first  appeared  to  be  very  accurate,  but  as  this  section  will  point  out, 
the  angular  uncertainty  associated  with  the  measurements  turned  out  to  be  rather  large. 

The  dots  produced  by  the  lasers  on  the  wall  of  the  tunnel  were  measured  and  found  to  have 
a  diameter  of  0.60  mm.  The  distance  between  the  dots  was  measured  as  1.98  mm.  When  setting 
the  angle  of  the  lasers,  it  was  assumed  that  the  dots  were  placed  along  the  edge  of  the  wall  within 
1/4  radius  of  the  center.  Figure  B.l  shows  a  graphical  representation  of  the  maximum  error 
possible  under  these  assumptions. 


Figure  B.l:  Probe  Angle  Alignment  Error 


With  this  accuracy,  the  angular  error,  ,  was  given  by 

^0.60/8^ 


A^  =  tan 


=  4.289“  =0.0749  rads. 


(B.5) 


V  1.98/2, 

Assuming  the  dots  produced  by  the  lasers  could  be  aligned  within  1/4  radius  of  the  center 
also  gave  an  estimate  for  the  possible  displacement  of  the  origin: 
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(B.6) 


Ayo  = 


=  0.074  mm 
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B.3  Velocity  Measurement  Error 

There  were  two  main  sources  of  error  associated  with  the  velocity  measurements:  biasing 
error  and  random  error.  Biasing  error  could  occur  from  poor  calibration  of  the  optics  and 
misalignment  of  the  lasers  in  the  z-direction.  Random  errors  which  might  have  influenced  the 
velocity  measurements  included  the  processing  performed  by  the  BSA’s,  fluctuations  in  the  wind 
tunnel  conditions,  traverse  positioning,  and  laser  angle  alignment. 

In  the  following  analysis,  the  symbol  u  was  used  to  represent  any  of  the  velocity-related 
terms  (m  ,  u' ,  v  ,  v' );  when  there  was  a  discrepancy  between  the  formulations  for  the  individual 
components  it  was  noted  accordingly. 

B.  3. 1  Biasing  Errors.  The  first  type  of  possible  biasing  error  resulted  from  imperfections 
in  the  optics:  slight  changes  in  the  laser  intersection  angle  due  to  possible  defects  in  the  optical 
glass,  changes  in  the  interference  pattern  due  to  the  fact  that  the  lasers  were  not  going  through  the 
optical  glass  normal  to  the  surface  (the  3.5°  offset  relative  to  the  z-coordinate  direction  caused  each 
beam  to  be  offset  at  a  slightly  different  angle),  etc.  Imperfections  such  as  these  would  cause  a  shift 
in  the  calibration  factor,  ,  used  to  convert  the  Doppler  frequencies  into  velocities.  However,  if 

the  measured  velocities  were  scaled  by  other  measured  velocities  (i.e.,  u/u^ ,  u'^Ju,  etc.),  the 
error  in  the  calibration  factor  would  be  scaled  out.  In  the  cases  where  the  measured  velocities  were 
not  scaled  by  other  measured  velocities  (-pw'v'/r,,, ,  etc.)  it  was  assumed  that  this  bias  was 
negligible  compared  to  other  sources  of  error. 

The  second  source  of  biasing  error  was  due  to  the  fact  that  the  lasers  were  aligned  at  a  3.5° 
angle  to  the  mean  flow.  This  alignment  caused  the  lasers  to  consistently  measure  the  w-velocity 
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component  as  being  lower  than  the  actual  value.  The  actual  value  of  m  ,  u^,  was  given  by  = 


My'cos^3.5°j  =  1.002m  .  Therefore,  the  biasing  error  due  to  the  laser  alignment  was  approximately 
0.2%  and  was  ignored  when  compared  to  the  other  sources  of  error. 

B.3.2  Random  Errors.  As  mentioned  above,  there  were  four  main  possibilities  for 
sources  of  random  error  in  the  velocity  measurements:  error  associated  in  the  processing 
capabilities  of  the  BSA’s,  error  associated  with  random  fluctuations  in  the  wind  tunnel  conditions, 
traverse  location  errors,  and  errors  from  the  streamwise  angular  alignment  of  the  lasers. 


B.  3.2.1  Error  from  BSA  Calculations.  The  error  incurred  from  the  BSA’s 
calculation  of  the  velocity  was  due  to  the  limitations  on  the  curve  fit  used  to  set  up  the  FFT  which 
was  used  to  calculate  the  Doppler  frequency,  .  In  general,  the  frequency  resolution  of  the  BSAs 
(i.e.,  the  accuracy  of  the  FFT)  was  given  by  [10] 


A/o  = 


n(BW) 

2*iV 


(B.7) 


For  the  current  experiment,  different  bandwidth  settings  were  used  on  each  BSA  as  given  in  Table 
5.1.  For  the  bandwidth  chosen  on  BSA  1,  n  =  12  ;  for  BSA  2  «  =  1.5 .  With  these  values,  the 
selected  bandwidths,  and  a  record  length  of  16,  the  errors  in  the  u  and  v  components  of  velocity 
were 


Aug^  -  0.095  m  /  s 
=  0.090  m  /  s 


(B.8) 


B.  3.2.2  Repeatability.  The  random  errors  associated  with  the  wind  tunnel 
conditions  were  determined  from  the  repeatability  test  results  (see  Appendix  C).  The  twenty  test 
runs  gave  these  errors  as  «0.19%,  «0.04%,  «0.07%,  £•„,  «0.03%  and 

^  ^ep  ^rep  ‘*rep  ^rep 

e-Tj  «  0.02% ,  where  and  s^,  were  defined  as  lu  and  =  Av' „  lu  . 

uWrtp  ’  ^np  ^rep  ^rep  /  ^rep  I 


B-5 


B.3.2.3  Traverse  Location  Error.  The  error  to  the  velocity  measurements  due  to 
the  traverse  location  error  depended  on  the  position  of  the  control  volume  within  the  boundary 
layer.  In  regions  where  the  velocity  gradients  were  high,  a  small  change  in  the  traverse  location 
caused  a  large  change  in  the  measured  velocity.  For  two-dimensional  boundary  layer  flow,  the 
error  in  the  velocity  depended  primarily  on  the  velocity  gradient  normal  to  the  wall: 

AV;,  AVj  (B.9) 

,  ^y  I 

where  the  subscript  I  once  again  denoted  the  data  collection  point  away  from  the  wall.  It  should 
be  noted  that,  as  I  increased,  the  traverse  error  increased  while  the  velocity  gradient  decreased, 
helping  to  minimize  the  error. 

B.3.2.4  Streamwise  Angular  Alignment  Error.  The  final  source  of  error  in  the 
velocity  measurements  was  due  to  the  angular  alignment  of  the  lasers.  Assuming  the  lasers  could 
be  misaligned  by  the  angle  A^,  the  relationship  between  the  actual  value  of  u,  u^,  and  the 
measured  value  of  u  is  shown  in  Figure  B.2  below. 


V. 


Figure  B.2:  Velocity  Error  Due  to  Angular  Misalignment 


Examination  of  Figure  B.2  showed  the  relationship  between  the  measured  and  actual 


velocity  components  was  given  by 


=  mcos(A^)  -  vsin(A^) 
=  Msin(A^)  +  vcos(A^) 


(B.IO) 
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Using  Taylor  expansions  for  the  trigonometric  terms,  the  associated  errors  were  given  as 

(B.ll) 

where  was  in  radians.  Normalizing  these  equations  by  w ,  the  percent  changes  were  found  to 


Am^  =u^-u  =  -  v(A^) 

-  V  =  u{A^)  -  v{A^f 


be 


(B.12) 


/“I  =  ^ 

Note  that  the  last  term  in  could  be  ignored  because  it  was  third-order  ( vju  «  1 .0). 


The  total  error  in  the  velocity  was  then  given  by  the  Euclidian  norm  of  the  independent 
error  terms: 


+  sl  +£l  +el 

^rep 


(B.13) 

The  error  in  the  velocity  components  was  calculated  at  each  point  across  the  boundary 
layer  and  was  found  to  be  maximum  near  the  wall  due  to  the  large  velocity  gradients.  In  order  to 
simplify  the  remaining  error  calculations,  the  error  in  the  velocity  was  taken  to  be  constant  across 
the  boundary  layer,  equal  to  the  maximum  error  possible.  In  this  case. 


£u=^  =  8-8%,  ^  =  8.5% 

u  u 

Aw  *  Av  ^ 

^  =  1.7%,  s,.=^  =  7.9% 
u  u 


(B.14) 


With  the  error  in  the  measured  velocities  known,  the  error  in  the  total  velocity  could  be 
computed  from  Eq.  (6.3)  and  was  given  by 


'  ("v’ 

+  [ir^v 
\  u 


(B.15) 


^£^=8.8% 


The  error  in  the  second-order  velocity  correlation,  u'v' ,  was 
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(B.16) 


Finally,  note  that  the  slope  of  the  ?/-velocity  profile  went  to  zero  at  the  fi-eestream;  this  led 
to  a  reduction  in  the  error  associated  with  the  edge  velocity.  The  error  was  only  dependent  on  the 
BSA’s  processing  capabilities,  the  repeatability  and  the  angular  alignment: 


=M  +^l  +^1  =1-3%  (B.17) 

Assuming  the  v-component  of  velocity  was  measured  just  as  accurately  in  the  fi'eestream  (v  did  not 
go  to  zero,  but  the  slope  was  reduced  and  the  error  should  have  been  less  than  that  given  by  Eq. 
(B.  15)),  the  error  in  the  total  velocity  at  the  edge  was 


=  13% 


(B.18) 


B.  4  Total  Pressure  and  Temperature 

Stagnation  values  of  the  pressure  and  temperature  were  also  measured.  The  error  in  the 
stagnation  pressure  was  estimated  at  1380  Pa  out  of  approximately  2.13  x  10^  Pa,  or 


Po  2.13  X 10^ 


(B.19) 


There  were  two  source  of  error  in  the  total  temperature:  the  error  in  the  measurement  of 
the  total  temperature  and  the  assumption  that  the  total  temperature  was  constant  across  the 
boundary  layer.  Recall  from  Eq.  (6.25)  that  the  assumption  of  constant  total  temperature  could 
lead  to  an  error  of  up  to  4%.  The  actual  temperature  measurement  was  estimated  to  have  an 
accuracy  of  ±02  K ;  with  total  temperatures  of  at  least  292  K,  this  gave  an  error  of  =  0.1%  . 

Combining  these  two  sources  of  error  gave  the  maximum  error  in  the  total  temperature  as 

=^  =  4.0% 

20  'T 

■*0,00 
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B.5  Mach  Number,  Temperature,  Pressure,  and  Density  Error 


The  error  in  the  Mach  number  was 


(^e)  +(2^0 


^9.0% 


(B.21) 


and  the  error  in  the  edge  Mach  number,  ,  was 


= 


(^s) 


2  (I 
^  I  ^'2 


=  1.5% 


(B.22) 


Using  the  error  in  ,  the  error  in  the  edge  temperature  and  pressure  were  given  by 


St  — 


+ 


(j'-i) 


-1 


=  1.8% 


(B.23) 


Pt 


Y 


-1 


MUm, 


=  6.4% 


(B.24) 


and  the  error  in  the  edge  density  was 

=6.6% 

The  errors  in  the  wall  pressure  and  temperature  were  from  direct  measurements  and  were  thus 
similar  to  the  error  in  the  stagnation  measurements:  Sj,^  =  0.6%  and  =  0.1% .  Using  these 


results,  the  error  in  the  wall  density  was  found  to  be 

and  the  error  in  the  wall  viscosity  (which  was  calculated  using  Sutherland’s  formula)  was 


s 


(3/2)(r,  +  iio.4)-7;^ 
(7; +  110.4) 


0.1% 


(B.27) 
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If  the  density  and  temperature  were  calculated  using  the  Crocco-Busemann  integral 
relation,  the  error  was  dependent  on  the  validity  of  the  assumptions  used  in  the  derivation  of  the 
relationship.  However,  the  relation  has  been  found  to  be  very  accurate  when  the  assumptions  hold 
[6],  and  the  error  was  assumed  to  be  less  than  the  error  produced  in  the  FPG  temperature  and 
density  calculations.  In  order  to  simplify  the  analysis  of  any  variables  which  follow,  the  errors  in 
the  temperature  and  density  were  taken  as  those  incorporated  from  the  assumption  of  a  linear 
pressure  distribution  through  the  FPG  boundary  layer. 

When  the  temperature  and  density  were  calculated  assuming  a  linear  pressure  distribution 
through  the  boundary  layer,  the  comparison  of  the  computed  density  to  CFD  and  experimental 
results  in  Figures  7.8  and  7.24  showed  that  the  assumptions  used  in  their  calculation  were  fairly 
accurate,  and  the  error  in  pressure  across  the  boundary  layer  could  be  estimated  as 
Sj,  »  Sp^  =  6.6% .  The  error  in  the  temperature  was  then  given  by 


S  'T  — 


inf 


+ 


(r  - 1) 


i+LAm'^ 


(B.28) 


<11.7% 


and  the  error  in  density  was 

Sp=^sl+el=UA% 


(B.29) 


B.  6  Boundary  Layer  Thickness  Error 

The  error  in  the  boundary  layer  thickness  was  determined  graphically.  The  process  used  to 
generate  ES  is  shown  in  Figure  B.3.  The  velocity  or  Mach  number  profile  (scaled  by  the  edge 
value)  was  plotted  versus  the  location  off  the  wall;  error  bars  were  included  to  give  the  possible 
values  of  the  velocity  or  Mach  number.  The  line  drawn  on  the  right  side  of  the  graph  indicated  the 
position  of  Q995u^ ,  the  value  used  to  determine  the  edge  condition.  The  two  horizontal  lines  gave 
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Figure  B.3:  Graphical  Derivation  of  AS 

the  position  of  the  minimum  and  maximum  value  at  which  the  non-dimensional  velocity  profile 
could  reach  0.995;  the  value  of  S  fell  somewhere  between  these  two  lines.  It  is  obvious  from 
Figure  B.3  that  the  error  in  S  was  relatively  large;  the  percent  errors  were  estimated  as 
£g  =  28%  for  both  the  velocity  and  Mach  number  boundary-layer  thicknesses. 

B.  7  Wall  Shear  Stress  Error 

The  error  in  the  wall  shear  stress  due  to  van  Driest’s  relationship  has  been  documented  as 
« 10%  [9];  however,  this  was  the  error  when  the  assumptions  used  to  generate  the  wall  shear 

stress  are  assumed  valid.  The  analysis  in  Chapter  7  showed  that  the  wall  shear  stress  was  not 
predicted  well  by  van  Driest’s  relationship  due  to  inaccuracies  in  the  assumptions  used  to  derive 
the  relationship.  As  a  result,  the  error  in  the  wall  shear  stress  was  estimated  by  comparing  the 
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computed  value  to  the  value  which  would  shift  the  data  points  in  Figure  7.16  on  top  of  the 
empirical  correlation.  The  calculated  value  was  70  Pa  and  the  value  needed  to  shift  the  data  was 
approximately  75  Pa,  an  error  of  7.1%. 


B.  8  Van  Driest  Velocity  Error 

The  error  in  the  van  Driest  velocities  was  due  to  the  error  in  the  effective  velocity  and  the 


error  in  the  scaling  parameter  u  ,  which  were  found  to  be 


=iois^}\(05sJ\(2sy^lS.9% 

=3.5% 

respectively.  Thus,  the  error  in  the  van  Driest  velocity  was 

=>«% 


(B.30) 

(B.31) 


(B.32) 


B.P  Turbulence  Statistics  Error 

The  error  to  the  second-order  velocity  correlation  u'v'  was  given  by 

^i;v=A/W+(^v')'=8.1%  (B.33) 

Using  the  variable  x  to  denote  either  the  «  or  v  component  of  velocity,  the  error  in  the  root- 
mean  square  was  given  by 


which  led  to 


1 


Act^=- 

\ZVf] 

2 

Yx'Ar'’ 

*  2 

n 

n 

^2 


S(-') 


n-ir 


A2 


S(-') 


=  s. 


(B.34) 


(B.35) 
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With  this  formulation,  the  error  in  the  RMS  values  of  the  u  and  v  components  of  velocity  were 


given  by 


£,,=7.9% 

Using  a  simitar  derivation,  the  error  in  the  skewness  and  flatness  were  found  to  be 


(B.36) 


£sk,  =V18£^, 

+(4^x'f  = 


(B.37) 


resulting  in 

=  12%  =  33.5%  (B.38) 

Sp,^=9.6%  £pi=U.l% 

The  error  in  the  intermittency  function  was  the  same  as  the  error  to  the  flatness, 

£  =Sp,  =9.6%  (B.39) 

and  the  error  in  the  correlation  coefficient  was  given  as 


(B.40) 


B.IO  Strain  Rate  Error 

The  error  in  the  strain  rate  was  dependent  on  the  method  of  calculation;  for  the  second- 
order  central  differences,  the  error  in  the  derivative  of  a  general  variable  u  with  respect  to  y  (or  z) 
was 

The  percent  error  in  the  position  was  dependent  on  whether  the  y-  or  z-derivative  was  being 
examined.  The  error  in  the  z  position  was  essentially  negligible,  while  the  error  in  the  y-direction 
increased  across  the  boundary  layer  as  given  in  Eq.  (B.4).  The  maximum  error  in  the  y  position 
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was  at  the  outer  edge  of  the  boundary  layer,  where  j  »  10  mm,  leading  to  «  1% .  Thus,  the 
errors  in  the  central-difference  strain  rates  were 


~  17-7%  ^SuIAz  —  12.4% 

—  17.1%  =  12.0% 

The  error  in  the  x-direction  derivatives  were  computed  from  Eq.  (6.54)  assuming  the 

position  error  in  the  x  direction  was  negligible.  Thus, 

=  =  15.2% 

^a./^  =  V3£„  =  14.7% 


B.ll  Error  Summary 

The  errors  calculated  above  are  summarized  below  in  Table  B.  1 .  For  the  sake  of  brevity, 
the  error  associated  with  each  figure  in  Chapter  7  has  not  been  presented;  it  should  suffice  to  say 
that  if  one  variable  is  scaled  by  another,  the  error  in  the  resulting  variable  is  the  Euclidean  norm 
(cf.  Eq.  (B.3))  of  the  constituent  parts. 
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C.  Compliance  With  Criteria  of  Settles  and  Dodson 

In  order  to  ensure  compliance  with  the  criteria  established  by  Settles  and  Dodson  [31] 
which  were  presented  in  Chapter  1,  a  list  of  the  criteria  and  how  they  were  met  is  presented  in 
Table  C.l.  If  it  was  found  that  one  or  more  of  the  criteria  still  needed  to  be  met,  they  were 
addressed  individually  in  the  remaining  sections  of  this  chapter. 


Table  C.l;  Compliance  With  Criteria  of  Ref  [3 1] 


Criterion 

Number 

Criterion 

Category 

HowAVhere 

Criterion  Met 

1 

Baseline  Applicability 

1 

M,  =  2.9  «  3.0 

2 

Simplicity 

Geometry  able  to  be  modeled  by  CFD 
(c;:  Ref  [17]) 

3 

Specific  Applicability 

Has  been  compared  to  CFD  in  order  to  validate 
code  {cf.  Ref  [17]) 

4 

Well-Defined  Experimental 

Incoming  conditions  tabulated  (Section  8.1) 

Boundary  Conditions 

2D  Check;  see  Section  C.2 

5 

Defined  Error  Bounds 

Appendix  B;  Section  C.l 

6 

Documentation  of  Data 

Section  C.3 

7 

Spatial  Resolution  of  Data 

Figures  in  Chapter  7  are  large  with  as  much 
resolution  as  possible. 

C.  1  Repeatability 

The  repeatability  of  the  test  results  was  determined  using  two  separate  methods.  The  first 
was  to  take  twenty  test  runs  at  a  given  location  to  find  the  effect  of  wind  tunnel  fluctuations  on  the 
measurements.  These  tests  provided  confidence  intervals  for  the  mean  flow  variables,  turbulence 
intensities,  skewness,  flatness  and  the  second-order  velocity  correlation.  A  confidence  interval  is  a 
measure  of  the  probability  that  the  error  in  the  measured  mean,  x  ,  will  be  less  than  a  certain  value 


C-1 


when  compared  to  the  actual  value,  That  is,  a  99%  confidence  interval  of  1%  means  that 


there  is  a  99%  probability  that 


<0.01 


I  I 

However,  the  value  of  was  not  known,  so  it  was  not  possible  to  directly  calculate  the  value  of 

the  confidence  interval  on  the  RHS  of  the  above  equation.  In  order  to  obtain  a  value  for  the 

confidence  interval,  it  was  necessary  to  use  Student’s  t  distribution  and  the  relation  [3] 

^  st„.ai2  (C.l) 

X  ~  x4n 

where  N  was  the  number  of  data  points  used  to  calculate  x  .  UN  was  large  enough,  s  = 
Student’s  t  distribution  was  a  tabulated  function  based  on  «  =  -  1  and  a  .  a  was  a  variable 

which  described  the  accuracy  of  the  confidence  interval  and  was  given  by  a  =  1- 1^/100% ,  where 
was  the  percent  accuracy  of  the  confidence  interval.  Therefore,  a  99%  confidence  interval 


meant  of  =  0.01. 

For  the  repeatability  study  in  the  present  experiment,  twenty  test  runs  were  performed  but 
two  were  discarded  due  to  the  fact  that  the  mean  velocity  was  extremely  low  when  compared  to  the 
other  test  runs,  suggesting  the  data  collection  was  initiated  before  the  tunnel  reached  steady  state  or 
was  not  cut  off  until  the  vacuum  was  depleted  and  the  flow  velocity  decreased.  In  any  event,  this 
same  practice  of  manually  evaluating  the  mean  flow  velocities  was  used  in  the  actual  data 
collection  tests,  so  discarding  these  two  tests  should  not  have  biased  the  confidence  intervals 
associated  with  the  presented  results.  After  discarding  the  two  flawed  data  sets,  there  were  18 
samples  used  to  calculate  the  99%  confidence  intervals;  from  Ref  [3],  =^170 os  =2.898. 

Using  this  value  in  Eq.  (C.  1),  the  99%  confidence  intervals  for  the  velocity  data  and  the  turbulence 
statistics  output  by  the  Burstware®  software  were  computed  and  are  presented  in  Table  C.2. 
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Table  C.2:  Confidence  Intervals 


u 

V 

U'rms 

u'v' 

0.15% 

0.076% 

1.53% 

0.99% 

13.59% 

SK 

SK 

Fh 

FK 

2.60% 

14.10% 

3.23% 

8.00% 

The  second  method  of  proving  the  repeatability  was  to  take  several  traverses  across  the 
boundary  layer  at  each  test  location.  The  test  results  for  four  separate  traverses  across  the  flat 
plate  boundary  layer  are  presented  in  Figure  C.l;  the  most  sensitive  variable,  u'v' ,  showed 
excellent  repeatability. 


1.6 

1.4 

1.2 

1.0 
t  0  8 
0.6 

0.4 

0.2 

0.0 

0.0000  0.0005  0.0010  0.0015  0.0020  0.0025 

-uV /  u^ 

Figure  C.l:  Repeatability  of  Second-Order  Correlation  Over  a  Flat  Plate 
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C.2  Two-Dimensionality  Check 

In  accordance  with  the  guidelines  set  up  by  Settles  and  Dodson,  a  two-dimensionality 
check  was  performed.  The  check  was  accomplished  by  comparing  data  taken  when  building  the 
finite-difference  grid.  Three  profiles  were  taken  along  the  z-axis;  the  velocity  profiles  generated 
are  compared  in  Figure  C.2.  Although  the  values  at  z  =  -5  mm  were  slightly  lower  than  those  at 
z  =  5  mm,  the  difference  was  essentially  negligible  and  the  flow  was  assumed  to  be  two- 
dimensional. 


Figure  C.2:  Velocity  Comparison  for  Two-Dimensionality  Check 


C.3  Tabulated  Data 

Tables  C.3  to  C.6  on  the  following  pages  contain  point-by-point  data  collected  by  the  LDV 
system  and  the  data  from  the  resulting  calculations. 
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Table  C.3:  Flat  Plate  16-Sample  Raw  Data 
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Table  C.3:  Flat  Plate  16-Sample  Raw  Data  (con't) 


Cl, 

o 

VO 

O 

VO 

d 

VO 

d 

VO 

d 

cn 

VO 

d 

VO 

d 

«n 

VO 

d 

VO 

d 

VO 

d 

VO 

d 

in 

VO 

d 

VO 

d 

r- 

so 

d 

VO 

d 

VO 

d 

VO 

d 

d 

r- 

d 

VO 

d 

CO 

r-* 

d 

d 

in 

d 

00 

d 

00 

d 

m 

00 

d 

m 

00 

d 

00 

d 

00 

d 

ON 

d 

cs 

ON 

d 

i960 

VO 

ON 

d 

ON 

d 

H 

Ov 

vq 

fn 

a^ 

00 

r— * 

m 

vd 

t—H 

o> 

5 

VO 

vq 

00 

vq 

id 

Os 

00 

fo 

d 

Os 

r- 

I-H 

vd 

os 

in 

00 

Os 

0\ 

Os 

o 

vd 

o 

<N 

<N 

vq 

os 

m 

fn 

in 

in 

O'! 

vq 

00 

m 

r4 

CNJ 

vd 

m 

<N 

00 

Os 

VO 

<N 

VO 

00 

Os 

00 

<N 

o 

cq 

*n 

Os 

<N 

00 

<N 

On 

cs 

(N 

O 

00 

VO 

CO 

r- 

p 

in 

p 

On 

CO 

CO 

<N 

in 

m 

p 

in 

00 

m 

00 

cs 

cs 

VO 

00 

cd 

CO 

VO 

cs 

VO 

VO 

in 

p 

cs 

CO 

VO 

VO 

ON 

o 

00 

VO 

od 

00 

r- 

00 

vd 

m 

ON 

o 

Os 

cs 

o 

in 

o 

I-H 

* 

t 

5? 

<s 

00 

00 

P"H 

00 

00 

ON 

VO 

O) 

00 

r- 

Os 

Os 

rs 

os 

os 

00 

00 

<N 

d 

<N 

OS 

Os 

r-H 

Os 

in 

Os 

o 

in 

Os 

•n 

o 

d 

<N 

<N 

00 

d 

m 

1 '  1 

d 

<N 

O 

•n 

d 

m 

cq 

CS 

<N 

Os 

d 

<N 

VO 

d 

(N 

m 

tq 

(N 

o 

VO 

CN 

<S 

o 

CN 

04 

(S 

in 

cd 

<s 

VO 

rq 

cd 

cs 

in 

cs 

cs 

p 

CS 

CS 

cs 

r^ 

CS 

00 

p 

in 

cs 

CO 

in 

cs 

ON 

p 

in 

cs 

ON 

Os 

in 

cs 

ON 

rq 

in 

cs 

A«  (m/s) 

T— H 

VO 

vd 

s 

ON 

o 

fd 

<N 

Os 

Os 

<N 

<N 

00 

<N 

00 

r-‘ 

<N 

r- 

Os 

r- 

00 

Os 

00 

Os 

00 

00 

r-H 

Os 

<S 

rq 

00 

<N 

00 

00 

CS 

Os 

00 

p 

od 

r-* 

r- 

CO 

od 

r- 

On 

00 

On 

K 

S 

CO 

0- 

s 

od 

VO 

CO 

K 

cs 

o 

S 

in 

Os 

vd 

On 

vd 

VO 

VO 

vd 

00 

tT 

vd 

CO 

4) 

1 

fo 

o 

fq 

<d 

1 

I 

r-j 

O 

W 

00 

ON 

<N 

ro 

O 

W 

o 

rq 

m 

o 

W 

Os 

<s 

m 

o 

W 

<N 

d 

1 

I 

s 

s 

1 

i 

I 

1 

I 

9 

w 

m 

00 

<d 

CO 

o 

1 

W 

in 

Os 

CO 

CO 

o 

1 

w 

CS 

00 

cd 

CO 

o 

pq 

<N 

<N 

CO 

o 

W 

CO 

00 

cd 

CO 

o 

1 

W 

00 

p 

cd 

I 

I 

CO 

p 

W 

CO 

0- 

cs 

s 

K 

i 

I 

I 

CO 

o 

1 

s 

o 

cs 

CO 

o 

1 

pq 

m 

p 

CO 

9 

p 

CO 

o 

1 

pq 

o 

1 

s 

w 

CO 

r' 

vd 

I 

ft; 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

O 

O 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

O 

O 

d 

O 

o 

d 

o 

o 

d 

O 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

O 

O 

d 

o 

o 

d 

o 

o 

d 

O 

O 

d 

O 

o 

d 

O 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

> 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

O 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

o 

o 

d 

> 

1 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

s 

o 

o 

d 

o 

1 

o 

d 

o 

o 

pij 

o 

o 

d 

o 

o 

o 

o 

d 

o 

1 

o 

o 

d 

o 

o 

ih 

o 

o 

d 

o 

o 

i 

o 

d 

o 

o 

i 

o 

d 

o 

o 

i 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

o 

d 

o 

o 

li 

o 

o 

d 

o 

s 

o 

p 

d 

o 

o 

pSj 

o 

o 

d 

o 

o 

+ 

pq 

o 

o 

d 

o 

o 

tS 

o 

o 

d 

o 

o 

o 

d 

o 

o 

li 

o 

o 

d 

o 

? 

s 

o 

d 

o 

o 

s 

o 

o 

d 

o 

o 

li 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

d 

o 

o 

+ 

g 

o 

d 

o 

o 

§ 

o 

d 

o 

o 

S 

o 

o 

d 

» 

t 

o 

1 

o 

o 

d 

o 

o 

^1 

o 

o 

d 

o 

o 

£ 

o 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

ui 

o 

o 

d 

o 

o 

s 

o 

o 

d 

o 

o 

4" 

W 

o 

o 

d 

o 

o 

^1 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

& 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

£ 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

+ 

pq 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

+ 

pq 

o 

o 

d 

o 

o 

pEj 

o 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

S 

o 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

X 

VO 

o 

d 

00 

o 

d 

C\ 

o 

d 

o 

r-H 

d 

<N 
«— 1 

d 

<N 

»— 1 

d 

d 

m 

d 

»-**4 

d 

in 

d 

d 

VO 

d 

VO 

d 

VO 

d 

00 

d 

00 

d 

o 

CN 

d 

CO 

CN 

d 

1-^ 

<N 

d 

»n 

CN 

d 

ON 

CS 

d 

CO 

d 

in 

CO 

d 

d 

CO 

d 

On 

d 

d 

cs 

in 

d 

in 

d 

cs 

VO 

d 

VO 

d 

cs 

d 

00 

d 

X 

o 

d 

ov 

o 

d 

o 
1— ^ 

d 

»— H 

d 

fO 

r*H 

d 

m 

r-H 

d 

*n 

d 

in 

I— ^ 

d 

«n 

d 

d 

in 

d 

o 

d 

00 

d 

00 

d 

o 

<N 

d 

o 

<N 

d 

cs 

<s 

d 

«n 

CN 

d 

CO 

<N 

d 

r-* 

<N 

d 

cs 

CO 

d 

CO 

d 

00 

CO 

d 

VO 

d 

00 

d 

CO 

m 

d 

ON 

d 

in 

d 

o 

VO 

d 

00 

VO 

d 

1— H 

r-' 

d 

o 

00 

d 

00 

00 

d 

VO 

d 

00 

r-* 

d 

00 

r- 

d 

Ov 

ov 

d 

o 

<N 

o 

<N 

o 

<N 

I-H 

o 

m 

1— 1 

cq 

<q 

o 

o 

in 

o 

vq 

CO 

vq 
»— • 

«n 

o 

00 

00 

Os 

00 

Os 

r-H 

o 

p 

<N 

o 

m 

<N 

00 

cs 

o 

p 

cd 

o 

in 

cd 

o 

p 

in 

cs 

in 

cs 

o 

p 

o 

p 

»n 

o 

p 

in 

o 

o 

vd 

o 

in 

vd 

o 

o 

ON 

o 

C-6 


Table  C.3:  Flat  Plate  16-Sample  Raw  Data  (con't) 
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Table  C.4:  Flat  Plate  16-Sample  Coincident  Data 
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Table  C.7:  Flat  Plate  8-Sample  Raw  Data  (con't) 
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Table  C.8:  Flat  Plate  8-Sample  Coincident  Data  (can't) 
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Table  C.9:  FPG  16-Sample  Coincident  Data  (x  =  70.98  cm) 


C-25 


Table  C.  10:  FPG  16-Sample  Coincident  Data  (x  =  70.44  cm) 
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